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Влияние динамических свойств взаимодействующих подсистем на эволюцию 
формирования избирательного переноса в узлах трения' 
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Введение. Типичным примером процессов самоорганиза- 
ции в трибосистемах является избирательный перенос. В 
этом случае в области контакта формируются сопрягаю- 
щие поверхности сервовитной пленки, что принципиально 
меняет условия трения и изнашивания. Для образования 
избирательного переноса в зоне сопряжения контактиру- 
ющих поверхностей необходима некоторая мощность 
необратимых преобразований подводимой энергии, кото- 
рая зависит от упруго-диссипативных свойств контакти- 
рующих поверхностей. 

Материалы и методы. Приводится математическая мо- 
дель динамической системы с учетом эволюционно изме- 
няющейся сервовитной пленки. Ее формирование зависит 
от фазовой траектории мощности необратимых преобра- 
зований в зоне сопряжения поверхностей и представляет- 
ся в виде интегрального оператора Вольтерры второго 
рода. 

Результаты исследования. Приводится анализ результа- 
тов исследования, в том числе зависимости формирования 
сервовитной пленки от динамических параметров взаимо- 
действующих подсистем. Впервые рассматривается про- 
блема математического моделирования эволюции свойств 
узла трения с образованием или разрушением сервовитной 
пленки. 

Обсуждение и заключение. В ходе эволюционного про- 
цесса образования сервовитной пленки меняются пара- 
метры динамической связи, формируемой в узле трения. 
Следовательно, меняются и динамические свойства си- 
стемы. Впервые проанализирована динамика трибосисте- 
мы в процессе образования избирательного переноса. 


Ключевые слова: динамическая система трения, избира- 
тельный перенос, сервовитная пленка, эволюция. 


Образец для цитирования: Заковоротный, В. Л. Влияние 
динамических свойств взаимодействующих подсистем на 
эволюцию формирования избирательного переноса в уз- 
лах трения/ В. Л. Заковоротный, В. Е. Гвинджилия, 
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Введение. После опубликования работ И. Пригожина [1-3] и Г. Хакена [4, 5] многие вопросы функци- 
онирования технических систем, взаимодействующих с различными средами, рассматриваются с точки зрения 
их самоорганизации [6-12]. Если следовать синергетической парадигме анализа, то при изучении таких систем 
прежде всего следует выполнить процедуру расширения размерности пространства состояния [7, 8]. Это озна- 
чает, что необходимо дополнительно рассмотреть модель среды в координатах состояния системы. Имеет место 
обмен со средой информацией, материалом, энергией и пр. Здесь возможны когерентные взаимодействия раз- 
личной физической природы, то есть явления синергетики. Типичным примером самоорганизации является 
эффект избирательного переноса, открытый И. В. Крагельским и Д. Н. Гаркуновым в 1956 году [13]. Идеи са- 
моорганизации были сформулированы в 60-е годы ХХ века Б. И. Костецким [14] и развиты в работах по струк- 
турной приспосабливаемости совместно с Л. И. Бершадским [15] и Н. А. Буше [16]. В ряде работ исследуется 
формирование сервовитной пленки в трибосистеме при ее выходе на режим избирательного переноса [17—22]. 
Показано, что для формирования избирательного переноса необходимы определенные трибохимические реак- 
ции, а также некоторая траектория мощности необратимых преобразований в области сопряжения «система — 
среда» [17—24]. При выходе на режим избирательного переноса наблюдаются автоколебания. В некоторых слу- 
чаях формируются хаотические притягивающие множества деформационных смещений контактирующих 
пар [17, 22—24]. 

Материалы и методы. Таким образом, при изучении перехода на режим избирательного переноса 
необходимо рассматривать динамику трибосистемы в единстве подсистем, взаимодействующих через трение, и 
динамической связи, формируемой узлом трения. Параметры такой динамической связи и формирование сер- 
вовитной пленки зависят от траектории мощности необратимых преобразований по совершенной работе. В ста- 
тье рассматривается математическое моделирование такой эволюционной системы и анализируется зависи- 
мость эволюции от динамических параметров взаимодействующих подсистем. 

Результаты исследования 

Математическое моделирование системы. При математическом моделировании воспользуемся ранее 
полученными результатами. В [17, 23, 24—26] показано, что основные динамические свойства системы трения 
можно раскрыть на основе следующих допущений: 

— образец является абсолютно жестким; 
— деформации индентора рассматриваются в плоскости, нормальной к контактирующей поверхности и прохо- 
дящей через направление скорости относительного скольжения. 

При изучении динамики можно ограничиться первыми формами колебаний. Тогда модель системы 
представляется уравнением 





о 
ах АХ 
т +й— +сХ = Ё; (1), (1) 
а2 Ге 
_|т 0 п № сл С21 
где т = В = 5. С = —ЩЩ Положительно определенные, симметричные матрицы 
О т й 2 22 С],2 С22 
инерционных, скоростных и упругих коэффициентов; Х ={Х1,Х>} — вектор деформационных смещений 


вершины индентора; 25 (0) =4Ё у (1), Е 5 [О — вектор сил, действующих на вершину индентора (рис. 1). 


Трибосреда 






Образец 
№ 
Рис. 1. Схема динамической системы трения 


Будем следовать синергетической концепции анализа, то есть представим силы Е; (1) в координатах 


состояния. Для этого введем понятие трибосреда — это третье тело, которое формируется в переходной зоне 
между контактирующими поверхностями. При сближении меняются условия их взаимодействий, что физиче- 
ски проявляется в изменении таких параметров, как фактическая площадь контакта, производство тепла, коэф- 
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фициент трения, условия молекулярного взаимодействия, процессы диффузии и пр. [27—30]. Однако первичны- 
ми причинами всех изменений являются взаимодействия, обусловленные механикой контакта. В этом случае за 
счет изменения фактической площади контакта непропорционально быстро возрастают силы нормального дав- 
ления, которые препятствуют сближению поверхностей. Моделирование этих сил [23, 24] удобно рассматри- 
вать в виде следующей функции сближения 

Е (Ху) = Е.оехр[-о( 1 )]-Ц, | в (0, №), (2) 


где а — коэффициент крутизны нарастания силы контактного взаимодействия В [мм , И — внешняя сила. 


В названных выше работах показано, что при формировании сервовитной пленки в функции сближения 
образуется потенциальный барьер, который фактически определяет несущую способность сервовитной пленки 
узла трения. Формирование потенциального барьера удобно представить в виде 


2 
Ф/ (Л) = Ф дехр[-ои (1 — 10) 1, (3) 
где О —_ Параметр, характеризующий крутизну потенциального барьера В мм? | В Ато —_ координата потен- 


циального барьера; Фо — эволюционный параметр в [кг]. Таким образом, функция сближения есть сумма 
Ну(М)=Ф(Х)+Е(Х1). Моделируя тангенциальную составляющую силы, примем во внимание, что при 
формировании сервовитной пленки силы трения на порядок меньше. Поэтому 
1 : . 
О (1), при: М (А1,0; 


К: Нх(,при: А) Хто, 


В ()= 
где ЕО уж | 


Например, в системе трения «сталь — глицерин — латунь» после формирования сервовитной пленки 
коэффициент трения уменьшается более чем на порядок. Ранее показано, что силы контактного взаимодействия 
характеризуются запаздыванием тангенциальной составляющей силы (силы трения) по отношению к силам 
нормального давления. В этом случае вместо (4) имеет место соотношение 


| 
ЗЕ (@(-Т®), при: Хи 


К Ну((-Т`^), при: ХХ о, 
где то р т“) —_ Постоянные времени запаздывания, зависящие от свойств трибосреды, сформированной В 


контактной области, и скорости относительного скольжения. 


В (5) всегда выполняется условие Г Фут 2. Кроме того, Г и т) уменьшаются при возрастании 


скорости относительного скольжения. 
В дальнейшем принята гипотеза, согласно которой потенциальный барьер формируется в зависимости 
от траекторий работы И МОЩНОСТИ необратимых преобразований при движении индентора относительно образ- 


ца в направлении скорости относительного скольжения. Причем, на Ф, 0 влияет не только текущее, но и пред- 


шествующие значения мощности. Поэтому справедливо 





1 
Фо (М) =В[(-ЭМ(ОаЕ, (6) 
0 
1 
где № => (У +АХ. / 4%) — мощность необратимых преобразований; и(#-—&) = {Вехр[ т (1-&)]} — ядро инте- 
грального оператора, учитывающее влияние предшествующих значений мощности; Г — параметр в [с], ха- 
рактеризующий сохраняющееся во времени влияние мощности на потенциальный барьер; В — параметр в 


и 


Уравнения (1)-(6) характеризуют математическую модель динамической системы трения с учетом эво- 
люции ее свойств при образовании сервовитной пленки. 

Изменение динамических свойств в зависимости от параметров потенциального барьера и дина- 
мических свойств индентора. Вначале рассмотрим свойства замороженной системы в предположении, что 


параметры сервовитной пленки являются заданными и неизменными. Положим также, что ти =0и те =0. 
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Рассмотрим пример изменения свойств системы в зависимости от стадий эволюционного преобразова- 
ния (рис. 2). 
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Рис. 2. Преобразование проекций фазовых траекторий на плоскости Х! — АХ / 


На начальной стадии преобразований (часть рис. 2, помеченная цифрой 1) в системе формируется 
единственная точка равновесия, которой соответствует трение с традиционным коэффициентом в пределах 0,2— 
0,3. Причем эта точка обладает свойством глобального притяжения. Затем (часть рис. 2, помеченная цифрой 2) 
при формировании сервовитной пленки на начальной стадии образуются две точки равновесия (Х\ 1, Х! 2). Об- 
ласть притяжения точки Х11, ограничена седлообразной сепаратрисой. Вся остальная область характеризуется 
притяжением к точке Х12. Важно подчеркнуть, что в точке равновесия Х1› уже сформирована сервовитная 


пленка, но ее величина и свойства (определяются Ф 0 ) не позволяют обеспечить притяжение траекторий всего 


фазового пространства. Поэтому в зависимости от начальных условий или возмущений (например, флуктуаций 
или формируемых притягивающих множеств деформационных смещений в вариациях относительно точки рав- 
новесия) возможно формирование свойств узла трения, которые принципиально отличаются друг от друга. В 
дальнейшем (часть рис. 2, помеченная цифрой 3) область, ограниченная седлообразной сепаратрисой, уменьша- 
ется и исчезает (часть рис. 2, помеченная цифрой 4). В последнем случае толщина сервовитной пленки и ее не- 
сущая способность возрастают. В модели это показано как увеличение потенциального барьера, и вторая точка 
равновесия приобретает свойства глобального притяжения. В этом случае в системе формируется стабильный 
избирательный перенос. 

Иллюстрации характеризуют динамику системы в предположении, что свойства сервовитной пленки 
являются замороженными, то есть не эволюционируют. Однако и в этом случае области притяжения точек рав- 
новесия, соответствующих трению с сервовитной пленкой и без нее, принципиально зависят от динамических 
параметров индентора, прежде всего от матриц его жесткости с и диссипации Й (1). В реальных условиях име- 


ет место эволюция системы во времени. Она зависит от скорости относительного скольжения и ‚ СИЛЫ и › ИС- 
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ходного состояния поверхности контактирующих тел и пр. На эволюцию, моделируемую интегральным опера- 
тором (6), влияют многие физические факторы, интегральный учет которых определяется ядром оператора. 
Поэтому при варьировании внешних условий (например, скорости относительного скольжения) необходимо 
учитывать влияние скорости на параметры интегрального оператора, а также на траекторию мощности необра- 
тимых преобразований. Экспериментально показано, что при увеличении скорости параметр Т;,1=1,2 в (6) 
уменьшается. Кроме того, силы трения (а следовательно, и траектория мощности необратимых В 
ний) зависят от скорости относительного скольжения. Заметим, что именно время, необходимое для установле- 
ния стационарной сервовитной пленки, определяет один из показателей, позволяющих идентифицировать па- 
раметры ядра интегрального оператора. 

Не останавливаясь на деталях, приведем пример изменения фазовых траекторий деформационных 
смещений в направлении Х| при варьировании динамических параметров индентора и свойств трибосреды 


(рис. 3). 
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Рис. 3. Преобразование проекций фазовых траекторий на плоскость А! — @Х/ 
в зависимости от параметров системы и условий трения 


Прежде всего, рассмотрим трение с достаточно большой скоростью относительного скольжения 


И, =2,0 м/с. В этом случае при формировании сервовитной пленки уменьшается мощность необратимых пре- 


образований и разрушение не наблюдается, так как скорость относительного скольжения достаточно высока 
(рис. 3, а). На начальном этапе при различных условиях траектории асимптотически притягиваются к точке 
равновесия без сервовитной пленки (точка А на рис. 3, а). После установления потенциального барьера все тра- 
ектории эволюционируют к точке равновесия А‘ на сервовитной пленке. На рис. 3, б, в иг приведены эволю- 
ционные кривые с одним начальным значением, соответствующим точке В. При этом деформационные смеще- 
ния и соответствующие им свойства системы трения стягиваются к одной из двух точек равновесия. Однако 
эволюционные кривые отличаются в зависимости от свойств трибосреды и добротности колебательных конту- 
ров подсистемы индентора без трения. На рис. 3, б приведен пример эволюционной траектории при условиях а, 


но при скорости относительного скольжения й = 1.0 м/с. В этом случае мощности, выделяемой в зоне трения, 


недостаточно для поддержания сервовитной пленки, и в системе периодически наблюдается трение с сервовит- 
ной пленкой и без нее. В деформационных смещениях это вызывает эффект низкочастотных колебаний. Харак- 
терно, что времена формирования и разрушения сервовитной пленки существенно отличаются. На рис. 3 это 
отображается значительным уменьшением скорости возвращения вершины индентора из точки А‘’ в точку А. 

Добротность подсистемы индентора можно увеличить (рис. 3, в). В этом случае в системе образуются 
сложные колебательные смещения. Некоторые их составляющие характеризуют низкочастотный материальный 
обмен между смазкой и поверхностью трения, а также колебательные смещения на частотах осцилляторов ин- 
дентора. В направлении, ортогональном к скорости относительного скольжения, формируются циклические 
силы. Благодаря этому, а также за счет динамического взаимодействия возможно кратковременное образование 
сил, превышающих величину потенциального барьера. Ситуация еще более усложняется, если дополнительно 
учесть запаздывание вариаций тангенциальных составляющих сил по отношению к нормальным их составля- 
ющим (рис. 3, г). В зависимости от запаздывания формируются различные притягивающие множества дефор- 
мационных смещений, а следовательно, и силы контактного взаимодействия. Замечено и формирование в си- 
стеме хаотической динамики, вносящей большую неопределенность в закономерности образования сервовит- 
ной пленки. 

Обсуждение результатов. Известно, что выход трибосистемы на режим избирательного переноса 
определяется трибохимическими реакциями и материальным обменом между контактирующими поверхностя- 
ми и смазкой. Для выхода трибосистемы на режим избирательного переноса с образованием сервовитной плен- 
ки и ее поддержания в процессе трения необходима некоторая мощность необратимых преобразований подво- 
димой энергии механической системы [17, 21-24]. В ходе эволюции динамической системы трения несущая 


способность формируемой сервовитной пленки меняется, что моделируется величиной потенциального барьера 


Фо › Препятствующего непосредственному контакту поверхностей трения. Если силы контактного взаимодей- 
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ствия в направлении, нормальном к этой поверхности, превышают Ф 0; ТО наблюдается динамическая пере- 


стройка системы трения с частичным разрушением сервовитной пленки. Поэтому условия трения в режиме из- 
бирательного переноса зависят не только от трибохимических реакций и мощности необратимых преобразова- 
ний, но и от динамического режима. Здесь принципиальное значение имеют устойчивость точки равновесия и 
формируемые в ее окрестности различные притягивающие множества деформационных смещений (предельные 
циклы, инвариантные торы, хаотические аттракторы). Кроме того, динамическую систему возмущают, напри- 
мер, биения при вращении заготовки. К тому же за счет принципиального изменения условий трения при непо- 
средственном контакте и при контакте через сервовитную пленку меняется мощность необратимых преобразо- 
ваний в сопряжении поверхностей, что изменяет эволюцию и условия поддержания сформированной сервовит- 
ной пленки. 

Заключение. Итак, можно сделать следующее заключение: на стабильность избирательного переноса 
влияют как динамические параметры взаимодействующих через узел трения подсистем, так и возмущения, все- 
гда имеющие место в реальной машине. Поэтому при изучении избирательного переноса с целью обеспечения 
безизносности узла трения необходимо достичь идентичности динамических систем на машине трения и в ре- 
альной машине. Для этого можно воспользоваться известными приемами [27, 31]. 

Приведенные данные показывают, что на формирование сервовитной пленки в узле трибоконтакта 
влияют не только трибохимические характеристики контактирующих поверхностей, но и параметры взаимо- 
действующих через трибосреду подсистем со стороны индентора и образца. Поэтому в каждом конкретном 
случае существует ограниченная область параметрического пространства динамических систем, а также воз- 


мущений, в которых образование избирательного переноса является правилом. 
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Введение. Функциональные свойства алмазных порошков 
обусловлены масштабным структурным фактором, по- 
скольку он влияет на формирование структурно- 
чувствительных механических свойств — пределов проч- 
ности и текучести. Учитывая качественную корреляцию 
между пределом текучести и твердостью, можно прогно- 
зировать повышение твердости, в том числе высокотвер- 
дых материалов. 

Материалы и методы. Рассмотрены, систематизированы 
и представлены в виде таблицы физические характеристи- 
ки основных типов наполнителей, входящих в состав 
упрочнителей. В качестве связки использована оловяни- 
стая бронза М2-01 (20 мас. % олова, 80 % меди). В нее 
добавляли ультрадисперсный природный алмаз (УДПА, 
0,5—4 мас. %), а также порошки природного алмаза (фрак- 
ции 3/2 мкм, 7/5 мкм, —40 мкм), полученные при перера- 
ботке алмазов на предприятии ОАО «Сахадаймонд». 
Названные материалы изготавливались на дробильно- 
классификационном оборудовании и вибростолах. Стадии 
получения порошков фиксировались с помощью растро- 
вой электронной микроскопии. Для классификации ал- 
мазных порошков по зернистости применяли вибросита. 
Физические и механические характеристики изготовлен- 
ных образцов испытывали по стандартным методикам. 
Для взвешивания использовали лабораторные электрон- 
ные весы четвертого класса ВЛТЭ-500. Плотность опреде- 
ляли микрометром МК 0-25 мм по ГОСТ 6507-78. 
Результаты исследования. Через фактическую и теорети- 
ческую плотности рассчитана пористость. Выяснилось, 
что с уменьшением размера фракции наполнителя наблю- 
дается улучшение физико-механических свойств связки, 
модифицированной алмазным порошком. 

Наилучшие показатели отмечены у образцов с наполните- 
лем из УДПА. 

Обсуждение и заключения. В результате проведенного 
исследования отмечено, что расчетные данные отличают- 
ся от экспериментальных, т.к. показывают увеличение 
упрочнения материала пропорционально количеству вво- 
димых в объем частиц алмаза. Выдвинуто предположение, 
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5гез$ ШИ ап уте!4 уаше. Сопу4ейп? фе дааШайуе согге- 
1аноп Бебмееп уле!4 уаше ап4 Багапез$, И 15 роз$1Ые ю рге- 
Ч1сё ап шсгеазе ш Вагапез$ шса те у 11214 таена|5. 
Миепта[5 ап Мефо4$. Рнуз1са| свагасетзНс$ оЁ фе Баз1с 
(урез оЁ ЯПегз ШФаё таКе пр Фе гешгсегз аге соп$14егеа, 
зузетанте4 апа 1абщае4. М2-01 ип Бгопте (20 \1. % ип, 
80% соррег) \аз зе аз а Бопа. ОШга1зретзе пага| Ча- 
топа (ООМЮ, 0.5— у/. %) уаз аадед о #, аз ме аз ром- 
Чегз оЁ пага! Ф1атопа (3/2 рла Насйоп, 7/5 ит, —40 рп) 
оБашей гои?Б ргосеззш» Фатоп4$ аё фе ешщегризе оЁ 
“ЗаКва Райпопа” ]$С. ТБе абоуе та{ета|$ \уеге та4е оп 1е 
сгазН ше ап 5сгеепше еддиртепе ап зВакш» {фа ез. ТБе 
арез ог оМашше ро\у4егз$ ууеге гесог4е4 изше Ше тгазег 
е!есёгоп псгозсору. У1гозсгееп$ \меге аррПе4 г Фе эташ- 
$12е с1аз Исай оп оф Фатопа роууаегз. Рнуз1са] ап плесБап1- 
са| свагасмет$Нс$ ог Фе ргодисе4 затр]ез \меге 1езеа Бу 
убапдаг4 тефоаз. УГТЕ-500 еесгот1с Юиг-с1азз 1абогаю- 
гу Ба|апсе ууаз изе4 ог \е1е те. ОепзИу \’аз деегтшеа Бу 
МК 0-25 шт пистотеег ассог4те №0 СОЗТ 6507-78. 
Вебеагсй Кези5. Рогозйу уаз сасшае4 @тгоиеН Фе асаа| 
апа Феогейса] депзНу. [Е \’аз юпа фай ИВ а десгеазе ш фе 
ВЦег $17е, ап паргоуетет ш фе рБуз1сотесВашса! ргоре!- 
(ез оЁ Ше Би4дег тофШеа у а Фатоп ро\м4ег 1$ о- 
зегуеа. ТЬе Ъезё реоппапсе \уаз обзегуе4 ш фе запарез 
у Ше ООМО ВПег. 

Бу5сизяюпт апа Сопсш51оп5. Аз а тез оЁ 1е заду, # уаз 
гесог4е4 фаё фе са|си]а{ёеа Чаёа @1ег от Бе ехрегилепа! 
аа зшсе Шеу зНо\ ап шсгеазе ш Фе таета| Ватдепте рго 
га{а ю фе атопие оЁР Фе Фатопа рагасез шеодасе4 шо Ше 
уоте. Ап аззитрНоп Ваз Бееп па4е а Фе сопз14еге4 
Баг4епте то! (Ого\уап по4е]) 4оез пой {аке шю ассоциЕ 


^" Е-тай: тайзаюп_2006@тай.га, аКапос!Кауа@уапаех.га, ЕедоЁ апагеу@тай.га 


^^" Тре гезеагсВ 15 допе у ве #гате оР\Ше ш4ерепдеги В &О. 
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что рассмотренная модель упрочнения (модель Орована) 
не учитывает образование углерода и агломерацию алма- 
зов в более крупные объекты в объе 

ме матрицы при повышении количества вводимых алма- 
зов. Если объем частиц УДПА достигает 3 %, в материале 
растет содержание углерода. В результате частицы напол- 
нителя полностью не окисляются, тем самым увеличивая 
количество пор в материале. 


Ключевые слова: связка, металлическая матрица, компо- 
зит, упрочнитель, ультрадисперсные частицы, механизмы 
упрочнения. 


Образец для цитирования: Сафонова, М. Н. К вопросу об 
особенностях механизмов упрочнения в металлической 
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Введение. Известно, что физико-механические свойства алмазных порошков определяются большей 
удельной поверхностью и дисперсностью зерен. В частности, это демонстрирует уравнение Холла — Петча, 
которое выполняется в широком диапазоне размера зерен (вплоть до 1 мкм). Функциональные свойства указан- 
ных материалов обусловлены масштабным структурным фактором, поскольку он влияет на формирование 
структурно-чувствительных механических свойств — пределов прочности и текучести [1, 2]. Учитывая каче- 
ственную корреляцию между пределом текучести и твердостью, можно прогнозировать повышение твердости, 
в том числе высокотвердых материалов: чем мельче наполнитель, тем меньше в нем дефектов — и, соответ- 
ственно, выше прочность. 

Материалы и методы. Механизмы упрочнения в высокотвердых материалах зависят от характера 
взаимодействия введенных частиц или волокон упрочнителя с материалом матрицы. Успешное применение 
эффекта дисперсного упрочнения показано в [3-6, 7]. При таком упрочнении в материалах создается структура, 
затрудняющая движение дислокации. Особенно сильно передвижение дислокаций тормозят дискретные части- 
цы второй фазы, характеризующиеся высокой прочностью и температурой плавления. С учетом двухфазной 
структуры и высокой твердости получаемых материалов следует ожидать, что их износостойкость также будет 
выше, чем у неупрочненных материалов матрицы. Кроме того, для износостойкости композитов важным фак- 
тором является прочность границы раздела компонентов. Прочная адгезионная связь на границе раздела обес- 
печивает получение композита с большой жесткостью и более высокой статической прочностью [3, 6-8]. 

Согласно механизму, описанному Орованом, распределенные в объеме связки ультрадисперсные ча- 
стицы алмаза влияют на деформационное упрочнение композиционного материала. Введенные в матрицу ча- 
стицы дисперсного наполнителя тормозят движение дислокаций в металле, увеличивая его прочность при нор- 
мальной и повышенных температурах. К тому же они представляют собой механическое препятствие на пути 
распространения трещины, которая может появиться в матрице, и повышают сопротивление разрушению ком- 
позиционного материала. 

Еще одной характеристикой, определяющей связь армирующего наполнителя и матрицы, является 
температурный коэффициент линейного расширения. Для твердых веществ при постоянных давлении и темпе- 
ратуре критерием термодинамического равновесия является минимум потенциала (или энергии) Гиббса. Эта 
величина показывает изменение энергии в ходе химической реакции и демонстрирует возможность протекания 
химических реакций между компонентами материала [5]. Таким образом, минимум изменения энергии Гиббса 
соответствует устойчивому равновесию между компонентами системы (табл. 1). 





























Таблица 1 
Физические характеристики основных типов наполнителей [9, 10, 11] 
Вещество АО’об,. 298,15, кДж/моль Микротвердость, х10* Температурная 
и состояние Изменение энергии Гиббса Мпа устойчивость, °С 
С (алмаз) 2,377 1000 650—700 
и 0 258 3300-3400 
АБОз —1582,3 180-220 1500-1700 
ВМ —226,8 800-900 1100-1300 
51С —60 300-320 1200-1300 
Вед —579,9 152 2500 
Ве›С —948 780 2150 
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Из таблицы видно, что оптимальный наполнитель — алмаз. У него достаточно низкое значение изме- 
нения энергии Гиббса, самое большое значение микротвердости, но наименьшая температурная устойчивость. 
Алмаз обладает большой адсорбционной способностью [12] и наименее химически активен в сравнении с дру- 
гими формами углерода. Названные качества являются важными преимуществами при использовании алмаза в 
качестве упрочнителя. 

Цель работы — исследование особенностей механизмов формирования структуры металломатричной 
композиции. 

Объекты исследования. В работе использовали стандартную связку — оловянистую бронзу М?2-01 
(20 мас. % олова, 80 % меди). В нее добавляли ультрадисперсный природный алмаз (УДПА, 0,5—4 мас. %), а 
также порошки (фракции 3/2 мкм, 7/5 мкм, —40 мкм). 

Методика проведения исследований. Порошки природного алмаза (ППА, фракции 3/2 мкм, 7/5 мкм, 
—40 мкм) и субмикропорошок УДПА получены при переработке алмазов на предприятии ОАО «Сахадаймонд». 
Порошки изготавливались при оптимальных режимах на дробильно-классификационном оборудовании и виб- 
ростолах. На рис. 1 показаны стадии получения порошков. Снимки выполнены с помощью растровой элек- 
тронной микроскопии (РЭМ). 


[72 $. 
'АссУ Эро! Маф \0 Г БО 
10.0КУ 5.0 500х 13.8 $ер4-1 
Ц 





Рис. 1. РЭМ-фотографии: мелкозернистого порошка от 7 мкм и менее (а); субмикропорошка (6) 


Для классификации алмазных порошков по зернистости применяли вибросита. Физические и механи- 
ческие характеристики изготовленных образцов испытывали по стандартным методикам. Для взвешивания ис- 
пользовали лабораторные электронные весы четвертого класса ВЛТЭ-500. Плотность р определяли микромет- 
ром МК 0-25 мм по ГОСТ 6507-78. 

Результаты исследования. Твердость измерялась на приборе Едиойр 3 фирмы Ргосед по соответству- 
ющей методике ГОСТ. Предел пропорциональности при сжатии и модуль упругости определялись по 
ГОСТ 25.503-97. Результаты испытаний приведены в табл. 2. 



































Таблица 2 
Влияние ППА на физико-механические свойства сплава М2-01 
Предел Фактическая МОДУ 
упругости 
ПА, % Пористость, % Твердость, НВ | пропорциональности плотность 
при сжатии оп, МПа р кг/м3 ыы 
ы ы ыы Ес, МПа 
а & < < < < < 
а юн аоозаюз в ао? Е Е 
Е 9 < [52 = т = [52] = ы г сп | ы г сб = ы ы г 
0 41,0 | 41,0 | 41,0 | 41,0 | 41| 41141141 | 9,7 | 9,7 | 9,7 | 9,7 | 7560 | 7560 | 7560 | 7560 5735,94 
1 29,0 | 30,0 | 30,0 | 28,0 | 43 | 47 | 48 | 48 | 10,5 | 11,5 | 11,0 | 11,5 | 7630 | 7620 | 7610 | 7670 6771,03 
2 27,0 | 28,0 | 28,0 | 26,0 | 43 | 43 | 51 | 53 | 10,5 | 11,5 | 12,0 | 12,5 | 7710 | 7700 | 7760 | 7750 6953,50 
3 25,0 | 28,0 | 27,0 | 26,0 | 47 | 43 | 51 | 53 | 10,2 | 10,5 | 12,1 | 12,2 | 7730 | 7720 | 7770 | 7750 6580,27 















































Чтобы определить модуль упругости и предел пропорциональности, образцы сжимали на прессах до 
относительной деформации, равной 15-16%, при скорости нагружения 0,2 кН/с. 
Пористость рассчитывали через фактическую и теоретическую плотности по формуле: 
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П= (1 — р/р:) х 100%, (1) 


где р, — теоретическая (расчетная) плотность беспористого материала; р — фактическая плотность исследуе- 
мого образца. 
Теоретическую плотность получали по формуле: 
р» = 100/(СИр1 + С>/р> + Су/рз + С4/ра), (2) 


где Су, С›, Сзи С. — концентрации меди, олова, ШТА и УДПА в порошковых смесях соответственно их плот- 
ности ру, Р›, Рзи ра. 

В расчет приняты плотности: 

— медь 8,96 х 103 кг/м°, 
— олово 7,28 х 103 кг/м®, 
— алмаз 3,5 х 103 кг/м, 
— УДПА 3,1 х 103 кг/м®. 

Объем образцов вычислялся по формуле: 

И= (т! — т?) / р (3) 
где Г — объем прессовки; т! — масса на воздухе; т› — масса в воде; р — плотность жидкости. 

В процессе исследования выяснилось, что с уменьшением размера фракции наполнителя наблюдается 
улучшение физико-механических свойств связки, модифицированной алмазным порошком. 

Значения модуля упругости определялись для образцов с добавлением УДПА и чистой связки. 
Наилучшие показатели демонстрируют образцы с наполнителем из УДПА. В то же время рассматриваемые 
физико-механические свойства ухудшаются, если объем добавленных частиц УДПА превышает 2 %. 

Металлографические исследования образцов позволили установить, каким образом частицы алмаза 
влияют на структуру матрицы. На рис. 2 представлены снимки микроструктур деформированных образцов. 





Рис. 2. Снимки шлифованной поверхности деформированного образца 
с добавлением 2 % частиц алмазных порошков при увеличении х1000 


На снимках четко различимы микроструктурные узкие и разветвленные объекты. Предположительно, 
это границы между зернами или микротрещины, образовавшиеся при деформировании. 

Внутри зерен также присутствуют точечные микрообъекты, образующие дисперсную субструктуру. По 
сравнению с гранично-распределенными точечными объектами их плотность значительно меньше, но ощутимо 
больше, чем у исходной матрицы, в которой отсутствуют алмазные наполнители. 

Основываясь на результатах металлографических исследований, можно утверждать, что при упрочне- 
нии материала матрицы действуют два механизма — дисперсный и зернограничный. 

Если речь идет о дисперсном механизме, объем вводимого в материал матрицы порошка природного 


алмаза можно рассчитать по уравнению Орована [13]: 


[67°] А 
ею Ш (4) 


где ^. — ближайшее расстояние между частицами; С — модуль сдвига матрицы; Ь — вектор Бюргерса; К, — 
коэффициент, характеризующий тип взаимодействующих атомов с дислокацией. 

Выбраны следующие значения: С = 0,367х10° МПа для бронзы; В = 2,564 А для меди; коэффициент Ао 
равен 0,85. 
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Ближайшее среднее расстояние между частицами в зависимости от содержания и дисперсности рассчи- 
тано по формуле, предложенной в [14]: 


= (ион _ 1 а, (5) 


1,91ЕН 
где [н — весовая доля наполнителя; 4 — диаметр или толщина частиц наполнителя. 
В табл. 3 приведены расчеты среднего расстояния между частицами наполнителя в зависимости от их 
объема и размера. 
Таблица 3 
Ближайшее среднее расстояние между частицами наполнителя в зависимости от их объема и размера и упроч- 
нение по уравнению Орована при введении частиц алмаза 





























Зернистость частиц 7/5 3/2 —40 УДПА 
Содержание частиц, % 1 2 3 1 2 3 1 2 3 1 2 3 
^, мкм 125,3 | 98,19 | 85,02 | 52,19 | 40,91 |35,42 | 417,56 | 327,30 |283,40 | 6,26 | 4,91 | 4,25 
С» МПа 0,13 | 0,16 | 0,18 | 0,28 | 0,35 | 0,40 | 0,04 | 0,05 | 0,06 | 1,91 | 2,38 | 2,70 
































Расчетные данные были подставлены в уравнение Орована, и таким образом определено упрочнение за 
счет введения в материал матрицы частиц алмазных порошков. 

Согласно расчетам, наибольшее упрочнение достигается при введении в матрицу УДПА, что в целом 
подтверждается экспериментальными данными. 

При изменении геометрии зерна за счет агломерации частиц наполнителей на границах раздела в мате- 
риале целесообразно рассчитывать изменение свойств материала согласно теории зернограничного упрочне- 
ния [15, 16]. 

Для определения количественного увеличения прочности материала при добавлении частиц алмазных 
порошков за счет зернограничного упрочнения были проведены расчеты по эмпирическому соотношению Хол- 
ла — Петча [17]: 

Дот = Каз ?, (6) 


где А — коэффициент Холла — Петча для данного материала; 4, — размер зерна. 

Для расчетов использовались образцы, показывающие наибольшее увеличение прочности по теории 
Орована. Расчеты выполнены согласно данным, полученным по итогам обработки снимков микроструктуры 
поверхности, по методике, предложенной в [18]. 

Коэффициент Холла — Петча применен к меди. Согласно [19] он является величиной переменной, за- 
висит от среднего размера зерна и варьируется в пределах 0,01-0,24 Мпа х м”. Расчеты показали, что 
наибольшее упрочнение дает введение в материал матрицы порошков ультрадисперсного природного алмаза. В 
целом, это подтверждается экспериментальными данными. 

Средний размер зерна рассчитан согласно металлографическим исследованиям поверхности образцов: 


Ч) 
аз = ——- (7) 

ГДе 5общ — общая площадь объектов; М№6бщ — количество объектов. 

При значении среднего размера зерен порядка 10' мкм коэффициент Холла — Петча равен примерно 
0,01 Мпа х м". 

Расчеты по соотношению Холла — Петча свидетельствуют об увеличении предела текучести материа- 
ла при добавлении частиц порошков природного алмаза. Предел текучести достигает максимального расчетно- 
го значения при содержании наполнителей 1 % (рис. 3). 


15 


Дот, МПа 
<л 


0 1 2 3 4 
Содержание УДПА ‚,% 


Рис. 3. Зависимость предела текучести от содержания УДПА 
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Обсуждение и заключения. Расчеты показали, что у образцов с добавлением порошков природных 
алмазов размер зерен меньше по сравнению с исходными. Этот факт можно объяснить следующим образом. 
Частицы алмаза, оседая на границах зерен материала, способствуют уменьшению их среднего размера. В ре- 
зультате меняется геометрия границ между зернами, образуются барьеры для дислокаций, и, как следствие, 
усиливается потенциальная способность материала сопротивляться пластической деформации. Предел текуче- 
сти увеличивается примерно на 12—13 МПа, что коррелирует с расчетными данными, полученными при исполь- 
зовании теории Орована для дисперсионного упрочнения. Если речь идет о поликристаллическом материале, 
то, в принципе, предел текучести увеличивается при уменьшении размеров зерен. Добавленные в объем матри- 
цы частицы алмаза увеличивая предел текучести, т. к. меняют геометрию зерен, уменьшая их среднюю пло- 
щадь и размер. 

Расчетные данные отличаются от экспериментальных, т. к. показывают увеличение упрочнения мате- 
риала пропорционально количеству вводимых в объем частиц алмаза. Можно предположить, что модель 
упрочнения Орована не учитывает образование углерода и агломерацию алмазов в более крупные объекты в 
объеме матрицы при повышении количества вводимых алмазов. 

Уменьшение количества пор при добавлении частиц УДПА в количестве 1-2 % можно объяснить вы- 
сокими сорбционными свойствами наполнителя. При спекании компактов, полученных методом порошковой 
металлургии, частицы УДПА абсорбируют кислород, содержащийся в порошковой смеси, с образованием га- 
зов-восстановителей СО и СО.. Эти газы разрушают окисную пленку, покрывающую частицы порошковой 
смеси, и препятствуют окислению при спекании, тем самым уменьшая общий объем газов в порошковой смеси. 
В то же время газы-восстановители ускоряют процесс спекания материала. Сочетание этих факторов в итоге 
уменьшает остаточную пористость в материале, что подтверждается и расчетными данными. Если же объем 
частиц УДПА достигает 3 %, происходит увеличение содержания углерода в материале. В результате частицы 
наполнителя полностью не окисляются, тем самым увеличивая количество пор в материале. 
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Введение. Статья разработке 
экспериментальных методов оценки характера распределе- 


посвящена расчетно- 
ния сил по ширине зубчатого венца и в окружном направ- 
лении зубчатого зацепления крупной волновой передачи. 
Исследование основано на результатах тензометрии кон- 
трольного зуба с использованием масштабного моделиро- 
вания опытных образцов. Цель работы — создание объек- 
тивной экспериментально-теоретической модели силового 
взаимодействия зубьев в зацеплении крупной волновой 
передачи. Такое решение предполагает преобразование 
деформаций модельных образцов в распределенную 
нагрузку между зубьями, что позволит исключить основ- 
ные неконтролируемые погрешности нелинейного характе- 
ра, повысить точность оценки силовых факторов в поле 
зацепления. 

Материалы и методы. Разработана уточненная методика 
силового анализа зубчатого зацепления крупной волновой 
передачи, оптимизированная по критериям точности. Точ- 
ность результатов исследования повышена за счет совер- 
шенствования физических и расчетных моделей. Такой 
подход позволил получить обоснованные зависимости рас- 
пределения силовых факторов в зубчатом зацеплении 
крупной волновой передачи. 

Результаты исследования. Упрощена конструктивная фор- 
ма контрольного зуба: по всей ширине зубчатого венца 
введен неизменный профиль. Таким образом исключены 
нелинейные искажения результатов экспериментов, вноси- 
мые переменной формы зуба по ширине зубчатого венца. В 
таком случае возможен монтаж тензорезисторов по всей 
ширине контрольного зуба. Кроме того, предлагаемое ре- 
шение позволяет установить зависимость деформации 
зубьев по всей ширине зубчатого венца, а не только на 
крайних участках, как предлагают известные методики. 
Разработка совершенных физических и математических 
моделей позволила повысить точность результатов теоре- 
тических и экспериментальных исследований силовых про- 
цессов в зубчатом зацеплении крупной волновой передачи. 
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Получены научно обоснованные двухпараметрические за- 
висимости распределения сил в зубчатом зацеплении. 
Обсуждение и заключения. Аппроксимация эвольвентного 
профиля зуба в трапецеидальный профиль упростила дока- 
зательство тождественности уравнений упругости и гра- 
ничных условий математических моделей. Полученные 
результаты применимы при математическом моделирова- 
нии плосконапряженного состояния зубьев с нелинейными 
профилями. Сравнительная оценка погрешностей, вноси- 
мых отклонениями геометрических форм и размеров физи- 
ческих моделей и математических аналогов, подтверждает 
корректность постановки эксперимента и обоснованность 
полученных количественных данных. Результаты работы 
могут быть использованы при уточненном расчете кон- 
структивных параметров элементов зубчатого зацепления в 
процессе проектирования крупных тяжело нагруженных 
волновых редукторов. 


Ключевые слова: волновая передача, зубчатое зацепление, 
силовой анализ, распределение нагрузки. 
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Введение. Промежуточное преобразование вращательного движения в непрерывную волновую дефор- 
мацию гибкого колеса изменило сложившиеся принципы и формы сопряжения элементов зацепления высших 
кинематических пар. При этом малая разность зубьев во внутреннем зацеплении волновой передачи (два зуба) 
не согласуется с традиционными представлениями об условиях геометрического синтеза и функционирования 
зубчатой пары. Поэтому с точки зрения интерференции зубъев зацепление волновой передачи находится за до- 
пустимыми пределами существования эвольвентной внутренней зубчатой передачи с жесткими зубчатыми ко- 
лесами. Возникающая в зацеплении интерференция зубьев второго рода усиливается за счет деформации гиб- 
кого колеса от передаваемой нагрузки, что ограничивает несущую способность волновой передачи. Негативное 
влияние масштабного фактора способствует возникновению интерференции зубьев в крупных волновых пере- 
дачах, что вызывает заклинивание и проскок зубьев в зацеплении. 

Базовая концепция волновой зубчатой передачи (5маш \ауе сеаппе, З\С) запатентована 
К. У. Массером в 1957 году. Решение остается неизменным для серийных образцов ограниченной мощности, 
несмотря на множество конструктивных усовершенствований. Подтверждением этому служат волновые редук- 
торы, выпускаемые ведущими фирмами США (Нагтошс Опуе Тесвпо]оз1ез шс., Опиед Звое МасЫпегу Сотр.), 
Японии (Нагтошс Опуе Зузетз шс., Германии Нагтопс Опуе АС) и других, в том числе лицензионная про- 
дукция Наппошс Опуе Зу\етз шс., производимая в КНР [1, 2]. 

В крупных волновых зубчатых редукторах реализуются большие мощности, вращающие моменты пре- 
вышают (0,3...1,5)х10° Нм. Такими деталями комплектуется горнорудное, металлургическое оборудование и 
большие машины, выпускаемые Новокраматорским машиностроительным заводом (НКМЗ). Типовая конструк- 
ция кулачкового генератора волн с гибким шарикоподшипником неработоспособна и неприменима в подобной 
продукции тяжелого машиностроения [3]. В крупных волновых редукторах используется дисковый генератор 
волн. В отличие от кулачкового генератора, одновременно контактирующего по всему периметру гибкого коле- 
са, дисковый взаимодействует с гибким колесом на диаметрально противоположных участках. Это создает бо- 
лее высокий уровень свободы деформации гибкого колеса и увеличивает возможные отклонения от заданного 
положения. 

Под действием дисков генератора волн и вращающего момента гибкое колесо приобретает искривле- 
ние в виде слегка закрученной конусности. Это обусловливает перекос зубьев по ширине зубчатого венца отно- 
сительно жесткого колеса. Нарушается равномерность величины зазоров в зубчатом зацеплении по ширине 
зубчатых венцов. При больших нагрузках деформации гибкого колеса превышают зазоры в зубчатом зацепле- 
нии и создают условия интерференции зубьев второго рода. При широких зубчатых венцах, малом модуле 
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зубьев в условиях больших нагрузок зазоры в локальных зонах деформированного зубчатого зацепления при- 
нимают отрицательные значения, активизируют интерференцию, вызывают заклинивание и проскок зубьев. 
Конструктивные отличия крупных волновых зубчатых передач обусловлены главным образом масштабом. Этот 
фактор усугубляет негативные процессы, наблюдаемые в зубчатом зацеплении высших кинематических пар, и 
снижает технические и эксплуатационные характеристики тяжело нагруженных волновых редукторов. 

Целью настоящей работы является повышение нагрузочной способности, технических и эксплуатаци- 
онных характеристик крупных волновых зубчатых редукторов, удовлетворяющих техническим условиям тяже- 
лого машиностроения. 

Теоретические и практические наработки по исследованию крупных волновых передач ограничены, 
полученные результаты бывают противоречивы, не всегда повторяются. Кроме того, необходимо отметить су- 
щественные конструктивные отличия, отсутствие геометрического подобия, в том числе в размерных рядах 
кинематических и крупных волновых передач. Все это не позволяет использовать известные расчетные методи- 
ки в проектировании волновых передач для целей тяжелого машиностроения [4—6]. 

Синтез зубчатого зацепления подчинен условию постоянства заданного передаточного отношения, ре- 
гламентируемому теоремой Эйлера — Савари. Перманентная деформация гибкого колеса не согласуется с ос- 
новной теоремой зубчатого зацепления, что способствовало возникновению противоречивых методик синтеза 
волнового зубчатого зацепления [7-11]. В одних работах за основу принимается общий случай зубчатого за- 
цепления [12], в других рассматривается клиновой механизм со сложным относительным движением звень- 
ев [13-15]. В некоторых методиках упрощается геометрия и кинематика волнового зубчатого зацепления [16, 
17]. Комбинированный подход представляет собой синтез общей теории зубчатого зацепления с элементами 
клинового механизма [1]. 

Для удовлетворения условиям основной теоремы зубчатого зацепления применительно к синтезу вол- 
новой передачи Е. Г. Гинзбург использовал кинематический подход [18]. Введенное им понятие «угловая ско- 
рость точки» ошибочно, основанные на нем доводы некорректны. 

В монографии Н. А. Ковалева [19] отмечены сложность силовых и кинематических процессов, проте- 
кающих в волновых передачах, отсутствие достоверных зависимостей, адекватно отражающих влияние внеш- 
них факторов на внутренние процессы в кинематических парах. Это не позволяет разработать эффективный 
метод синтеза волновой зубчатой передачи, оптимизирующий основные параметры несущих звеньев. По 
Н. А. Ковалеву, при нерастяжимой срединной поверхности равномерно движется лишь точка пересечения оси 
зуба и срединной поверхности, что противоречит выводам Е. Г. Гинзбурга [18]. 

Величину максимального радиального упругого смещения и’, срединной поверхности гибкого колеса 
задает генератор волн при отсутствии крутящего момента. Смещение и’, должно превышать половину рабочей 
высоты зуба. Размеры зацепления следует корректировать по величине у», не наоборот [19]. 

Относительно значения параметра у. не существует единого мнения. Ряд авторов связывает макси- 
мальную радиальную деформацию ии». с передаточным отношением волновой передачи [20] и дисковым гене- 
ратором волн. С целью повышения КПД рекомендуется принимать и’. = т. Е. Г. Гинзбург рекомендует опре- 
делять И’мах Через условный модуль т, отличающийся от модуля т [18]. 

Материалы и методы. Для предупреждения интерференции зубьев в серийно выпускаемых волновых 
редукторах применяют известные приемы: 

— корригирование зубьев, 
— увеличение ширины впадины между зубьями, 
— увеличение угла исходного контура зубьев. 

Это наряду с малой шириной зубчатых венцов и сравнительно небольшими передаваемыми крутящими 
моментами во многом решает практические вопросы приближенного синтеза зубчатого зацепления в части 
устранения интерференции зубьев для серийно выпускаемых малонагруженных волновых редукторов. Однако 
теоретические исследования этих вопросов далеки от завершения. Пока нет удовлетворительного решения от- 
носительно зацепления высших кинематических пар крупных волновых редукторов с шириной зубчатых вен- 
цов 100-200 мм, модулем зубьев 1,5-=3 мм и дисковым генератором волн. Спроектированные по известным ме- 
тодикам крупные волновые редукторы, изготовленные на НКМЗ, оказались неработоспособными из-за интер- 
ференции зубьев. Решение подобных задач требует новых технических решений, основанных на достоверных 
результатах теоретических и экспериментальных исследований. 

Для оценки граничных условий проявления интерференции и проскока зубьев в зацеплении, а также в 
целях предотвращения этих негативных явлений проведены научные изыскания. Их результаты позволили со- 
здать обобщенную модель взаимодействия зубьев в зацеплении крупной волновой передачи. Разработана уточ- 
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ненная методика силового анализа такого зацепления. Точность результатов исследования повышается за счет 
совершенствования физических и расчетных моделей. Это позволяет получить обоснованные зависимости рас- 
пределения силовых факторов в зубчатом зацеплении крупной волновой передачи. 

Деформацию напряженных элементов конструкций, в том числе зубъев, принято измерять методами 
тензометрии. Однако на зубьях волновой передачи с небольшим модулем невозможно закрепить тензорезисто- 
ры. По известным методикам [20, 21], в ободе жесткого колеса выполняют две параллельные прорези, удаляя 
при этом два смежных зуба (рис. 1). 





Рис. 1. Локализованные элементы зуба жесткого колеса редуктора наклона передвижного миксера МП-600АС 


По ширине колеса прорези могут быть сквозными или глухими, выполненными с обоих торцов обода 
таким образом, чтобы средняя часть осталась нетронутой. Согласно утверждениям авторов распространенных 
методик, сквозные глубокие прорези могут существенно уменьшить жесткость колеса, что недопустимо. В из- 
вестных работах применяется экспериментальная модель, где средняя часть зубьев не удалена и работает наря- 
ду со всеми другими. 

Утверждения о недопустимом снижении жесткости колеса сквозными прорезями [22] неубедительны, 
т. к. толщина обода жесткого колеса не регламентируется и в условиях эксперимента ее можно принять доста- 
точной, чтобы прорези не оказывали заметного влияния на жесткость. Кроме того, жесткое колесо запрессовано 
в конструкцию корпуса и при оценке жесткости рассматривается совместно с корпусом редуктора. 

По методикам [20, 21] прорези жесткого колеса увеличивают деформацию контрольного зуба относи- 
тельно его срединной части. Изменение жесткости контрольного зуба по длине повышает деформацию выде- 
ленных элементов у торцов, искажает деформированное состояние контрольного зуба и не позволяет с доста- 
точной точностью оценивать силовые характеристики зубчатого зацепления по ширине венца и в окружном 
направлении. Многократное снижение жесткости выделенных элементов зуба не учитывается известной мето- 
дикой [22], что нарушает объективность результатов эксперимента и осложняет их анализ. 

Типовая методика измерения деформации зубьев волновой передачи [8] мало отличается от методики 
обычных зубчатых передач [23]. В [22] показано, что перепад жесткости модели, обусловленный глухими про- 
резями на жестком колесе, искажает деформацию контрольного зуба при нагрузке и не учитывается при прове- 
дении эксперимента и обработке полученных данных. Таким образом снижается точность и достоверность по- 
лученных результатов. 
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Тензорезистор не крепится посередине контрольного зуба, поэтому невозможно получить полное пред- 
ставление о его деформации, качественные и количественные характеристики распределения сил в зубчатом 
зацеплении. Кроме того, при рассмотрении ступенчатого изменения жесткости физической модели по извест- 
ной методике не учитываются граничные условия, влияние которых усиливается фактором масштаба [18]. 

Рассмотрим волновой редуктор ВЗ-1120 привода наклона передвижного миксера МП-600АС емкостью 
600 т расплавленного металла. Примем допущение, что нагрузка в зацеплении волновой передачи приложена 
на середине высоты зуба, т. к. во внутреннем зацеплении одной волны находится - 20-25 % зубьев. 

Определим возможную ошибку допущения в объеме эксперимента. Рассмотрим жесткое зубчатое ко- 
лесо с внутренними зубьями: 

— модуль зубьев т = 1,5 мм; 

— число зубьев жесткого колеса 2 = 762; 

— число зубъев долбяка 1 = 68; 

— диаметр делительной окружности колеса 4 = 1143 мм; 

— высота зуба колеса Н = 2,985 мм; 

— диаметр окружности вершин зубьев колеса О = 1155,12 мм; 

— диаметр окружности впадин зубьев колеса О,= 1161,09 мм; 

— угол профиля в точке зуба на делительной окружности 4, о = 20°; 
— коэффициент смещения исходного контура инструмента х = +4,953 мм; 
— ширина зубчатого венца Б = 100 мм; 

— нормальный исходный контур по ГОСТ 13755-81; 

— степень точности 70 (ГОСТ 1643-81). 

Относительная погрешность допущения К, касающегося условного приложения нормального усилия на 
середине высоты зубьев: 

_ 2Н 100% — 2.2,985 
р-р, 1155,12-+1161,09 


Принятое допущение не оказывает заметного влияния на передачу усилий В зубчатом зацеплении, т. к. 


К 








-100% = 0,3%. 


незначительно изменяет диаметр приложения нагрузки. Для аппроксимации эвольвентного профиля зуба в 
прямолинейную форму определим толщину зуба 5, жесткого колеса по произвольной окружности радиуса 7’: 


с05 п р ? в 
5х =т—бб+А, + 25 (туа,, — туб }, с05 4, = —с05 90, 
с05 9 7; 


где о. — угол профиля в точке профиля зуба, в которой радиус — вектор г; А› — коэффициент изменения 
толщины зуба по любой окружности колеса в нормальном сечении при нарезании зубъев долбяком средней 
изношенности: 


Д> = (2 — 74) (ту ад — ту ас). (1) 
Относительная погрешность определяется как модуль отношения абсолютной погрешности к значению 
приближенной величины. Максимально возможное значение абсолютной погрешности в рассматриваемом слу- 


чае равно г Приближенной величиной является радиус окружности середины зуба жесткого коле- 


са(р+2, ) / 4. Косинус угла станочного зацепления ас в выражении (1) при нарезании жесткого колеса дол- 
бяком средней изношенности определяется с учетом коэффициента смещения исходного контура долбяка: 


22 — 24 
22—24 + 252 | 





С050С2 = СО50\) 


Значения ТОЛЩИНЫ зуба 5х по дуге произвольного радиуса г, жесткого колеса приведены в табл. 1. 


Таблица 1 
Толщина зуба 5х жесткого колеса в различных сечениях, вычисленная по дуге окружности радиуса гх, мм 





т: 579,239 579,426 579,612 579,799 579,985 580,172 580,358 580,545 
































5, 2,613 2,765 2,918 3,071 3,223 3,377 3,530 3,685 





Толщина зуба по дуге окружности отличается от толщины зуба по хорде при заданных параметрах 
жесткого колеса на величину 5:10 мм, поэтому с достаточной точностью можно приравнять толщины зуба по 
дуге окружности радиуса 7’, и по соответствующей ей хорде. 
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Определим толщину зуба 5, аппроксимированного в прямолинейный профиль на радиусе 
ух = 579,799 мм. Она равна полусумме крайних значений толщины зуба по радиусам г, = 579,053 мм и 
г, = 580,545 мм: Зах = 3,074 мм. 

Толщина зуба эвольвентного профиля жесткого колеса на окружности радиуса г’. = 579,799 мм равна 
5, = 3,071 мм. Абсолютная погрешность аппроксимации профиля зуба жесткого колеса в прямолинейный про- 
филь: 

А =5, -5х =-0,003 мм. 

Относительная погрешность прямолинейной аппроксимации эвольвентного профиля зуба жесткого ко- 

леса: 


— ох — Зах 





В -100 = 0,1%. 
ах 
Аппроксимация эвольвентного профиля зуба жесткого колеса прямолинейным профилем не вносит за- 
метную погрешность в расчетную схему. Уравнение прямолинейного профиля модели зуба в координатах ХОУ 
представлено на рис. и 
ТИ. _И- 
Х-х, Х-хХ, 





(2) 


1, = 15,947 

















Рис. 2. Расчетная модель контрольного зуба жесткого колеса с аппроксимированным прямолинейным профилем 


Индексами | и 2 на рис. 2 указана принадлежность координатных функций Х и У крайним точкам пря- 
молинейного профиля зуба 1 и 2. В уравнение (2) подставляются значения Х, У в точках | и 2 согласно рис. 2, и 
записывается уравнение прямолинейного профиля зуба с линейными размерами в миллиметрах: 

У =7,767-—0,410 Х. 

Толщина и жесткость аппроксимированного прямолинейного зуба жесткого колеса в произвольном се- 

чении, мм: 
5х =2У =15,534—0,820Х. (3) 


3 
ЕТ, =. (4) 


Удельная жесткость аппроксимированного прямолинейного зуба жесткого колеса в произвольном се- 
чении определяется подстановкой в выражение (4) значений (3) при [= 17,5 мм: 


ЕТ, =0,169-105-(18,943—Х)". 
Прогиб упругой линии модели контрольного зуба под силой №, определяется при Х = [. = 15,947 мм: 
У, = 0,819-10* №, [мм]. (5) 
Напряжения растяжения с’ и сжатия о` модели контрольного зуба в области установки тензорезисто- 


ров на расстоянии а, = 7,2275 мм от консольного закрепления (защемления) модели: 


Хх 
а АМИ 20 олвбтм,- 9, 
И’, Е мм (6) 
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5. М Н 
с =—И+- 09 = (,1983М,,—_, (1) 
и, Е мм? 


Здесь М, — изгибающий момент в области установки тензорезисторов; Ми = М. - У - ат; Й’, — осевой момент 
сопротивления сечения линейной модели зуба: 


2 
Й’. = = = 39,61 мм. 


С помощью соотношения 
о = РЕ, 
где = — относительная деформация тензорезистора, а также формул (5-7) устанавливается зависимость между 
относительной деформацией тензорезисторов на растянутых и сжатых поверхностях модели контрольного зуба 
и прогибом упругой линии под силой: 
У, = 5013,33-=*, У, =472,07-=. (8) 

Результаты исследования. Теоретические и экспериментальные исследования, практические наработ- 
ки связаны преимущественно с освоением, производством и эксплуатацией сравнительно небольших волновых 
редукторов [25-32]. Накопленный опыт невозможно использовать без учета масштабного фактора для крупных 
волновых редукторов с крутящими моментами 5-10” Нм и более. Гибкое звено усиливает отрицательное влия- 
ние масштабного фактора, искажая обычный характер функционирования зубчатого зацепления высших кине- 
матических пар. При больших нагрузках интерференция второго рода может вызвать заклинивание и проскок 
зубьев. Интерференция также активизирует силовые процессы в области генератора волн, повышая энергетиче- 
ские потери. 


Фрагмент осциллограммы соответствует обороту генератора волн 2л (рис. 3). 


2 





Рис. 3. Осциллограмма деформации контрольного зуба жесткого колеса при нагрузке М. = 3: 10° Нм и частоте вращения на 
входе ; = 1,667 с; 1,2, 3, 4, 5 — осциллограммы деформации локализованных элементов контрольного зуба 


Число пар зубьев в зацеплении 7; определяется продолжительностью зацепления контрольного зуба в 
течение полного оборота генератора волн (© + а"): 


# 


" 
2 =“ 75 .100%. (9) 
2 

Число одновременно зацепляющихся пар зубьев и характер распределения нагрузки по ширине зубча- 
того венца зависят от конструктивных параметров ВОЛНОВОЙ передачи, определяющих условия зубчатого зацеп- 
ления, и представляют некоторые функции вращающего момента. 

Разработана методика силового анализа зубчатого зацепления крупной ВОЛНОВОЙ передачи. От извест- 
НЫХ решений ее отличают два преимущества. 

1) Она позволяет устранить искажения деформированного состояния контрольного зуба при нагрузке 
за счет постоянства жесткости по длине. С этой целью выполняются сквозные параллельные прорези на жест- 
ком колесе симметрично оси контрольного зуба. 
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2) Открывается возможность существенно повысить точность оценки силовых характеристик зубчатого 
зацепления по ширине венца и в окружном направлении благодаря учету снижения жесткости контрольного 
зуба. Для этого разрабатываются физические модели контрольного зуба и масштабированного зуба жесткого 
колеса. 

Обсуждение и заключения. Форма контрольного зуба принята неизменной по всей ширине зубчатого 
венца (это отличает представленный подход от методик [20-22]). Таким образом исключается искажение ре- 
зультатов эксперимента, вносимых переменной по ширине формой зуба. Решение обеспечивает возможность 
монтажа тензорезисторов по всей ширине контрольного зуба, что устанавливает полную зависимость деформа- 
ции зубьев по всей ширине зубчатого венца, а не только на крайних участках, как предусматривает известная 
методика. 

Совершенствование физических и математических моделей позволило повысить точность результатов 
теоретических и экспериментальных исследований силовых процессов в зубчатом зацеплении крупной волно- 
вой передачи [8-10]. Получены научно обоснованные двухпараметрические зависимости распределения сил в 
зубчатом зацеплении. 
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Математическое моделирование температуры в зоне контакта инструмента и изделия 
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Введение. Сопоставлены два фактора токарной обработки 
металлов резанием: рассеиваемая температура и мощность 
необратимых преобразований в материале изделия и ин- 
струмента. Статья посвящена вопросам математического 
моделирования их связи. 

Материалы и методы. Математический аппарат основан 
на модификации уравнения Вольтерры, которая предпола- 
гает использование двукратного интеграла. Он показывает, 
как выделенная ранее при резании тепловая энергия влияет 
на текущее состояние температуры в зоне контакта ин- 
струмента с обрабатываемой деталью. В работе помимо 
предложенной новой базовой математической модели ис- 
пользованы результаты обработки экспериментальных 
данных о мощности необратимых преобразований и об 
измеряемой температуре в зоне контакта инструмента с 
деталью при токарной обработке металла. Эксперименты 
проводились на станке 1К625 и стенде ЗТО.201-1. В пакете 
МаЧаб была создана специализированная программа для 
обработки массивов информации, описывающих процессы, 
протекающие при резании (силы реакции, вибрации ин- 
струмента и мощность необратимых преобразований). В 
этой же программе выполнен расчет температуры в зоне 
контакта инструмента и детали. 

Результаты исследования. Проведена процедура парамет- 
рической идентификации предлагаемой базовой математи- 
ческой модели. Полученная модель показала высокую сте- 
пень близости экспериментальных данных о состоянии 
температуры в зоне резания и моделируемого уровня теп- 
ловой энергии, однако на начальном участке измеряемой 
температурной зависимости результаты двух этих подхо- 
дов не вполне согласуются. Это можно объяснить погреш- 
ностью экспериментального измерения температуры на 
основе оценки вырабатываемой термоэлектродвижущей 
силы (термоЭДС), которая формируется в результате обра- 
зования естественной термопары в зоне контакта инстру- 
мента и детали. 
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Обсуждение и заключения. Предлагаемая математическая 
модель позволяет адекватно описывать процесс преобразо- 
вания механической составляющей энергии резания в теп- 
ловую составляющую через показатель суммарной выделя- 
емой мощности механического взаимодействия в зоне ре- 
зания за все время обработки. 
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мых преобразований, обработка металлов, теплопровод- 
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Введение. Тепловые процессы, протекающие при обработке металлов на металлорежущих станках (как 
в металлах, так и в других средах), являются сложноописываемыми. Их анализ основан на уравнении тепло- 
проводности Фурье [1]. В случаях, связанных с обработкой металлов резанием, нужно иметь представление о 
температуре в зоне контакта инструмента и обрабатываемой детали [2—5]. При этом нет необходимости кон- 
тролировать распространение тепла за пределами данной зоны, так как на условиях обработки это не сказыва- 
ется. Аналогичные рассуждения можно применить к другому широко распространенному методу описания 
тепловых процессов, который учитывает формирование истоков и стоков тепла в детали [6]. Такие подходы не 
позволяют использовать динамику процесса резания в описании возникающих тепловых явлений в зонах кон- 
такта и создать систему уравнений, представляющих преобразование механической энергии в тепловую. 

При операциях с металлами на металлорежущих станках необратимые преобразования в детали и ин- 
струменте вызывают производство тепла, которое, согласно второму закону термодинамики, рассеивается в 
пространстве. Тепловые процессы в металлах более инерционны, чем динамика процесса обработки, поэтому 
рассеивание продолжается достаточно длительное время. Само производство тепла в рассматриваемом случае 
связано с мощностью необратимых преобразований [11], которая зависит от величины сил реакции на формо- 
образующие движения инструмента [2—5] в зоне контакта с деталью. При этом производство тепла является 
основным механизмом рассеивания выделенной энергии преобразования материалов [12]. 

Динамика процесса резания и вопросы ее моделирования напрямую связаны с представлением сил и 
реакций в координатах состояния процесса. В этом случае, как правило, учитывается только механическая со- 
ставляющая, что видно из зарубежных [13—18] и отечественных работ [19—21]. Таким образом, целесообразно 
ввести в моделирование динамики процессов обработки металлов координаты, описывающие изменение тем- 
пературы в зоне контакта инструмента и детали через учет необратимых преобразований. Такой подход позво- 
лит качественно повысить адекватность используемого математического аппарата. 

Материалы и методы 

Обоснование математической модели. Силы, препятствующие внедрению инструмента в материал 
обрабатываемой детали, имеют сложное пространственное расположение в системе координат, связанной с 
осями деформации инструмента. В литературе, посвященной анализу сил реакций, возникающих в зоне реза- 
ния, принято разложение такой реакции на составляющие [4—6, 10, 19—21]. В рассматриваемом нами случае 
будут справедливы следующие соотношения между составляющими сил реакции [5]: 


2 2 2 
В=Р2+Р2+Р2, (1) 


где в — общая силовая реакция со стороны процесса резания на формообразующие движения инструмента; Р. 
—Щ проекция силовой реакции на ось главного вращательного движения резания в направлении подачи; Р, === 


проекция силовой реакции, направленная по радиусу главного вращательного движения резания в вершине ре- 
зания, В, —Й проекция силовой реакции на ось, совпадающую со скоростью главного движения в вершине ин- 


струмента. 
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Соотношения между введенными силами Р,,Е,,Ё, зависят от многих факторов: геометрия резца, сте- 


пень износа резца и т. д. [4]. Так, в [5] при обработке острым резцом с параметрами у = 15°, ф= 45° и/,= 0° со- 
отношение в среднем равно 
РРР. = (0,3 — 0,4), (0,4 — 0,5). (2) 
Аналогичные рассуждения можно провести по анализу скоростей относительного движения инстру- 
мента вдоль обрабатываемой детали и получить общую скорость такого движения: 


И. +И2+И.2. (3) 


С учетом этих рассуждений определим мощность необратимых преобразований: 


[ы 8 2 2 [/2 2112 
№М=уЬ/ +В +Е/ 7; +, +Т,. (4) 


Допустим, что производство тепла при резании в металлорежущих станках линейно зависит от значе- 
НИЯ МОЩНОСТИ необратимых преобразований: 
Т, = М. (5) 


Здесь ТГ, — величина, описывающая приращения температуры в зоне резания (°С); Кг — коэффициент связи 


между мощностью необратимых преобразований в зоне резания и величиной приращения тепловой энергии в 
зоне контакта инструмента и детали размерностью с°С/Нм. Однако общее значение температуры в заданный 
момент времени и в заданной зоне контакта инструмента и детали будет определяться не только текущим зна- 
чением прироста тепла. В связи с рассеиванием тепла в материале детали и инструмента следует учитывать 
влияние: 

— мощности необратимых преобразований по всему предыдущему пути, 

— времени обработки изделия. 

Представим себе обработку в виде некоторого дискретного графика процесса, где (и) — текущее зна- 
чение мощности необратимых преобразований; Г[(п) текущее значение пути, пройденного инструментом при 
обработке; #(п) — текущее значение времени обработки. Согласно (5), в каждой рассматриваемой точке дис- 
кретного графика процесса обработки будет выделено приращение тепла, пропорциональное мощности необра- 
тимых преобразований в этой точке пространства и времени. Текущее значение температуры в точке и будет 
определяться как приростом тепла в этой точке, так и влиянием приростов температуры, происходивших ранее 
на пути, пройденном инструментом (Г) во время обработки (1). С увеличением времени и пройденного пути 
ослабевает влияние этих факторов на процесс выработки тепла в рассматриваемой зоне контакта, поэтому для 
математического описания удобен подход, основанный на использовании оператора Вольтерры второго рода 
[10]. Его невозможно применить напрямую из-за сложности описания распространения тепла в металлах, по- 
этому в качестве базовой модели примем мультипликативный критерий оценки влияния предыдущего прироста 
тепла на текущее значение в виде двойного интеграла: 

7, = ©, + Аг [[,@-Ви(п-ОмМтачам. (6) 

р 
Здесь Т, — значение тепла в зоне контакта инструмента и обрабатываемой детали; ®; — температура окру- 
жающей среды; иг (у-Г) — ядро оператора, характеризующее влияние выделенной ранее мощности необра- 
тимых преобразований по пути обработки на текущее значение температуры; и",(п-#) — ядро оператора, ха- 
рактеризующее влияние выделенной ранее мощности необратимых преобразований по времени обработки на 


текущее значение температуры; р — область интегрирования, характеризующая пространственно-временное 
рассеивание тепла; ‘у имеет размерность расстояния, м; | — имеет размерность времени процесса обработки, 


с. Сами ядра интегрального оператора представляют собой безразмерные коэффициенты усиления. 
Интегральный оператор, представленный в виде двойного интеграла (6), можно свести к кратному ин- 
тегралу следующего вида: 
(и) 1 


7. -0; + | ма-В ау [м (т-0мая. (7) 
0 0 


Пройденный инструментом путь является функцией времени, т. к. каждому значению времени можно 
поставить в соответствие значение пройденного пути. Следовательно, подход (7) справедлив. 
В качестве ядер предложенного в (7) интегрального оператора удобно использовать экспоненциальные 
функции следующего вида [10]: 
р-р) =е* 2 
и, (п-1) = г’ 


где 01, 06 — масштабирующие параметры интегрального оператора, подлежащие идентификации. 


(8) 
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Таким образом, предложенный в (7) интегральный оператор примет следующий вид: 
770) 1 
Т, 9,4 | с“ Эду ЭмМтат. (9) 
0 0 
Этот оператор имеет решение для стационарного случая, когда мощность необратимых преобразований 
является постоянной величиной № = М(®. 


Ат № 
0195 


т. =©, + (-е “4-е. (10) 


Как видно из (10), при стационарном процессе обработки температура в зоне резания экспоненциально 
нарастает с начального значения ©, до некоторого конечного значения, равного: 


Ат. № 
0192 | 


Т. =©, + (11) 


Проиллюстрируем эти рассуждения, задав ®, =25 °С; Аг= 0,0026 с°С/Нм; ви =0,03; @> =0,01; 
№ =20 Нм/с (рис. 1). 


Т2(С°) 


Им) 





Рис. 1. Реализация процесса стационарного нарастания температуры в зоне резания 


Как видно из рис. 1, стационарный тепловой процесс в зоне резания нарастает экспоненциально и, до- 
стигнув некоторого равновесного состояния, больше не изменяется. Реальные процессы, протекающие в ме- 
таллорежущих станках, не являются стационарными [11-12]. Исходя из этих соображений, процесс, представ- 
ленный на рис. |, можно рассматривать как движение, к которому в устойчивом случае стремится температура 
зоны контакта инструмента и детали. При этом уравнение (10) является порождающим уравнением для инте- 
грального оператора (9). 

Результаты исследования 

Параметрическая идентификация модели по экспериментальным данным и итоги моделирова- 
ния. Для параметрической идентификации модели проведен натурный эксперимент на токарном станке 1К625 
с установленным стендом 5ТО.201-1. Стенд предназначен для изучения динамических и тепловых процессов, 
протекающих при резании металлов на токарных станках в различных режимах. 

Основной целью эксперимента было исследование влияния динамики резания на тепловые процессы в 
зоне контакта инструмента и детали. Под изменением динамики резания понимается взаимное влияние износа 
инструмента и процесса резания, связанное с возникающей в зоне резания силовой реакцией. Измерение этой 
реакции и температуры через термоэлектродвижущую силу (термоЭДС), формируемую в зоне контакта ин- 
струмента и детали, позволяет оценить влияние динамики резания на тепловые процессы в указанной зоне. 

Во время эксперимента выполнялось непрерывное точение детали типа вал. При этом сохранялись все 
параметры обработки и осуществлялся контроль изменения силовой реакции со стороны процесса резания. 
Условия и основные результаты эксперимента подробно описаны в [22]. 
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Данные, полученные в результате натурного эксперимента, носят дискретный характер, следовательно, 
к ним нельзя напрямую применить модель, представленную интегральным оператором (9). По этой причине 
приходится видоизменить базовую математическую модель, описывающую динамку тепловых процессов в 
контакте инструмента и детали. Для начала оставим в рассуждениях только член интегрального оператора, за- 
висящий от мощности необратимых преобразований: 


Ци) 1 
ТА =Ю | ау е 0 М(мат. (12) 
0 0 


Так как эксперимент уже проведен, то время его окончания | и соответствующее этому времени зна- 
чение пути Ца.) являются постоянными величинами. Следовательно, будет справедливо вынести экспоненты 


с этими значениями за пределы интегралов: 
ЩЕ) | 
ТМ = те ео, | ета | е- 1 М(т)ат . (13) 
0 0 
Примем, что существует некоторое дискретное множество расчетных значений мощности необратимых 
преобразований, полученных в результате анализа дискретной экспериментальной информации о динамике 


процесса резания, то есть №(#) = М, М, №, №4......... Исходя из этого, интегральный оператор (13) примет вид 
следующей суммы: 
Г. 1 Г. я 
ТМ = те ей [м | су | е°-Пам + № | су | °-Пат+ > 
2. й ° . й 1, . - (14) 
> +М; | е* у | с°-Пат+ №4 | ету [ с°Пат+ вы 1] 
Г, 1, Г, ь 


Сделав допущение ЛМ; =0, получим окончательное решение для случая, описываемого выражением 
(14): 
ти - АВА [5 (с) — 9) (ебыь — еб + «> 
о: (15) 
<> +М№;3 (Е) — еб.) еб еб)... ] 
Выражение (15) очень удобно для обработки больших массивов дискретной информации о процессе 
резания. Фактически речь идет о взвешенной относительно конечного значения сумме мощностей необратимых 





преобразований. 

Рассмотрим только один эксперимент из серии, при котором износ по задней поверхности инструмента 
отсутствовал. В следующих случаях менялась изношенность инструмента, а значит, и динамика процесса (как с 
механической, так и с термодинамической точки зрения). Отметим, что полученные в этом эксперименте ре- 
зультаты обобщаются для всех остальных случаев, однако объем статьи не позволяет их рассмотреть. 

Реакция со стороны процесса на формообразующие движения инструмента вдоль оси х представлена 
на рис. 2. 
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Рис. 2. Реакция со стороны процесса на формообразующие движения инструмента вдоль оси х на первом шаге эксперимента 
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Как видно из рис. 2, реакция на формообразующие движения инструмента при токарной обработке но- 
сит существенно нестационарный характер, но при этом наблюдается определенное притягивающее многообра- 
зие, вокруг которого и происходит возмущенное движение инструмента. Аналогично представленному на 
рис. 2 случаю в эксперименте получены реакции на формообразующие движения инструмента вдоль остальных 
осей. Кроме того, доступна информация об ускорениях инструмента относительно обрабатываемого изделия. 
Интегрирование этих данных с учетом известных элементов резания позволило получить значения скоростей 
обработки и пути, пройденного инструментом вдоль обрабатываемой детали. 

В пакете МаНаб была создана специализированная программа для обработки массивов информации, 
описывающих процессы, протекающие при резании (силы реакции, вибрации инструмента и мощность необра- 
тимых преобразований). В этой же программе выполнен расчет температуры в зоне контакта инструмента и 
детали по (15). Параметры модели представлены в табл. 1. 








Таблица 1 
Параметры идентифицированной модели 
1 02 Ат с°С/Нм ©, °С 
0,00099 0.00078 0.000159 24,8 




















На рис. 3 представлены результаты моделирования выражения (15) с учетом параметров табл. 1, а так- 
же экспериментально измеренная характеристика. 
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Рис. 3. Сравнение результатов, полученных в эксперименте и рассчитанных по модели 


Сопоставление рис. 3 и рис. 2 показывает, что измеренное значение температуры при врезании инстру- 
мента в деталь увеличивается с 25 °С (температура окружающей среды) до 330-340 °С. Согласно данным рис. 3, 
сначала наблюдается более быстрый рост измеряемого значения термоЭДС, а в конце процесса термоЭДС ста- 
билизируется. Это объясняется тем, что на стенде 5ТО.201-1 измеряется не сама температура в зоне контакта, а 
термоЭДС естественной термопары, образуемой в контакте инструмента и детали. В начале эксперимента из-за 
резкого скачка температурных градиентов термоЭДС показывает слишком большой и быстрый рост темпера- 
туры, а затем, по мере расхождения температурных градиентов стабилизируется. Этот эффект широко известен. 
Например, в [5] сказано: «Метод естественной термопары ненадежен, так как температура резания стабилизи- 
руется в течение 2-3 с». 

С точки зрения динамики, изменение температуры (рис. 3), полученное в результате моделирования 
выражения (15), близко к процессам, протекающим при резании (рис. 2). На графике виден продолжающийся 
рост температуры, обусловленный продолжающимся ростом сил реакции (рис. 2). Также график отражает все 
динамические особенности процесса, в том числе малоамплитудные скачки сил. Из этого можно сделать сле- 
дующий вывод: предлагаемая модель более точно описывает изменения температуры резания, чем данные 
естественной термопары. Также из рис. 3 ясно, что моделируемая характеристика, описывающая изменение 
температуры в зоне резания, в целом качественно совпадает с экспериментальной зависимостью. Эксперимен- 
тальная и моделируемая характеристики наиболее полно совпадают после стабилизации процесса резания, в 
пределах 15—24 секунд. 
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Обсуждение и заключения. Итоги исследования дают ответ на вопрос, каким образом можно добиться 
адекватного экспериментальным данным описания наследственного характера передачи и накопления тепла в 
процессе обработки металлов точением. С этой целью необходимо смоделировать температуру в зоне контакта 
инструмента и детали на основании предложенного модифицированного оператора Вольтерры. 

С научной точки зрения, наибольший интерес представляет предлагаемая базовая математическая мо- 
дель и результаты ее идентификации на основе эксперимента по реальной обработке металлов. Здесь наблюда- 
ется высокая степень совпадения моделируемых значений тепловой энергии в зоне контакта инструмента и об- 
рабатываемой детали в период установившегося процесса обработки. 
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Введение. Статья посвящена исследованию микротвердости 
и износостойкости алюминиевого сплава (АК5М7) после 
обработки его поверхности электроискровым легированием 
(ЭИЛ) электродами А (Си) и Р (Си +Р) с разной приложен- 
ной энергией. 

Цель работы — исследование взаимосвязи физико- 
механических и прочностных свойств поверхности сплава 
АК5М7, полученной после ЭИЛ медным и медно- 
фосфорным электродами. 

Материалы и методы. Использованы методы рентгеноди- 
фракционного и рентгенофлуоресцентного анализа для 
определения фазового состава, областей когерентного рас- 
сеяния (ОКР) и микродеформации поверхности. 
Результаты исследования. Установлено, что при использо- 
вании обоих видов электродов с увеличением энергии воз- 
растают микротвердость и износостойкость. Образец был 
модифицирован электродом Р с энергией 0,79 Дж. В этом 
случае максимальное значение микротвердости увеличи- 
лось в 5,3 раза, износостойкости — в 1,6 раза. Установлено, 
что при ЭИЛ на поверхности исследуемого сплава образу- 
ются новые интерметаллические фазы А|5Си и СизР. Если 
значение используемой энергии максимально (0,79 Дж), то 
концентрация фазы А]5Си увеличивается в 5 и 9 раз с элек- 
тродами А иР соответственно, а концентрация фазы СизР 
увеличивается более чем в 4 раза при использовании элек- 
тродаР. Полученные экспериментальные данные свиде- 
тельствуют об уменьшении размеров ОКР и увеличении 
значений микродеформации для всех имеющихся на по- 
верхности фаз А|, А Си и Си3Р. 

Обсуждение и заключения. Увеличение энергии при ЭИЛ 
приводит к повышению дефектности и изменению фазово- 
го состава вновь полученных поверхностей. Этим можно 
объяснить повышение износостойкости данной поверхно- 
сти. 


Ключевые слова: электроискровое легирование (ЭИЛ), 
область когерентного рассеяния (ОКР), микродеформация, 
алюминиевые сплавы, интенсивность изнашивания. 
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Введение. Для повышения износостойкости деталей используют различные покрытия и методы моди- 
фикации поверхности. Один из них — электроискровое легирование (ЭИЛ) токопроводящих материалов (авто- 
ры: Б. Р. Лазаренко и Н. И. Лазаренко [1-4]). 

Сплавы А]-$1-Си широко используются в автомобильной промышленности. Они хорошо отливаются, 
обладают всеми необходимыми механическими свойствам, характеризуются низким коэффициентом теплового 
расширения и высоким соотношением прочности и веса. Однако твердость и износостойкость алюминиевых 
сплавов относительно низкие. Этим обусловлена актуальность задачи улучшения названных свойств за счет 
поверхностного упрочнения [5-10]. Речь идет о необходимости создания упрочненных слоев с улучшенными 
физико-механическими свойствами. В представленной работе показаны результаты ЭИЛ различными электро- 
дами при разных энергиях с целью упрочнения поверхности сплава АК5М7. 

Материалы и методы. Для проведения исследований были изготовлены образцы размером 
15х15х4 мм из алюминиевого сплава АК5МУ. Их поверхности обрабатывались методом ЭИЛ медным и медно- 
фосфорным электродами. Элементный состав сплава и электродов определялся методом рентгенофлуоресцент- 
ного анализа на приборе АВГ, Ре{огт?”Х 4200. Результаты приведены в табл. 1. 














Таблица 1 
Элементный состав сплава и электродов 
А, % Си, % 51, % Р, % 
Сплав АК5М7 88 7 5 
Р-электрод 92,9 7,1 
А-электрод 0,1 99,9 























Электроискровое легирование проводилось на установке АМЕК-31 (фирма ЭСПМТТ, Молдова), которая 
позволяет проводить ЭИЛ в широком диапазоне энергий. 

Рентгенофазовый анализ выполнялся с целью определения фазового состава покрытий, нанесенных ме- 
тодом ЭИЛ на алюминиевый сплав. Для съемки использовался рентгеновский дифрактометр РАМа[уйса1 
ЕМРУКЕАМ на СиКо-излучении, укомплектованный никелевым фильтром с автоматической щелью расходи- 
мости первичного пучка. Рентгенограмма расшифровывалась с помощью программы Н1ВЗсогеРав. Изучены: 
— количественное содержание фаз, обнаруженных на поверхности; 

— средний размер областей когерентного рассеяния (ОКР); 
— величины микродеформаций после ЭИЛ образцов. 

Микротвердость модифицированного слоя измерялась на поверхности, полученной после ЭИЛ. В ра- 
боте использовали исследовательский комплекс на базе микротвердомера НУ5-1000 и цифровую видеокамеру. 
Измерения выполнялись по ГОСТ 9450-76 с индентором в виде пирамиды Виккерса (нагрузка Р = 25 г). По 
каждому образцу проводили не менее десяти измерений. Основной сплав и поверхности его образцов после 
обработки ЭИЛ подвергались трибологическим испытаниям (предварительно на поверхности образцов нано- 
сился смазочный материал). 

Результаты исследования. Полученные результаты представлены в табл. 2 и на рис. 1-4. 

Таблица 2 


Концентрация фаз основного сплава и поверхностей, полученных после ЭИЛ в зависимости от энергии 


























Концентрации фаз 
Энергия, Дж Сплав АК5М7 А-электрод Р-электрод 
А] АБСи 51 А] АБСи $1 А] АБСи 51 СизР 
0 92 2 4 
0,07 97 1 94 3 2 1 
0,2 96 3 1 70 22 5 3 
0,39 94 1 65 25 5 5 
0,79 89 10 1 59 28 5 8 
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В табл. 2 приведены концентрации фаз, обнаруженных на поверхности основного сплава и образцов 
после ЭИЛ при различных энергиях. 

В основном сплаве имеются три фазы: А|, АБСи и $1. С повышением энергии ЭИЛ наблюдается увели- 
чение концентрации АБСи при использовании обоих электродов и уменьшение фазы А|, а фаза $1 остается по- 
стоянной для А-электрода при всех энергиях. При этом концентрация 51 возрастает в 2,5 раза с ростом энергии 
для Р-электрода. В образцах, обработанных Р-электродом, появляется фаза СизР. Ее концентрация с повышени- 
ем энергии импульса до 0,79 Дж возрастает в 8 раз. 
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Рис. 1. Зависимость размеров ОКР фазы на основе А! от энергии для электродов АиР 
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Рис. 2. Зависимость величины микродеформации фазы на основе А от энергии для электродов АиР 


Из рис. 1, 2 видно, что при использовании обоих электродов с повышением энергии ЭИЛ уменьшается 
ОКР и увеличиваются микродеформации. Области когерентного рассеяния представляют собой практически 
бездефектные участки поверхности. Следовательно, увеличение ОКР свидетельствует об улучшении микро- 
структуры модифицированного слоя, уменьшении его дефектности. Вместе с тем известно, что изменение де- 
фектности структуры металлических материалов приводит к изменению их физико-механических 
свойств [11, 12]. Величина микродеформации указывает на степень искаженности кристаллической решетки, и 
таким образом по ее величине можно судить о степени дефектности структуры [13, 14]. Аналогичные зависи- 
мости ОКР и микродеформации получены для фаз АБСиа и СизР после ЭИЛ поверхности сплава АК5М7 элек- 
тродами А иР. При этом обнаружено, что закономерности изменения ОКР и микродеформации такие же, как 
для фазы на основе А|. Следовательно, данные рис. 1, 2 позволяют предположить, что с увеличением энергии 
сильнее искажается структура полученной поверхности. 

На рис. 3 представлены результаты измерения микротвердости поверхностей после ЭИЛ медным и 
медно-фосфорным электродами. 
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Рис. 3. Зависимость микротвердости от энергии для А- и Р-электродов 


Видно, что микротвердость возрастает с увеличением энергии, причем для электрода А при макси- 
мальном значении энергии отмечается почти 3-кратный рост, для электрода Р — 5-кратный. 

Известно, что при ЭИЛ происходит плавление материалов анода и поверхности образца, а затем быст- 
рое охлаждение (закалка). В результате на поверхности могут образоваться новые фазы (в нашем случае СизР), 
дислокации, точечные дефекты, высокая степень напряженного состояния поверхности, мелкозернистая струк- 
тура и пр. В связи с этим можно предположить, что с ростом энергии микротвердость увеличивается из-за вы- 
шеперечисленных искажений, возникших в процессе ЭИЛ (см. рис. 3). 

На рис. 4 представлены результаты измерения интенсивности изнашивания сплава АК5М7 и поверхно- 
стей образцов после ЭИЛ медным и медно-фосфорным электродами при энергии 0,79 Дж. 
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Рис. 4. Значения интенсивности изнашивания основного сплава АК5М7 и его поверхности после покрытия при энергии 0,79 Дж 


Из рис. 4 видно, что поверхности образцов, обработанных медно-фосфорным и медным электродами, 
изношены значительно меньше, чем поверхность сплава АК5М7. 

Обсуждение и заключения. Полученные результаты свидетельствуют об эффективности ЭИЛ поверх- 
ности сплава АК5М7 различными электродами и с разной энергией. Поверхности, полученные таким методом, 
изнашиваются меньше. При использовании медного электрода изнашивание сокращается в 1,37 раза, при ис- 
пользовании медно-фосфорного — в 1,57 раз. 

При достижении минимального значения износа после ЭИЛ поверхности имеют максимальные значе- 
ния микротвердости и микродеформации, а также наименьшие значения ОКР. При этом наблюдается значи- 
тельный рост концентрации образующихся после ЭИЛ интерметаллических фаз: А15Си — от 5 до 9 раз с мед- 
ным и медно-фосфорным электродами соответственно; СизР — в 8 раз с медно-фосфорным электродом. При- 
чиной повышения износостойкости поверхностей после ЭИЛ могут быт различные искажения и дефекты, обра- 
зующиеся после плавления, а также формирование новых интерметаллических фаз. 
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Исследование зависимости от частоты констант поляризованной пьезокерамики 
в схемах замещения при слабых электрических полях (часть Ш)” 
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Введение. Принято считать, что электроупругие модули 
не зависят от амплитуды и частоты колебаний. Это под- 
ход отражен в российских и зарубежных стандартах для 
определения полного набора электроупругих модулей 
пьезокерамики. Например, для определения пьезомодуля 
Яз1 образца в форме диска необходимо провести изме- 
рения в трех частотных областях: в области первого и 
второго резонансов, в области антирезонанса и на часто- 
тах много ниже 1 кГц. В связи с этим предполагается, 
что при определении 431 модули исследуемой керамики 
в диапазоне частот от 1 КГц до второго резонанса не 
зависят от частоты. 

Целью данной работы является исследование зависимо- 
сти от частоты электроупругих модулей керамики. При 
этом используется образец в форме диска из ЦТС (цир- 
коната-титаната свинца). 

Материалы и методы. Использованы методы постанов- 
ки и решения задач стационарной электроупругости и 
разделы теоретических основ электротехники. Для реа- 
лизации метода конечных элементов применены метод 
возмущений и пакет прикладных программ АМЗУ5. Ре- 
зультаты экспериментов обработаны в среде МАТГАВ. 
Результаты исследования. Для пьезокерамики ЦТС 
исследованы зависимости от частоты различных моду- 
лей (пьезоэлектрических 431, диэлектрических =1. и 
упругих модулей гибкости $11, 512, 513). Рассмотрены 
радиальные колебания образца в форме диска с электро- 
дами на торцах. Толщина образца — 1 мм, диаметр — 
40 мм, диапазон колебаний — до 700 КГц. Сначала ис- 
следовалась частотная зависимость для упругих модулей 
керамики из определения десяти резонансных частот. 
Затем по измеренным значениям проводимости образца 
была определена зависимость от частоты модулей 
431 И #13. С этой целью использовалось полученное в 
работе выражение для электрической проводимости из 
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решения радиальных колебаний диска с учетом его тол- 
щины. 

Обсуждение и заключения. Разработана методика для 
определения зависимости от частоты модулей пьезоэлек- 
трической керамики ЦТС. Образец в форме диска иссле- 
довался при диапазоне частот 15-650 КГц. Показано, что 
в диапазоне до 650 КГц упругие модули $11, $12, 513 © 
верхним индексом ЕЁ (он опущен) или измеренные при 
постоянном электрическом поле практически не зависят 
от частоты. В указанном диапазоне константы #23, 41, ре 
для рассматриваемых радиальных колебаний имеют не- 
значительную частотную зависимость. 


Ключевые слова: пьезоэлектрическая керамика, элект- 
роупругие модули, схема замещения, зависимость от 
частоты, диск, метод конечных элементов, АМЗУ5, ме- 
тод возмущений. 
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Введение. Значительное число работ посвящено исследованию и развитию математических методов 
решения задач о колебаниях пьезоэлектрических тел. В составе пьезоэлектрических устройств пьезоэлементы 
служат для возбуждения и регистрации колебаний, вызванных внешними воздействиями. Выбор пьезоэлектри- 
ческого материала для измерительных преобразователей и анализ их характеристик требует большого объема 
сведений о параметрах материалов. К таким сведениям относятся: 

— полный набор электроупругих модулей [1], 
— потери, 
— эквивалентные электрические схемы, или схемы замещения пьезоэлементов [2]. 

Реактивные динамические параметры Г, С схем замещения определяются упругими диэлектрическими 
и пьезоэлектрическими модулями, а также плотностью пьезокерамики [3]. Принято считать, что электроупру- 
гие модули — параметры схем замещения постоянны, они не зависят от амплитуды (слабые электрические по- 
ля) и от частоты колебаний. Все это нашло отражение в российских [4] и зарубежных [5, 6] стандартах для 
определения полного набора электроупругих модулей пьезокерамики. Например, для определения пьезомодуля 
Яз1 образца в форме диска необходимо провести измерения в трех частотных областях: в области первого и 
второго резонансов, в области антирезонанса и на частотах много ниже резонансов 1 кГц. В связи с этим пред- 
полагается, что при определении 431 модули исследуемой керамики не зависят от частоты в достаточно широ- 
ком диапазоне от 1 КГц до второго резонанса. 

Целью данной работы является исследование зависимости от частоты электроупругих модулей керами- 
ки. Рассматривается образец в форме диска из ЦТС (цирконата-титаната свинца), т. к. именно такой тип кера- 
мики наиболее хорошо известен. Для проверки корректности представленной методики используется програм- 
ма АМЗУ$ [7]. 

Материалы и методы. Рассмотрим пьезоэлектрический диск толщиной 2 и радиусом а. Введем ци- 
линдрическую систему координат (и, ©, 2), причем ось 2 совпадает с направлением оси поляризации. Коорди- 
натная плоскость 2 = 0 совпадает со срединной плоскостью диска. 

Основываясь на известных линейных пьезоэлектрических соотношениях, уравнениях динамики 
сплошной среды [1] и Максвелла [8], систему уравнений осесимметричных колебаний пьезоэлектрического 
диска можно записать в следующем виде: 


И + роз? = 0, 


.Т,„ + 051, + + ро? =0; д.0, +050, += =0. (р 
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Здесь и далее введены следующие обозначения и определения: И, Й’ — механические смещения по координат- 
НЫМ ОСЯМ Г, 2 соответственно; © — круговая частота; р — плотность; Ти„ — механические напряжения; 2.., О» 
— компоненты вектора электрической индукции; д/ и дз — операторы производных по ги 2. 

В случае осевой поляризации линейные уравнения пьзоэффекта для слабых электрических полей в ци- 
линдрических координатах можно записать [1]: 


Тек = с11 010 + с12 =  с1здз И + ез1д3 Ф, 

Тев = с1204 И + сл - + с1з03 И/ + ез1 03 Ф, 

Т.» = с1з910 + С1з- ой Сзздз И/ + езздз Ф, (2) 
р, = ез1 (010 + =)  езздз И/ — =. ф; О, = е15(0зИ + д.И/) — Е11 д: Ф. 


В соотношениях (2) и далее введены следующие обозначения: сии — упругие постоянные в матричном 
обозначении, измеренные на образцах с закороченными электродами или при постоянном электрическом поле 
Е (верхний индекс Е у констант керамики здесь и далее опущен); еи„ — пьезоконстанты; =и„ — диэлектриче- 
ские проницаемости для зажатого образца; ф — электрический потенциал, причем Ё = —ота4 ф [2, 8]. 

Будем считать, что на электродированных торцах и на боковых поверхностях диска заданы граничные 
условия [1, 2]: 











при 2 ==й Т,„„ = 0; Т,, = 0; ф=2И, (3) 
приг=а Т,., = 0; Т,.. =0; О, =0. (4) 
В (3) величина 27 — подведенная к торцам электрическая разность потенциалов [2, 8]. 
Введем безразмерные координаты и величины по формулам: 


г 2 в р 13 
=-; ЕЕ; ==; О = ов |2; 0. = ха а= 11 = —^ 
5 ые в’ а: Сад’ 11 сз 


Решение краевой задачи (1-3) состоит из суммы двух решений: 
— однородного решения при нулевых граничных условиях при 2 = -Й; 











— частного решения, удовлетворяющего только ненулевым условиям на торцах (3). 

Построенная система однородных решений позволит в дальнейшем удовлетворить и граничным усло- 
виям (4) на боковой поверхности (как правило, с помощью вариационных методов). 

Не составляет труда построить частное решение 02 = соп5{ и 08 = 0, которое автоматически удовле- 
творяет третьему уравнению из системы (1) и граничным условиям на торцах (3). Механические и электриче- 
ские составляющие частного решения согласно первым двум уравнениям системы (1) равны: 


ЗС 
То = ТФ = Або +) изтОО + ен; И° = 0; 1% =0; 
33 
фо = КК + АЙзт(ХО; 02 = —ЕззК; То = еззК + ез.Ах соз(хО; р 
у р у 1 © 


е вы | 
В езз (0500 — КЁ 51100)’ 1 _ 2 ав 00 
Хх 





В (5) и далее введены следующие обозначения: 


Е С33Е с 
с3з 


с3з' 





Если вектор внешних усилий и о потенциал равны нулю на торцевых поверхностях, то 
построение однородных решений связано с определением корней дисперсионного уравнения [9]. Для симмет- 
ричных колебаний дисперсионное уравнение имеет вид: 

а„М, вап 1(Ви) = 0, (и=1,2, 3). (6) 

В (6) введены следующие обозначения: 

ав = а? с1зК1п + СззК2п + @ззКзп; Ви = Киви + Кап; 
М, = Б-Кзз — ВзКз2; М» = БзКз1 — В1Кзз; Мз = В1Кз2 — Б>ЁКза. 
Здесь Кии — алгебраические дополнения элементов третьей строки определителя системы для симметричных 
колебаний (1); В„ — корень бикубического уравнения из [9], и — безразмерное волновое число колебаний по 
ОСИ 7. 
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Для нахождения корней о при заданных значениях © необходимо совместно решить дисперсионное (6) 
и бикубическое уравнения. Подробный анализ корней дисперсионного уравнения симметричных колебаний без 
потерь для пьезоэлектрического слоя приведен в [9], с учетом потерь — в [10]. Именно сумма частного и одно- 
родного решений позволит удовлетворить граничным условиям как на торцах, так и на боковых поверхностях 
диска. 

При произвольном соотношении размеров диска обратная задача о его вынужденных колебаниях (1-4) 
является очень сложной для анализа и имеет конечное аналитическое решение лишь в некоторых частных слу- 
чаях (например, колебания тонкого диска или без учета его толщины, когда # «< 1). Поэтому при определении 
модулей керамики обратную краевую задачу удобнее решать приближенными методами с учетом толщинных 
поправок для низкочастотных радиальных колебаний диска. В этом случае можно получить аналитическое ре- 
шение в виде конечных формул. В данной работе ищем решение в виде разложения по малому параметру Е: 

42 = =202 + 7=* 0% + 1=8 06... Ц =Ц + #2020), + =040, -.... (7) 








Здесь И„ — вектор-функция с составляющими И(0, И’, ф); п — порядок построенной приближенной теории 
рассматриваемых в работе симметричных колебаний диска; 7, | — неизвестные постоянные, которые зависят 
от модулей пьезокерамики, определяются из (1, 2) и удовлетворяют однородным (нулевым) граничным услови- 
ям (3). 

Опускаем громоздкие выкладки. Ограничимся членами с =6 для определения волнового числа а, чле- 
нами с =? для вектор-функции И в (7) и приведем конечный результат рассматриваемой краевой задачи с вы- 
полненными граничными условиями (3) на торцах [11]: 











у= С: , — 912. 10 2431 (в _ 613 _ `` 913. 
Зе: ° 51’ Г 23 (5аа + 512), Сз3 511 +512, 
о г +у] +7: 
=у|— 1е20; 
ПУ 15+ 5. У] +Р 
12 аз1 





а о | 
{е20 = — с. (24:1 +а . 
р 45СР. Е13(1-К2) 14912 11( 31 з3)| 


Здесь 5„„ — модули гибкости при постоянном Е; 431 — пьезомодуль; К» — планарный коэффициент связи; у — 
планарный коэффициент Пуассона; =13 — свободная диэлектрическая проницаемость диска. 

Далее введем следующие определения и обозначения: &5 = =202 — приближенное волновое число ну- 
левого порядка; 42 =Е?202 +1='0% — —  приближенное волновое число второго — порядка; 
02 = =20 7+ =“ + 1=б Од ° — приближенное волновое число четвертого порядка; С — емкость диска. 

В табл. 1 представлены результаты точного решения волнового числа о керамики РИТА в зависимости 








от частоты из дисперсионного уравнения (6) и приближенный расчет из (7) при # = 0,033 для различных при- 
ближений о. 

















Пример расчета волнового числа а при различных частотах для рассматриваемого а —— 
КГц о а, из (6) 60 02 94 с р1его = 0 04 
50 0,085 0,045055 0,045049 0,045055 0,045055 0,045055 
250 0,425 0,2262 0,2252 0,22603 0,22605 0,2261 
500 0,85 0,4578 0,4505 0,4568 0,4574 0,4575 
700 1,1903 0,6522 0,6307 0,64777 0,6514 0,6518 





























Приведенные в табл. | результаты показывают, что для диска толщиной | мм и частот до 700 КГц 
расчет волновых чисел из дисперсионного уравнения (6) и приближенный расчет для 04 практически совпада- 
ют, а поправкой рЁего можно пренебречь. В этом случае разложение (7) для а? практически зависит только от 
геометрических размеров диска, плотности и модулей гибкости $11, 512, 513. 

Опуская относительно громоздкие выкладки, приведем выражение проводимости У для пьезоэлектри- 


ческого диска. Оно получено из приближенного решения (7) для различных =. 
7000 — нулевое приближение # = 0, или известное уравнение радиальных колебаний диска нулевой толщины: 


зы 96: оббы Воее 
Е. В и хо (х) — 1-1 (х)/ 


7040 — приближенное волновое число 4-го порядка: 
ча. _ И. ло 
х=-; У040 = «С (1 Кр + Кр (1 + у) ны) (8) 


7042 — приближенное волновое число четвертого порядка, частное решение второго порядка: 
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02е2. х 
нова = с (1 Вена Л ) 


З= 333 Сзз О Ло(х) — (1 — у)/1 (<) 


7242 — приближенное волновое число четвертого порядка, однородное и частное решение второго порядка: 





02 е2. И ) т @ззС1з 
узна = ес (1 Чена + Ка + У) (1 тих 2 ›) 
2зз с3з е31С3з 


1: 
эт = Боб — а - Уре [1 + 208 (5+ : = аду хо 69; (9) 


ЕззС1з _ 431 
езз 511 + 512 
В табл. 2 представлены результаты численно-аналитических расчетов. Проводимости точного реше- 





Р _ 
ез1 = ез. + 


ния (5, 6) и приближенных решений (8, 9) по =2 найдены с помощью программы МАТЕАВ [12]. Для диска из 
РИТА [13] толщиной | мм и диаметром 30мм численные расчеты проводились в системах АМЗУЗ и 
АСЕГАМ [14]. Последний пакет программ разработан сотрудниками Южного федерального университета 
(ЮФУ) и ориентирован на расчеты пьезоэлектрических устройств. 














Таблица 2 
Расчет проводимостей для различных частот 
Частота, КГц |АМЗУ5, См/м |Решение (5,6), См/м | 7000, См/м | 7040, См/м |042, См/м | 7242, См/м 
350 5,03 5,032 4,9688 5,014 5,0457 5,028 
554,4 8,74169 8,741029 8,478 8,5864 8,73237 8,735 
700 10,98 10,97 10,376 10,682 10,957 10,96 





























Из табл. 2 видно, что при построении приближенных теорий типа (8, 9) желательно использовать раз- 
ложение (7) для расчета волнового числа как минимум второго порядка. 

Результаты исследования 

1. Исследование зависимости упругих модулей керамики от частоты. Для достаточно хорошо изу- 
ченной пьезокерамики ЦТС проведем исследование зависимости от частоты следующих модулей: пьезоэлек- 
трических 431, диэлектрических #13 и модулей гибкости $11, $12, 51з. Рассмотрим низкочастотные радиальные 
колебания образца в форме диска с электродами на торцах. Толщина образца — 1 мм, диаметр — 40 мм, диапа- 
зон колебаний — до 700 КГц. Сначала исследуем частотную зависимость для упругих модулей. Для этого, со- 
гласно [15], предварительно измеряем три первые резонансные низкие частоты {. (основной резонанс и два его 
обертона). Упругие константы определяют из решения частотного уравнения (9) радиальных колебаний диска с 
учетом относительной толщины : три уравнения для трех неизвестных Именно относительная толщина # отли- 
чает приведенное выше частотное уравнение (9) для радиальных колебаний диска конечной толщины от из- 
вестного частотного уравнения для радиальных колебаний диска с «нулевой толщиной» [15]: 

В.Л(®) = 1-У) л(®. 

Введение толщинной поправки & в решение известного уравнения радиальных колебаний образца в ви- 
де диска повышает точность и информативность измерения упругих модулей. 

Определенные по [11, 15] методом трех резонансов упругие податливости для рассматриваемой кера- 
мики оказались равны: 

511 = 12,29е- 12, 512 = -4,05е- 12, 53 =—5,28е - 12. 
В табл. 3 приведены первые десять резонансных частот для керамики ЦТС19. В стоке «Эксперимент» 





указаны частоты, измеренные на установке для измерения проводимости И’К 6500В в Институте высоких тех- 
нологий и пьезотехники при ЮФУ. В строке «Расчет» — рассчитанные по формуле (9) для независимых от ча- 
стоты упругих модулей $511, 512, 51з, определенных по [11, 15]. Погрешности измерений не превышали значений, 
рекомендуемых стандартом [4]. 











Таблица 3 
Первые десять резонансных частот для керамики ЦТС19 
Резонансные частоты /,, КГц 
Эксперимент 50,8 132,8 |211 2884 |364,56 |4394 |513,2 |5854 |6544 |722,8 
Расчет 50,8 132,8 |210,95 |288,1 |364,8 |440 |516 590,4 |663 735,4 









































Из табл. 3 видно, что в диапазоне частот до 600 КГц различие между расчетными и экспериментальны- 
ми данными при постоянных упругих модулях не превышает 1 %, поэтому можно сделать следующий вывод: 
«Для рассматриваемой керамики упругие модули $11, $12, $1з с верхним индексом Ё (он опущен) или измерен- 
ные при постоянном электрическом поле практически не зависят от частоты в диапазоне до 650 КГц». 
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2. Исследование зависимости от частоты модулей 431 и #213. Используем выражение (8) для прово- 
димости У040 из решения радиальных колебаний диска с учетом его толщины, а также измеренные при ком- 
натной температуре значения мнимых частей комплексной проводимости исследуемой пьезокерамики. В этом 
случае можно исследовать зависимость модулей 431 и =13 от частоты в диапазоне от 10 КГц до 600 КГц. Для 
определения двух неизвестных 41 и #1. измерим У при двух значениях частот /, №. Разница между № и Л вы- 
брана так, что можно пренебречь зависимостью модулей 41 и =1з от частоты в диапазоне /- р. 

В настоящей работе введены два условия: 

 Б-Лл=200 Гц; 

2) проводимости рассчитаны для идеального пьезоэлектрического тела или без учета потерь, поэтому 
не учитываются частоты резонансов или в окрестности резонансов. 

В итоге получим легко решаемую систему двух линейных уравнений относительно двух неизвестных 

2: Г 
Кри 233. 

На рис. 1 представлены зависимости К ‚Чи 213 /Е0 от частоты 10-650 КГц. 
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Рис. 1. Зависимости пьезоэлектрических параметров от частоты 


Здесь & — диэлектрическая проницаемость вакуума. В диапазоне от 15 до 650 КГц коэффициент связи 
| сначала увеличивается с ростом частоты от 0,28 до 0,34, а затем убывает до 0,26. Пьезомодуль 431 имеет ана- 
логичную частотную зависимость: при 15 КГц с —1,50е '° до максимума 1,661е '°, а затем убывает до 1,31е ®. Отно- 
сительная диэлектрическая проницаемость с ростом частоты монотонно уменьшается с 1766 до 1455. 

На рис. 2 представлены зависимости от частоты в диапазоне 1-10 КГц измеренной и рассчитанной по 
формуле (9) емкости керамики при постоянных модулях, определенных согласно ОСТ [4, 15]. 
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Рис. 2. Зависимость емкости пьезокерамики от частоты 


Мадорский В. В. и др. Исследование зависимости от частоты констант поляризованной пьезокерамики в схемах замещения 
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Видно резкое отличие двух емкостей керамики ЦТС для низких частот. Вероятно, этим и объясняется 
зависимость параметров керамики от частоты ниже 15 КГц (см. рис. 1). Это предполагает наличие более высо- 
кой степени дисперсности — в первую очередь, у пьезомодуля 431. Подробнее этот вопрос будет рассмотрен в 
следующей статье (константы #зз^Т, 41, А» исследуются в низкочастотном диапазоне). 

Обсуждение и заключения. Разработана методика для определения зависимости от частоты модулей 
пьезоэлектрической керамики ЦТС. Образец в форме диска исследовался при диапазоне частот 15-700 КГц. 
Показано, что в диапазоне до 650 КГц упругие модули $11, 512, 51з с верхним индексом Е (он опущен) или изме- 
ренные при постоянном электрическом поле практически не зависят от частоты. В указанном диапазоне кон- 
станты 233, 41, К» для рассматриваемых радиальных колебаний имеют незначительную частотную зависимость. 

Изучена зависимость только реальной части модулей керамики от частоты, потери не учитывались. По- 
этому экспериментально полученные зависимости мнимой части проводимостей керамики от частоты были 
измерены вдали от резонансов, где влияние потерь ничтожно. Тем не менее в дальнейшем необходимо учиты- 
вать как реальные, так и мнимые части модулей. Это значит, что потери являются включенными. Данную про- 
блему предполагается изучить с учетом зависимости комплексных модулей исследуемой керамики от частоты. 

В представленной работе рассмотрены радиальные колебания образца для получения более полной ин- 
формации при измерении большего набора констант (кроме пьезоэлектрических и диэлектрических модулей 
измерены упругие модули $11, 512, 51з ). Данный метод может быть распространен и на другие формы образцов 
(стержни, пластины и т. п.), поскольку в этих одномерных задачах происходит простая замена функций Бесселя 
на тригонометрические функции. 
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Сравнительная характеристика пространственно-затрудненных фенолов и 


нитроксильных радикалов как стабилизаторов фотодеструкции полиэтилена" 


И. Ю. Жукова', И. И. Кашпаров?, С. В. Кучеренко*, Е. Ш. Каган“”” 


1,2,3 


Донской государственный технический университет, г. Ростов-на-Дону, Российская Федерация 


24 Южно-Российский государственный политехнический университет (НПИ) им. М. И. Платова, г. Новочеркасск, Российская Федерация 


Сотрагайуе сВагасет1$Яс$ о? емсаПу Бтдегеа рвепо[5 апд пйгох14е 


га@!са[5 а; баб Шхег$ оЁ роГуе У епе рвофо4есгадайоп ыы 


1. Уи. Ивикоуа', 1. 1. Казпрагоу?, $. У. Киспегепко?, Е. 5В. Казап“”" 


123 роп аще Тесыса! Ошуегзйу, Возюу-оп-Ооп, Виззап Еедеганоп 


24 Рау Зои -Взз1ап Зе Ро[убесвис ОтуегзНу (МР, МоуосвегказзК, Визз1ап ЕедегаНоп 


Введение. Полимерные материалы (ПМ) все шире исполь- 
зуются в различных отраслях промышленности и сельско- 
го хозяйства. Под действием УФ-света ПМ разрушаются. 
Для защиты ПМ от фотодеструкции применяются УФ- 
стабилизаторы. Их действие основано на поглощении 
фотохимически активной компоненты солнечного света 
или на деактивации возбужденных молекул, уже погло- 
тивших квант света, а также на торможении темновых 
реакций, индуцированных светом. Цель работы — срав- 
нительный анализ соединений из ряда пространственно- 
затрудненных фенолов (ПЗФ) и нитроксильных радикалов 
(НР) в качестве ингибиторов фотодеструкции ПМ. 
Материалы и методы. Для исследования использовали 
реактивы квалификации «ч», полиэтиленовую пленку 
ПЭВД в/с 15803-020. Стабилизатор наносили на пленки 
методом окунания. Карбонильные группы в полиэтилене 
определяли методом ИК-спектроскопии. ИК-спектры ре- 
гистрировали на приборе Уапап-640. 

Результаты исследования. В качестве стабилизаторов 
фотоокислительной деструкции полиэтилена испытаны 
2,4,6-три-трет-амилфенола (1), 
бутилфенола (2), 
тетраметилпиперидин-1-оксила (3), 


2-метил-4,6-ди-трет- 
4-ацетиламино-2,2,6,6- 
З-карбоксамидо- 
2,2,5,5-тетраметилпирролин-1-оксила (4). 
Известно, что при облучении в образцах ПМ в процессе 
фотодеструкции методом ИК-спектроскопии фиксируется 
накопление карбонильных и гидроксильных групп. В ИК- 
спектрах окисленного полиэтилена появляется полоса 
поглощения карбонильной группой при 1720 см '. Анализ 
ИК-спектров показывает, что содержание карбонильных 
групп в контрольных образцах значительно выше, чем в 
образцах, обработанных растворами стабилизаторов. 


* о 
Работа выполнена в рамках инициативной НИР. 


Ттгодиспоп. Ро]утейс таепа]з (РМ) аге шсгеазту изе4 
ш уапоиз пез апа ахиси те. Оп4ег Фе асноп оЁ ЧУ 
По, РМ аге дезкоуе4. ДУ заб Штегз аге и5е4 №0 ргоесЕ РМ 
Нот рБоюдезтаданоп. ТБеш асНоп 15 Базе оп Фе абзогрНоп 
оё Фе рБоюасйуе зипИ2{ сотропеп{ ог оп Фе 4еаснуаНоп 
оф ехсИе4 по]есшез Вай Вауе атеа4у абзогбе4 а ИР диап- 
сбит, аз ме аз оп Фе шоп оЁ дак Пзоб-штаисе4 геас- 
@01$. ТВе \уогК оБеснуе 1$ №0 ргоу1е а сотрагануе апа[уз1$ 
ОР сотроипа$ оЁ Фе зепсаПу Блт4ете4 рКепо]5 (НР) зепез 
апа пигох14е га@ са! (МВ) аз Ше РМ рпоюдезтаданоп шЫЪ- 
Пот. 

Маепта[5 апа Мешо4. Веасепз ог е “ригит” отаде, 
ГОРЕ ргетпит стаде 15803-020 роуевйУепе Ят \уеге изе4 
Юг Фе шуезНганопз. ТВе за БШтег ууаз аррПе4 ю фе НИ Ъу 
Ч1рри®. Сагфопу! эгопрз ш ро]уефУепе \еге де{егишед Бу 
18 зресвозсору. 18 зресфа \уеге тесог4е4 оп Уапап-640 


штягатеп(. 
Кезеатсй Кезийз. 2,4,6-и1-е-атурБепо] (1), 2-теу!-4,6- 
1-4ег-БибурВепо] (2), 4-асе{УЛатито-2,2,6,6- 


3-сагбохап19о-2,2,5,5- 
{енатефу!руггоПп-1-оху1 (4) \уеге 1е5еЯ аз за БИтег$ Юг 


{енатефу!ррепат-1-оху! (3), 


р®Вою-ох14айуе 4еотаданоп оЁ ро!уеУепе. [ 1$ Кпо\уп а 
фе асситшаНоп оЁ сагфопу| ап4 ВудгохУ| этопре 1$ гесотде4 
ш РМ затр/ез ипдег птаФаноп ш Фе ргосез$ оЁ рБою4езта- 
аНоп гой Фе Г зресёгозсору. ТВе абзогрНоп Бапа оЁ 
Фе саопу! эгоир аё 1720 сшг ' арреагз п Ше 1В. зресфа оЁ 
ох11те ро]уепу[епе. Тре Г зресга апа1уз1$ зНо\з аё фе 
сотепё оЁ сагропу! этоирз ш Ше спесК запр[ез 1$ $1етИЙсап{- 
[у ЫзВег Фап ш Фе запар]ез {геаёе4 Бу заб Шхег зоаНопз. 


^* Е-птай: {ул @тай.та, {уап-Казбрагоу@уапдех.ги, Кив-зу82@тай.та, Казап2 9 @тай.га 


** Тве \огК Уаз регогте4 ш пе Нате\могК оЁ фе шшайуе В & О. 
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Обсуждение и заключения. Эксперименты показали, что 
лучшими фотостабилизаторами полиэтилена являются 
З-карбоксамидо-2,2,5,5- 
4-ацетиламино-2,2,6,6- 


нитроксильные 
тетраметилпирролин- 1-оксил и 


радикалы 


тетраметилпиперидин-1-оксил. Причем нет существенной 
разницы между 2,2,6,6- 
тетраметилпиперидина и 2,2,5,5-тетраметилпирролина. 
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Введение. Полимерные материалы (ПМ) все шире используются в различных отраслях промышленно- 
сти и сельского хозяйства. При достаточно высокой устойчивости к воздействию климатических факторов ра- 
бочие характеристики ПМ снижаются под влиянием влажности и температуры, а также ультрафиолетовой (УФ) 
составляющей спектра солнечного излучения: в этом случае в полимере образуются радикальные частицы и 
происходят трансформации, которые ведут к разрушению материала [1, 2]. Если в ПМ есть продукты его окис- 
ления, например кетоны, то они являются фотоинициаторами процесса разложения полимера. Для защиты ПМ 
необходимо обеспечить взаимодействие радикалов с высокоэффективными светостабилизаторами [3]. С этой 
целью применяются УФ-стабилизаторы, которые могут действовать следующим образом: 

— поглощают фотохимически активные компоненты солнечного света; 

— деактивируют возбужденные молекулы, уже поглотившие квант света; 

— тормозят темновые реакции, индуцированные светом; 

— разрушают или деактивируют фотохимически активные примеси и продукты фотореакции [4, 5]. 

В качестве УФ-стабилизаторов полимеров широко применяются пространственно-затрудненные фено- 
лы (ПЗФ), пространственно-затрудненные амины (ПЗА) и нитроксильные радикалы (НР) [6, 7]. ПЗФ действуют 
как агенты обрыва растущей цепи. Поскольку присоединение подвижного атома водорода к первичным ради- 
калам фотолиза полимеров является окислительной реакцией, сами стабилизаторы при этом легко окисляют- 
ся [3]. ПЗА при фотодеструкции образуют стабильные НР, взаимодействующие с алкильными радикалами по- 
лимеров, что приводит к обрыву цепи фотодеструкции [2]. 

Целью работы является сравнительный анализ соединений из ряда ПЗФ и НР в качестве ингибиторов 
фотодеструкции ПМ. 

Материалы и методы. Для исследования использовали реактивы квалификации «ч» фирмы АШпсв, 
пленку из полиэтилена высокого давления (ПЭВД 15803-020, ГОСТ 16.337-77) толщиной 0,2 мм. Стабилизато- 
ры фотоокислительной деструкции соединения: 

— из ряда НР — 3-карбоксамидо-2,2,5,5-тетраметилпирролин-1-оксил (1), 
тетраметилпиперидин-1-оксил (2); 
— из ряда ПЗФ — 2-метил-4,6-ди-трет-бутилфенол (3) и 2,4,6-три-трет-амил-фенол (4). 

Подготовку поверхности пленки для нанесения растворов стабилизаторов проводили следующим обра- 
зом: образцы размером 5х5 см закрепляли в кассете, не повреждая рабочую часть образца, промывали в теплой 
воде, промокали хлопчатобумажной тканью, обезжиривали в ацетоне и в течение трех минут сушили на возду- 
хе при температуре 20-25 °С. 

Для приготовления растворов стабилизаторов в 12 мерных колб емкостью 100 мл наливали по 50 мл 
толуола и добавляли стабилизаторы (1-4) по 0,01; 0,005; 0,0025 моль. Затем содержимое колб доводили до мет- 
ки растворителем. Полученные растворы использовали для обработки образцов пленки. 

Подготовленную пленку на 30-40 мин окунали в стакан с раствором УФ-стабилизатора [8, 9] нужной 
концентрации. Затем образец извлекали из раствора и 3—5 мин выдерживали над стаканом, чтобы избыток рас- 
твора стек с обработанной поверхности. Образец полностью высушивали в подвешенном состоянии. 

Для каждой концентрации раствора стабилизатора испытывали по 5 образцов пленки. В целом было 
испытано 120 образцов: 60 штук — в естественных условиях и столько же — в искусственно созданных. 
10 контрольных образцов (КО) не обрабатывались стабилизаторами. 

После проведения испытаний были сняты ИК-спектры контрольных и экспериментальных образцов. 
ИК-спектроскопия определила содержание в образцах карбонильных групп [10, 11]. ИК-спектры регистрирова- 
ли на приборе Уапап 640. 
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Результаты исследования. Наиболее достоверным методом определения стойкости ПМ к воздействию 
УФ являются натурные испытания. Для сокращения сроков эксперимента применяют ускоренные испытания в 
лабораторных условиях. Образцы подвергали фотодеструкции с УФ-стабилизатором (соединения (1-4)) и без 
него. При этом на образцы воздействовали: 
— в течение 365 дней естественный источник УФ-лучей — солнце (натурные испытания согласно ГОСТ 9.708- 
83); 
— в течение 504 часов искусственный источник УФ-лучей — ксеноновая лампа (плотность потока УФ- 
излучения 68 Вт/м* в интервале длин волн 300—400 нм, ГОСТ 9.708-83). 

Известно, что в процессе фотодеструкции методом ИК-спектроскопии в образцах ПМ фиксируется 
накопление карбонильных и гидроксильных групп [1, 12]. Но выход гидроксильных групп в несколько десятков 
раз меньше, чем карбонильных [12], поэтому скорость и степень деструкции полимера определяется концен- 
трацией карбонильных групп по соответствующим полосам поглощения в ИК-спектре [10]. 

Аналитическое приложение ИК-спектроскопии к определению степени окисления пленки основано на 
различии спектров обработанного и не обработанного (контрольного) образцов. После УФ-облучения в ИК- 
спектрах полиэтилена (ПЭ) появляется полоса поглощения карбонильной группы при 1720 см ' (валентные ко- 
лебания карбонильной группы С = О). Эта полоса является аналитической для проведения количественного 
анализа. В максимуме полосы поглощения при 1720 см ' методом базовой линии определяли значение оптиче- 
ской плотности. Расчеты содержания карбонильных групп в образцах выполнены по известной методике [10]. 
Результаты исследований представлены в табл. | и 2. 

Таблица 1 
Содержание карбонильных групп в образцах полиэтилена”, обработанных растворами стабилизаторов, 
при облучении искусственным источником УФ-лучей — ксеноновой лампой 






























































Концентрация раствора 
стабилизатора, 
моль/Л 01 М 0,05 М 0,025 М 
ь. Наименование 
п/п 
стабилизатора 
Содержание С = О групп, % 
Й З-карбоксамидо-2,2,5,5-тетраметил-пирролин-1- __ __ 0.02 
оксил (1) ? 
2 4-ацетиламино-2,2,6,6-тетраметил-пиперидин-1- __ __ 0.025 
оксил (2) 
3 2-метил-4,6-ди-трет-бутилфенол (3) 0,10 0,25 0,40 
2,4,6-три-трет-амилфенол (4) 0,10 0,30 0,45 
5 Контрольный образец” 2,50 2,50 2,50 
*До испытаний карбонильные группы в образце не обнаружены. 
**Не обрабатывался раствором стабилизатора. 
Таблица 2 
Содержание карбонильных групп в образцах полиэтилена", обработанных растворами стабилизаторов, 
при облучении естественным источником УФ-лучей — солнцем (период облучения 365 дней) 
Концентрация раствора 
стабилизатора, 
№ моль/Л 0,1 М 0,05 М 0,025 М 
_ | Наименование 
пт 
стабилизатора 
Содержание С = О групп, % 
1 З-карбоксамидо-2,2,5,5-тетраметилпирролин-1- __ = 0.016 
оксил (1) ? 
2 4-ацетиламино-2,2,6,6-тетраметилпиперидин-1- __ = 0.021 
оксил (2) | 
3 2-метил-4,6-ди-трет-бутилфенол (3) 0,07 0,095 0,16 
2,4,6-три-трет-амилфенол (4) 0,06 0,09 0,17 
5 Контрольный образец” 0,35 0,35 0,35 




















*До испытаний карбонильные группы в образце не обнаружены. 
**Не обрабатывался раствором стабилизатора. 
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Анализ ИК-спектров показал, что содержание карбонильных групп в контрольных образцах значитель- 
но выше, чем в обработанных. На рис. 1 приведены ИК-спектры КО и образцов, обработанных 0,025 М раство- 
рами стабилизаторов (2) и (4) после облучения искусственным источником УФ-лучей. 
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в) 
Рис. 1. ИК-спектры образцов ПЭ-пленки: обработанной 0,025 М раствором 4-ацетиламино-2,2,6,6-тетраметилпиперидин-1- 
оксила (2) (а); обработанной 0,025 М раствором 2,4,6-три-трет-амилфенола (4) (6); контрольный образец (в). 
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Начало процесса деструкции в образце, обработанном ПЗФ (4), характеризуется появлением в ИК- 
спектре полосы поглощения при 1720 см ' незначительной интенсивности (рис. 1, 6). Такая полоса практически 
отсутствует в образце, обработанном стабилизатором (2) из ряда НР (рис. 1, а). Следовательно, вклад НР в тор- 
можение процессов фотодеструкции значительно выше, чем ПЗФ. 

УФ-излучение разрушило КО. Об этом свидетельствует эффект накопления карбонильных групп в об- 
разце, выражающийся в увеличении интенсивности полосы поглощения при 1720 см ' (рис. 1, в). 

С помощью ИК-спектроскопии показано, что при УФ-облучении ПЭ-пленок происходит накопление 
карбонильных групп, связанное с фотодеструкцией углеродных цепей ПЭ. На рис. 2 представлены кинетиче- 
ские кривые накопления карбонильных групп при искусственном фотоокислении пленок в зависимости от вре- 
мени облучения. 
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Рис. 2. Зависимость скорости накопления карбонильных групп от времени искусственного УФ-облучения образцов, 
обработанных: 1 — 0,025 М раствором 3-карбоксамидо-2,2,5,5-тетраметилпирролин-1-оксила; 2 — 0,025 М раствором 
4-ацетиламино-2,2,6,6-тетраметилпиперидин-1-оксила; 3 — 0,025 М раствором 2-метил-4,6-ди-трет-бутилфенола; 

4 — 0,025 М раствором 2,4,6-три-трет-амилфенола; 5 — контрольный образец не обрабатывался 


Анализ полученных данных показывает, что НР (1, 2) и ПЗФ (3, 4) существенно замедляют фотоиници- 
ированную окислительную деструкцию ПЭ и уменьшают количество разрывов макромолекул. Об этом свиде- 
тельствует появление индукционного периода с увеличением времени экспозиции (см. рис. 2). 

На начальном этапе эксперимента отмечаются изменения в КО. Постепенно растет концентрация кар- 
бонильных групп в пленке, не обработанной стабилизатором. Это связано с увеличением скорости окислитель- 
ных процессов, ведущих к разрушению структуры углеродных цепей ПЭ и образованию низкомолекулярных 
активных компонентов (радикальных частиц). В образцах ПЭ, обработанных стабилизаторами (1-4), изменений 
не наблюдается. Карбонильные группы фиксируются примерно через 150 часов в образцах, обработанных рас- 
творами ПЗФ (3, 4); через 220 часов — в образцах, обработанных растворами НР (1, 2). Скорость накопления 
карбонильных групп линейно растет в течение 150—200 часов. Далее зависимость сохраняется в течение всего 
процесса фотоокисления (до 504 часов). 

Индукционный период в присутствии НР и ПЗФ может быть объяснен способностью стабилизаторов 
(1-4) с первого момента количественно улавливать радикалы, в результате на этом этапе процесс деструкции 
ингибируется. При этом НР быстро взаимодействуют с алкильными радикалами, при участии полимера пре- 
вращаются в гидроксиламины, способные восстанавливать радикальные частицы [6, 7, 12]. ПЗФ могут превра- 
щаться в хиноны или моно- и диэфиры гидрохинонов [6], которые также являются ловушками для радикальных 
частиц. 

Образцы, обработанные и не обработанные стабилизаторами, облучались также естественным источ- 
ником УФ-лучей согласно ГОСТ 9.708-83 в течение 365 дней (см. табл. 2). Контроль показателя естественного 
фотоокисления ПЭ пленок проводился через 1, 3, 6, 9 и 12 месяцев (рис. 3). 
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Рис. 3. Контроль показателя концентрации карбонильных групп в образцах ПЭ-пленки, обработанной 0,025 М раствором 
стабилизатора (1—4) и контрольного образца (5) 


Сравнение образцов, обработанных НР (1, 2) ПЗФ (3, 4), позволило отметить следующие результаты: в 
первом случае концентрация карбонильных групп повышается несущественно, во втором случае содержание 
карбонильных групп на порядок выше. В КО значительный рост концентрации карбонильных групп фиксиро- 
вался через 1 месяц. 

Следует также отметить, что в процессе испытаний КО ПЭ-пленки становились хрупкими, а по окон- 
чании эксперимента разрушались, в отличие от образцов, обработанных растворами стабилизаторов. 

Обсуждение и заключения. Установлено, что лучшими фотостабилизаторами полиэтилена являются 
нитроксильные радикалы. Не обнаружена существенная разница между радикалами ряда 2,2,6,6- 
тетраметилпиперидина и 2,2,5,5-тетраметилпирролина. ПЗФ в условиях фотодеструкции оказывают меньшее 
стабилизирующее влияние, уступая нитроксильным радикалам. Ингибирующее действие НР и ПЗФ связано с 
их способностью взаимодействовать с радикальными частицами, образующимися в результате фотодеструкции 
полимера. 

Испытанные соединения из класса НР эффективно замедляют процесс фотоокислительной деструкции 
и могут быть рекомендованы для использования в качестве стабилизаторов полимеров. 
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Введение. Статья посвящена исследованию синтезирован- 
ной принципиальной схемы фрикционного контакта твер- 
дых тел в кузнечно-штамповочных машинах. Установлена 
возможность получения максимума нагрузочной характе- 
ристики фрикционного контакта внутри интервала изме- 
нения коэффициента трения. Выявлены две следующие 
возможности сил трения фрикционного контакта: на гра- 
ницах указанного интервала они будут равны при наличии 
максимума равенства; при данных условиях они достига- 
ют наибольшей стабильности. 

Материалы и методы. При изменении величины угла 
меняется положение точки максимума. Это приводит к 
нарушению равенства сил трения на границах интервала 
изменения коэффициента трения. В таком случае коэффи- 
циент точности должен определяться отношением макси- 
мума функции к наименьшему граничному значению. Для 
этого необходимо установить тенденции изменения гра- 
ничных значений функции, связанные с варьированием 
величины угла. Для решения этой задачи новую величину 
тангенса угла давления представили в виде произведения 
коэффициента варьирования на базовое значение тангенса 
угла. 

Результаты исследования. Полученные результаты пока- 
зывают высокую стабильность силы трения при проскаль- 
зывании тел фрикционного контакта (ФК). Однако при 
больших величинах угла давления чувствительных эле- 
ментов датчика-преобразователя максимальная сила тре- 
ния кратковременно может быть пропорциональна теку- 
щему значению коэффициента трения. 

Обсуждение и заключения. Модернизированная принци- 
пиальная схема ФК позволяет теоретически получить 
очень высокую стабильность силы трения. ФК не должен 
обращаться в ноль в интервале изменения коэффициента 
трения выходного параметра основной фрикционной 
группы (ОФГ) и при наличии максимума функции нагру- 
зочной способности ФК. Необходимым условием этого 
является передача чувствительными элементами дополни- 
тельной фрикционной группы (ДФГ) ее полной нагрузки. 
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Введение. Исследована синтезированная принципиальная схема фрикционного контакта (ФК) твердых 
тел в кузнечно-штамповочных машинах. В результате анализа установлена возможность получения максимума 
нагрузочной характеристики ФК внутри интервала изменения коэффициента трения. Выявлены две следующие 
возможности сил трения фрикционного контакта: на границах указанного интервала они будут равны при 
наличии максимума равенства; при данных условиях они достигают наибольшей стабильности. 

Полученные результаты показывают высокую стабильность силы трения при проскальзывании тел ФК. 
Однако при больших величинах угла давления чувствительных элементов датчика-преобразователя максималь- 
ная сила трения кратковременно может быть пропорциональна текущему значению коэффициента трения. 

Материалы и методы. Указанного недостатка лишена схема ФК, приведенная на рис. 1. 








Рис. 1. Принципиальная схема фрикционного контакта 


Основную фрикционную группу ГФ1 составляют тела 1, 2 и 4; дополнительную фрикционную группу 
ГФ2 — тела 3, 5 и 10. Между телами 1 и 2, 2 и 3 размещены в профилированных гнездах чувствительные эле- 
менты в виде тел качения 7 и 8. 

Отличием рассматриваемой схемы является разделение тела 2 на две части, благодаря чему связь меж- 
ду телами 2 и 5 осуществляется посредством тела 10, взаимодействующего с телом 2 через пружины 11. Это 
позволяет передавать половину общей нагрузки ГФ2 от тела 2 непосредственно на тела качения 8, сделав их 
ведущим элементом в составе ГФ2. Вторая половина общей нагрузки ГФ2 передается за счет пружин 11 от тела 
2 телу 10 и далее посредством трения — телу 5. Кроме того, место приложения движущей силы Ё)„ перенесено 


с тела | на тела качения 7. Указанные отличия позволяют исключить влияние величины угла давления тел ка- 
чения 7 на распределение нагрузки между поверхностями трения в ГФ1 и ГФ2 при изменении величины коэф- 
фициента трения [1]. 

Для сведения к линейному характеру трения между телами | и опорной поверхностью 9, а также между 
телами 2 и 10 установлены катки 6 и 12. 

Формула для определения величины силы трения ФК в зависимости от возмущающего воздействия 
имеет следующий вид (при одинаковых значениях коэффициента усиления обратной связи ГФ1 и ГФ2): 


4Е 
я. = = , (1) 
(Иво) 
где УЕ, — суммарная предельная сила трения между телами 1, 2, 4 и 3, 5, 10; ЁР, — начальное усилие замы- 


кания фрикционных пар; Г — текущее значение коэффициента трения между упомянутыми телами; а — угол 
давления между телами качения 7, 8 и гнездом. 

Установлено также, что функция (1) внутри интервала значений Гш...Глах Не имеет максимума, так 
как сила трения ГФ1 обращается в нуль при значении коэффициента трения /, =1/ 19. (Здесь Риш и Лиах — 


соответственно наименьшее и наибольшее значения, которые может принимать коэффициент трения в реаль- 
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ных условиях эксплуатации ФК для принятого сочетания материалов фрикционных пар.) Это же значение со- 
ответствует мнимому максимуму функции (1), поскольку в интервале Х....Г пах нагрузка ФК передается фрик- 
ционной группой ГФ2, сила трения которой возрастает с увеличением /. Это не обеспечивает стабильность 
силы трения ФК, несмотря на то, что она несколько выше, чем у ФК первого поколения (при значениях 
ль = 0,1, Ллах = 0,8 и 12а = 1,125 коэффициент точности соответственно равен К. =2,5 и К. = 3,67). 
Реальный максимум нагрузочной характеристики ФК внутри интервала изменения возмущающего воз- 
действия возможен путем видоизменения ГФ2 согласно схеме (см. рис. 1). Данное видоизменение заключается 
в уменьшении числа фрикционных пар. Для этого из схемы необходимо исключить элементы 10, 11 и опереть 


тело 5 непосредственно на тело 2 через катки 12. 
В соответствии с этим найдем: 


Е =(Е-Ер›)У, 


где Е,> — сила трения ГФ2; Ё,› — управляющее воздействие ГФ2 (распорная сила на телах качения 8). 


Но 
Ер2 = 2,2180, 
Поэтому 
ий 
Е; = а = (2) 
1--Аоо 
Сила трения между парами ГФ]1 равна: 
Е 
Е = 2 [Анк е и 
Учитывая соотношение (2), найдем 
ЯОЙ. / (3) 


а 
(1+ зо.) 

Соотношение (3) не содержит разности в числителе, поэтому данная функция НИ при каких условиях не 
обращается В нуль. Она лишь асимптотически приближается К нулю при следующих условиях: 
= теоретически неограниченное возрастание коэффициента трения; 


— максимум в точке, соответствующей значению /. =1/19, [2—5]. 


Суммируя силы трения по соотношениям (2) и (3), получим 


3+ Даа 
не (4) 

(Ива) 
Продифференцировав функцию (4) по аргументу Г и приравняв нулю производную, найдем значение 

коэффициента трения, соответствующее максимуму функции: 
'_ 3 
Ле (5) 
120, 


Подставив в выражение (4) последовательно значения Риш, Ллах И приравняв друг другу полученные 


| 2 
в (т-+Т)+ т +12т | (6) 


Здесь т — относительная ширина интервала изменения коэффициента трения: т = Гаах / Лат: 


соотношения, найдем: 





При найденной величине 12а функция (6), имея максимум в точке, соответствующей значению (5), 


принимает одинаковые значения на границах интервала изменения коэффициента трения. 
Коэффициент точности в этом случае определяется отношением максимума функции (4) к ее любому 
граничному значению (при значениях Лии или Гиах ). На основании этого получим: 


___ 9+ ах)? 
8 ах 93+ Лах 0) 


Установим значение параметра 120, ‚› При котором величина Ку минимальна. При изменении величины 





(7) 


т 


{20 изменяется положение точки максимума Г. Только соотношение (6) устанавливает равенство сил трения 
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на границах интервала изменения /, поэтому изменение положения точки максимума функции (4) приводит к 


нарушению упомянутого равенства. В таком случае коэффициент точности должен определяться отношением 
максимума функции (4) к наименьшему граничному значению. Для этого необходимо установить тенденции 
изменения граничных значений функции (4), связанные с варьированием величины 100 [6-10]. 


Для решения данной задачи представим новую величину тангенса угла давления в виде произведения 
коэффициента варьирования К’ и базового значения тангенса угла в соответствии с выражением (6), т. е.: 
120; = Юра, . (8) 

На основании этого можно записать, с учетом (4): 

(Зт+ КУ пах 9.) Лпах 5 ЗК ах) Ушах 
тк хо’ Ч+К/ик од? 


Здесь левая часть соответствует силе трения ФК при значении Х = И › Правая часть — силе трения 





при значении / = Д‚ах. Решение последнего неравенства относительно коэффициента варьирования имеет вид: 


—_ 2 
к 5 (т+1)—\/(@т+1)7 +12т а 
2 пах 80! 





Так как 120>0, дробь в скобках левой части решения отрицательная, поэтому при К >1 составленное 
неравенство, безусловно, выполняется, т. е. при смещении точки максимума /’ в область меньших значений 


для вычисления коэффициента точности необходимо использовать величину силы трения ФК, соответствую- 
щую значению Диах, и, наоборот, при 1> К >0 необходимо принимать величину силы трения, соответствую- 


щую значению Ди [11-13]. 
Используя данный вывод, составим неравенство вида К.| >К., где К.| — коэффициент точности, 
вычисленный с учетом равенства (8). Имеем: 
2 2 
0+ Как 180) (ТА тах 189) 
К(3+К/ нах а!) 3 + ах а 


Неравенство составлено с учетом того, что К >1. 
Преобразование полученного неравенства к виду 


3+(1+К) аква + Каха > 0 


показывает справедливость предположения о том, что Ки > Ку . 





Исследуем соотношение величин коэффициентов точности при смещении точки максимума функции 
(4) в область больших значений, т. е. при К <1. Тогда К. > Ку, или 


и. 2. 
(т+ К/ пах 29.) 5 (тах 20%) 
К(Зт-+К/ пах 0) 3+ Лпах (80 
Преобразование составленного неравенства к виду 


Зт+(т- К) нах = а & >0 





и решение последнего дает 
т(ЗА ах 20.) 
(пах 15 — ИЛах 12а | 
Дробь В правой части полученного решения равна единице при значении 


РЕ (тт? +12т 
2 тах 


что полностью соответствует решению (6). Следовательно, неравенство К. > К. выполняется в случае, если 
К <. 








2 


Результаты исследования. Примем окончательное суждение о стабильности силы трения ФК при раз- 
личных формах нагрузочной характеристики. Для этого определим коэффициент точности ФК, когда функция 
(4) имеет максимум при значении Х = Г ‚ах . В этом случае функция монотонно возрастает в интервале измене- 


ния коэффициента трения. При этом справедливо равенство Ги„х =3/ 420. Получим: 
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УЕ.) _ (т+3)? 
УЕ) 8(т + 19) 





Кл 


Учитывая Куз > К. , найдем: 


Зт(3— тихо) + т(3—т) Г? ака +9(1— Г лахе00 <0. 
Для всех фрикционных материалов, применяемых в качестве пар трения ФК, т>3. Согласно (6), 


лак >1 › Поэтому очевидно, что разности в скобках полученного неравенства отрицательные и сделанное 


выше предположение справедливо. 

Таким образом, проведенный анализ показывает, что наибольшая стабильность выходного параметра 
ФК будет в том случае, когда функция (4) имеет максимум внутри интервала изменения коэффициента трения и 
принимает одинаковые значения на его границах. 

При т =8и Лиах = 0,8 получим 19 =14. Тогда РГ =0,214. В этом случае максимум функции (3) будет 
при значении /.=0,071, что практически идентично нижней границе интервала / = /ш= 0,1. Функция (3) 
убывает внутри интервала изменения коэффициента трения. При таких исходных параметрах К. =1,04, а 
Ктз=1,68. 

Обсуждение и заключения. Как видим, модернизированная принципиальная схема ФК позволяет тео- 


ретически получить очень высокую стабильность силы трения. Однако вследствие относительно большого зна- 
чения параметра 120 усилие РЕ. используется неэффективно. ФК не должен обращаться в ноль в интервале 


изменения коэффициента трения выходного параметра основной фрикционной группы (ОФГ) и при наличии 
максимума функции нагрузочной способности ФК. Необходимым условием этого является передача чувстви- 
тельными элементами дополнительной фрикционной группы ДФГ ее полной нагрузки. Дополнительное усло- 
вие существования максимума можно сформулировать следующим образом: при равном числе пар трения обе- 
их фрикционных групп чувствительные элементы ОФГ передают часть ее полной нагрузки; при меньшем, чем 
в ДФГ, числе пар трения чувствительные элементы передают полную нагрузку ОФГ. 
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Анализ эффективности системы внешнего обогрева мерной емкости насосной установ- 
ки от выхлопной системы двигателя внутреннего сгорания" 


С. О. Киреев', В. Н. Степанов?, М. В. Корчагина}, А. В. Ефимов" 
р р 


1,2,3,4 
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Введение. Статья посвящена исследованию эффективно- 
сти обогрева или поддержания рабочей температуры жид- 
кости, находящейся в мерной емкости передвижной 
насосной установки, от выхлопной системы двигателя 
внутреннего сгорания. Цель работы — совершенствование 
процесса проектирования и снижение затрат на натурные 
эксперименты. 

Материалы и методы. Для решения поставленной задачи 
были использованы методы численного анализа при рас- 
чете систем обогрева мерной емкости с внешним распо- 
ложением трубопровода от выхлопной системы двигателя 
внутреннего сгорания (ДВС) на примере широко распро- 
страненных установок ЦА-320, УНБ, АЧФ и др. Применен 
нестационарный нелинейный решатель газодинамических 
процессов (З1етеп$ ЗТАВ-ССМ+), позволяющий оценить 

правильность постановки задачи и значительно сокраща- 
ющий затраты на полномасштабные натурные испытания. 

Результаты исследования. Исследование проводилось для 
условий работы цементировочного агрегата на шасси 
«Камаз-43118» с насосом «СИН-32» и приводом от короб- 
ки отбора мощности на коробке переключения передач 
двигателя шасси. 

В расчетной модели применен конвекционный тип тепло- 
передачи между корпусом мерной емкости и окружающим 
воздухом; между корпусом мерной емкости и жидкостью; 
между трубопроводом выхлопной системы шасси и воз- 
духом окружающей среды; между трубопроводом вы- 
хлопной системы и выхлопными газами. 

В результате исследования получены графические зави- 
симости изменения температуры жидкости мерной емко- 
сти в контрольных точках; распределение поля темпера- 
тур жидкости в мерном баке; распределение поля темпе- 
ратур окружающего воздуха; линии потоков и поля скоро- 
стей окружающего воздуха и выхлопных газов ДВС. 


Гигодиспоп. ТЬе Веайп> еРйслепсу ог гесшаНоп оф Фе орег- 
айпе {етрегаеге о# Ве Йи14 1осайе@ ш Фе гацое 1апК оРа 
шобе ритрше чпй (РО) и5ше Ве ехБаи5Ё зузет оЁ фе 
пиегпа! сотбазНоп епеше (СЕ) 1$ зад1ед. Тве рарег обес- 
Пуе 15 Ю пиргоуе Ше 4ез1еп ргосезз ап4 10 гефиасе соз8 оЁ 
Не!4 ехрегитеп. 

Маепта[5 ап Мефо4з. То зо]уе 1$ рго ет, питейса] 
апа[уз1$ пео4$ ууеге изе4 \уВеп сасШайпе Фе Беание зуз- 
{еп15$ ОР Ше гапое {апК \ИЬ Фе ежегпа| 1осаНоп оЁ фе р!ре- 
Попе Нот Ше ТСЕ ехВаи${ зу$ет Бу фе ехатр!е ое улаеу 
изе4 шаПайотз$ (ТзА-320, ОМВ, АСКВЕ, ек.). З1етепз 
ЗТАВ-ССМ+, а поп-баНопагу поп-Шпеаг зо]уег оЁ газ- 
упапус ргосеззез, \’аз изе4, умсН еуаайе$ Ше соггесвтез$ 
о# Фе рго Мет зайетепЕ ап4 гефисез 1ет1Исап Ну Фе со 
о# № -5са]е Не! 1е545. 

Везеагсй Кезий5. Тре заду уаз сопдисе4 Юг Фе орегайп 
сопоп$ оЁ фе сетепипе ипй оп фе Ката7-43118 сВаз$15 
УИ Фе $П\-32 ришр ап4 а апуе Вот Фе ро\уег 1аКе-оР 
аНасЬтлепе оп Фе реагфох оЁ Ве сВаз$15 епрше. ш Ще са[си- 
1аноп плоде], а сопуесйоп-6уре Беай шапзег \’аз аррШеа Ъе- 
Смееп Ше Бо4ду оЁ 1е сайсое 1апК апд Ше зитоипате аи; 
Бебуееп Ше Бо4у о# Ше гаисе {апК апа Фе Паи1а; Бебмееп 
Фе сВаз$15 ехВаизЕ зу$бет 4асё ап4 ат еп аш; Бебмееп Ве 
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Обсуждение и заключения. В ходе исследования выясни- 
лась недостаточная эффективность рассматриваемой кон- 
струкции. Анализ данных позволяет предложить усовер- 
шенствованную конструкцию обогрева мерной емкости. 
Результаты, полученные в работе, могут быть использова- 
ны при расчетах подобных устройств, применяемых в 
установках гидроразрыва пласта и цементирования сква- 
жин. 
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Введение. Залежи трудноизвлекаемых углеводородов в северных широтах России приурочены к низко- 
проницаемым, слабодренируемым, неоднородным и расчлененным коллекторам. В связи с необходимостью 
освоения таких нефтегазоносных территорий возникает потребность в гарантированной плюсовой температуре 
жидкости, находящейся в мерной емкости установки. Эта проблема сохраняет актуальность при строительстве 
и капитальном ремонте скважин. 

На качество жидкости гидроразрыва пласта (ГРП), бурового и тампонажного растворов влияют следую- 
щие факторы: 

— типы пород, взаимодействующие с названными системами; 
— глубина скважины (давление); 
— температура на забое. 

Примем два утверждения стандарта ТНК-ВР «Соблюдение контроля качества при проведении ГРП»: 

1) жидкие химические реагенты, предназначенные для смешивания на кустовой площадке, должны храниться 
или перевозиться при температуре более 15 °С; 

2) хранение, перевозка, испытания и закачка всех жидких реагентов должна производиться при температуре 
>15 °С, независимо от погодных условий [1]. 

Температура системы не должна быть ниже 8-10 °С при всех видах работ по приготовлению жидкости 
затворения, продавочных систем и систем бурового и промывочного растворов [2]. 

В литературе по нефтегазовому оборудованию 1979 года упоминается использование для обогрева си- 
стемы выхлопных газов двигателя шасси или палубного двигателя (цементировочный агрегат ЦА-320А) [3]. 
Ссылки на продувочные системы осушения гидроблока насоса высокого давления есть в двух руководствах по 
эксплуатации: «Установки насосные УНБ-1000, УНБ-800, УНБ-630» и «Цементировочный гидроразрывный 
агрегат АЧФ-1050. ОРЕТКОМ». В руководствах по эксплуатации «Установка двухнасосная передвижная УНП- 
320%40» и «Установка насосная УНБ-125-50 СО» отдельно выделены пункты по работе с системой осушки и 
обогрева. 

В представленной работе проведен анализ эффективности обогрева мерной емкости с внешним располо- 
жением трубопровода от выхлопной системы двигателя внутреннего сгорания (ДВС). 

Материалы и методы. Использованы методы численного моделирования. В качестве примеров рас- 
смотрены широко распространенные установки ЦА-320, УНБ, АЧФ и др. 

Применение численных методов исследования позволяет оценить правильность постановки задачи (с по- 
грешностью выбора идеализированной модели), значительно сокращает затраты на натурные полномасштаб- 
ные испытания или полностью доказывает конструкторско-технологическую несостоятельность задачи. При 
отсутствии статистических данных по исследуемой проблеме полученные результаты численных методов 
сужают спектр решаемых натурных задач. В качестве инструмента численного метода был применен нестацио- 
нарный нелинейный решатель газодинамических процессов (51етепз ЗТАВ-ССМ-.) [4]. 

Для численного исследования была выбрана задача нагрева мерной двухсекционной емкости. С этой це- 
лью использован цементировочный агрегат на шасси «Камаз-43118» с насосом «СИН-32» и приводом от ко- 
робки отбора мощности на коробке переключения передач двигателя шасси. В задней зоне тележки шасси 
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установлена мерная емкость объемом 2х2 м”. Она представляет собой цельносварную конструкцию из листово- 
го металла, имеющую два отсека с наклонным днищем и донными клапанами. В предлагаемой для расчета и 
анализа схеме выхлопная труба (труба обогрева) проходит под днищем мерной емкости и по стенке левой ем- 
кости поднимается вверх (рис. 1). 
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Рис. 1. Расчетная идеализированная схема цементировочного агрегата на шасси «Камаз-43118» с насосом «СИН-32» 


Результаты исследования. Исследование проводилось для следующих условий работы агрегата на ме- 

сторождении: 

— насосный агрегат неподвижен; 

— температура окружающей среды (воздуха) плюс 5 °С; 

— скорость ветра 1,0 м/с; 

— плотность жидкости в обеих секциях мерной емкости 1100 кг/м°; 
— температура жидкости в мерной емкости плюс 5 °С. 

В расчетной модели применен конвекционный тип теплопередачи между следующими средами: 

— корпус мерной емкости и окружающий воздух; 

— корпус мерной емкости и жидкость; 

— трубопровод выхлопной системы шасси и воздух окружающей среды; 
— трубопровод выхлопной системы и выхлопные газы. 

Для упрощения задачи в качестве выхлопных газов был использован воздух, чтобы исключить примене- 
ние решателя многофазных процессов. 

Параметры воздуха: плотность — 1,184 кг/м; теплопроводность — 0,026 Вт/м-К; удельная теплоемкость 
— 1003,6 дж/кг-К. 

Параметры жидкости в мерной емкости: плотность — 1100 кг/м"; теплопроводность — 0,569 Вт/м-К; 
удельная теплоемкость — 4217 дж/кг-К. 

Параметры материала стального корпуса мерной емкости и трубопровода выхлопной системы: плотность 
— 7832 кг/м?; теплопроводность — 63,9 Вт/м-К; удельная теплоемкость — 434,0 дж/кг-К; толщина листа кор- 
пуса мерной емкости — 4 мм; толщина трубы выхлопной системы — 3 мм [5]. 

Температура выхлопных газов на выходе из коллектора ДВС — 450 °С [5, 6]; скорость выхлопных газов 
на выходе из коллектора ДВС — 35 м/с [5]. 

Расчет (см. рис. 1) был выполнен из условия работы ДВС 1800 секунд (30 минут). Контрольные точки 
температуры жидкости были получены для левой и правой частей мерной емкости у донных клапанов на рас- 
стоянии 200 мм от днища, на расстоянии 1000 мм от днища (Т1) и на расстоянии 1000 мм от днища и 200 мм от 
задней стенки (Т2). Графики изменения температуры жидкости за время работы ДВС в течение 30 минут пред- 
ставлены на рис. 2. Поля температур и линии тока представлены на рис. 3—6. 
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Рис. 2. График изменения температуры жидкости мерной емкости в контрольных точках: 1 — температура в верхней 
точке левого отсека мерника; 2 — температура в левом отсеке (точка 1); 3 — температура в левом отсеке (точка 2); 
4 — температура в левом отсеке мерника у донного клапана; 5 — температура в правом отсеке (точка 1); 
6 — температура в правом отсеке (точка 2); 7 — температура в правом отсеке мерника у донного клапана 
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Рис. 3. Распределение поля температур жидкости в мерном баке (вид слева) 
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Рис. 4. Распределение поля температур жидкости в мерном баке (вид справа) 
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Рис. 6. Линия тока и поле скоростей окружающего воздуха и выхлопных газов ДВС 


Обсуждение и заключение. Результаты расчетов наглядно показывают, что за 30 минут работы ДВС и 
конвекционного теплообмена между горячей выхлопной трубой, воздухом и жидкостью в мерной емкости в 
контрольных точках у донных клапанов жидкость прогрелась на 0,1 °С в левом отсеке и на 0,11 °С в правом 
отсеке. 

Это позволяет сделать вывод о том, что при рассматриваемом варианте конструктивного решения боль- 
шая часть тепловой энергии направлена на обогрев кабины оператора. Это положительный фактор, но не реше- 
ние поставленной задачи. Проведенные исследования и построенные графические зависимости показывают, 
что для более эффективного распределения тепловых потоков внутри мерной емкости и в окружающей среде 
необходимо разместить выхлопную тубу непосредственно в мерном баке. Учитывая особенности конструкции, 
вход трубы целесообразно расположить внизу правой части емкости и, пропустив через левую часть, направить 
вдоль бака вверх в предохранительном канале. 

С целью определения возможности обогрева мерной емкости необходимо провести расчеты для новой 
конструктивной схемы. 
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Введение. Статья посвящена некоторым теоретическим и 
технологическим аспектам обеспечения герметичности 
соединений на основе магнитострикции, а также поляри- 
зации герметизируемой среды в состоянии внешнего ин- 
дуцирования. Исследовано управление шероховатостью 
поверхностей стыкуемых деталей для повышения плотно- 
сти стыка при их индуцировании внешним магнитным 
полем. Рассмотрено создание электромагнитных барьеров 
для перемещения молекул герметизируемой среды через 
герметизатор. Цель работы — обоснование технологиче- 
ских условий обеспечения герметичности подвижных 
соединений в описанных выше случаях. 

Материалы и методы. Условия обеспечения плотности 
стыка показаны как результат решения контактной задачи 
и как фактор, определяемый положениями молекулярно- 
механической теории трения. Приняты геометрические, 
эксплуатационные и трибологические условия герметич- 
ности соединений. Демпфирующие свойства неподвижно- 
го фрикционного контакта определены молекулярной 
составляющей. Представлено теоретическое и расчетное 
обоснование факторов, влияющих на плотность стыка. В 
качестве целевых результатов технологической подготов- 
ки поверхностей деталей соединения указаны уменьшение 
глубины сглаживания, уменьшение отношения шагов по- 
перечной и продольной шероховатости и увеличение 
площади контакта. Потеря герметичности определена как 
специфический трансфер молекул. Они переносятся в 
зону стыкуемых поверхностей или свободно перемещают- 
ся через герметизатор на стадиях сорбционного поглоще- 
ния, диффузии и десорбции. Преобладание какой-либо 
стадии происходит при изменении энтропии и обусловле- 
но температурой и давлением. Визуализированы схемы 
герметизации соединений в управляемом магнитном поле, 
зависимости магнитострикции и магнитострикционных 
напряжений от напряженности магнитного поля. 
Результаты исследования. Экспериментально исследова- 
на устойчивость герметизаторов в легколетучих и газовых 
средах при их поляризации и намагничивании во внешнем 
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поле. В первом случае сначала вектор магнитной индук- 
ции ориентировался перпендикулярно продольной оси 
соединения. Установлено падение величины магнитного 
потока при работе соединения под нагрузкой в течение 
268 часов. Общая наработка соединения составила 
1070 часов. Если же вектор магнитной индукции ориенти- 
ровался продольно оси вала, наработка до момента кор- 
рекции величины напряженности поля составила 87 часов. 
В газовой среде продолжительность работы соединения до 
момента корректировки напряженности составила 187 
часов при общей наработке 935 часов. 

Обсуждение и заключения. Проникающая способность 
герметизируемых сред уменьшается в ряду «газ — пар — 
жидкость». Она зависит от температуры на контакте со- 
единения. Разгерметизацию можно отследить по измене- 
ниям магнитного потока, определяемого собственной 
магнитной проницаемостью молекул герметизируемой 
среды при их проникновении на поверхность стыка. Что- 
бы повысить герметичность, необходимо подавить актив- 
ность молекул. Для этого применяются ионизация и инду- 
цирование в постоянном и переменном магнитном поле 
напряженностью <60 кА/м. 


Ключевые слова: герметичность, соединение, детали 
машин, индуцирование, плотность, магнитное поле, кри- 
сталлографическое упорядочение, контактная жесткость, 
анизотропные эффекты. 
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Введение. Теоретическая основа герметологии изложена в трудах Г. А. Голубева, А. В. Чичинадзе, 
В. П. Тихомирова, Л. А. Кондакова, Г. В. Макарова, Х. Х. Вахтера, Э. Майера, Р. Х. Уоринга и др. Известно, что 
герметичность обеспечивается главным образом двумя факторами: 

— стабильностью технологической плотности стыкуемых шероховатых поверхностей деталей герметизатора; 
— стабильностью свойств материала герметизатора и герметизируемой среды, исключающей проникновение и 
перемещение через уплотнитель ее молекул за пределы стыка. 

Условия обеспечения плотности стыка обосновываются с помощью решения контактной задачи и в 
существенной степени определяются положениями молекулярно-механической теории трения [1, 2, 3]. Отме- 
ченный подход остается одним из основополагающих в инженерии герметичности. 

Проникновение герметизируемых сред (жидкостей, газов) через герметизатор описывается физической 
моделью С. Слихтера, И. Козени, Т. Кармана и Д. К. Коллерова (для пористых полимеров). Из ее анализа сле- 
дует, что объемный поток среды через уплотнение определяется длиной свободно проходимого пути моле- 
кул [4, 5, 6] и зависит от: 

— высоты микронеровностей А, твердой поверхности, 
— плотности стыка Р, 

— относительного сближения &, 

— геометрических размеров уплотнения. 


Характерно, что процесс фреттинга инициируют граничные структуры, в качестве которых выступают 
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различные полимерные уплотнители, жидкости и пленки на поверхностях стыка при амплитудах контактных 
вибраций менее 0,05 мм. Указанные среды могут выступать катализаторами, что приводит к снижению техно- 


логически обеспеченной проектной величины удельного давления на уплотнитель и увеличивает утечки. 
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Амплитуду и частоту относительных смещений деталей соединения можно снизить за счет управления 
плотностью стыка, например, при увеличении номинального давления на уплотнитель. Исследования, прове- 
денные на полимерах в узлах трения [7], позволяют отметить, что плотность стыка и контактная жесткость в 
значительной степени зависят от усилия сжатия поверхностей и от их межмолекулярной активации: 

Вир =@5+6Р,, 


где Р„„ — сила трения; а — средняя интенсивность молекулярной составляющей силы трения; 5 — фактиче- 
ская площадь контакта; 6 — коэффициент механической составляющей силы трения; Р, — сила сжатия поверх- 
ностей. 

Из этого следует, что молекулярная составляющая а5 пропорциональна: 

— площади фактического контакта, 
— интенсивности взаимной адгезии контактирующих материалов. 

Оба названных показателя можно увеличить за счет пластического насыщения по плоскости сопряже- 
ния, в том числе при приработке деталей соединения. Аналогичные условия могут создаваться на основе маг- 
нитострикционных эффектов, в большей степени выраженных для текстурированных материалов. Такие усло- 
вия применимы и в случаях, когда повторные нагружения деталей соединения затруднены или не предусмотре- 
ны нормами эксплуатации. 

Уменьшение проникающей способности герметизируемых сред прогнозируется при формировании по- 
движной сферолитной структуры герметизатора или его модифицирования порозаполнением низкомодульны- 
ми материалами [8]. Однако требуют особого рассмотрения ситуации, складывающиеся при существенно отри- 
цательных температурах и растрескивании полимерного герметизатора. В этих случаях отмечаются различия в 
температурных коэффициентах линейного расширения с металлами, поэтому представляет интерес возмож- 
ность управлять проникающей способностью молекул герметизируемой среды путем создания для них элек- 
тромагнитных барьеров. 

Итак, цель настоящей работы — изучение технологических условий обеспечения герметичности со- 
единений на основе индукционных эффектов, создаваемых в материалах деталей соединения и в герметизируе- 
мой среде. 

Материалы и методы. Условия герметичности соединений следуют из выводов задачи Ляме, в кото- 
рых граничными факторами выступают: 

— геометрические (например, длина / контакта поверхностей и диаметр 4 деталей сопряжения); 

— эксплуатационные (продольная осевая сила Ко. или крутящий момент М,» стремящиеся сдвинуть одну де- 
таль относительно другой); 

— трибологические (коэффициент трения /), определяющие момент распрессовки или проворачивания. 

Фрикционное стопорение обеспечивается при минимизации амплитуды относительных смещений по- 
верхностей, при контактных напряжениях, не превышающих область демпфирующей способности материа- 
лов [7]. 

Демпфирующие свойства неподвижного фрикционного контакта определены, в частности, молекуляр- 
ной составляющей. Причем сопряжения с текстурированной (анизотропной) структурой в направлении обеспе- 
чения высокой демпфирующей способности допускают более высокое значение критических вибрационных 
нагрузок в широком диапазоне частот [8]. 

Кроме того, при внешнем индуцировании изменения параметра шероховатости А„ вследствие магнито- 
стрикционного эффекта способны повлиять на плотность стыка. Это объясняется сопоставимостью величин 
магнитострикции и размеров блоков структурных составляющих (около 10 “10° м). Причем технологическое 
текстурирование по кристаллографическим направлениям с выраженными максимумами упругих или пласти- 
ческих свойств материалов позволяет при однонаправленном положении вектора магнитной индукции добиться 
максимального магнитострикционного эффекта. Это приводит к тому, что в разборных соединениях с натягом 
величина номинального давления возрастает, т. к. часть напряжений, расходуемых на податливость деталей, 
компенсируется магнитострикционными напряжениями в пределах упругих свойств. В неразборных, наоборот, 


часть напряжений может расходоваться на пластическое насыщение контакта. По этой же причине фактический 
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натяг оказывается больше расчетного. 

Коэффициент плотности стыка Р при магнитострикционном эффекте: 

Р=(Итм) + тв) Жем) / Гобщ- (0 

Здесь Грин) — объем материала шероховатого слоя с текстурой 1/1; Г›мии, — объем материала несжи- 
маемого уплотнителя, выдавленного в зазоры; А — коэффициент магнитострикции; Г»„ — общий объем шеро- 
ховатого слоя. 

Известно, что наличие волн приводит к уменьшению опорной площади контакта в 5—10 раз [9]. Следо- 
вательно, решение контактной задачи в аспекте обеспечения герметичности соединения требует учета волни- 
стости поверхностей. 

Высотные и шаговые параметры шероховатости взаимосвязаны, поэтому технологическая подготовка 
поверхностей деталей соединения должна быть нацелена на следующие результаты: 

— уменьшение глубины сглаживания К,; 
— уменьшение отношения шагов поперечной и продольной шероховатости; 
— увеличение площади контакта и отношения К, / Ко. 

Прямой зависимости между значением допуска и высотой неровностей, как правило, нет, однако следу- 
ет заметить, что впадины на профилограмме поверхности оказываются резервуарами, в которых скапливаются 
побочные продукты, инициирующие поверхностное разупрочнение. 

При обеспечении плотности стыка за счет натяга податливость 5» соединения представляется в виде: 

5м = + (Ему (ВК М, (2) 
где х (1) — нормальное контактное смещение в функции времени; & (В) — тензор магнитострикции при ре- 
гламентированном кристаллографическом упорядочении, учитывающий пластические деформации в покрытии 
и изменяющийся с магнитной индукцией; А` — коэффициент кристаллографического упорядочения; №* — ди- 
намическая нагрузка. 

При работе соединений в условиях гидростатического давления герметизируемой среды (жидкости или 
газа) потеря герметичности выражается в определенной форме трансфера молекул. Они переносятся в зону 
стыкуемых поверхностей или свободно перемещаются через герметизатор на стадиях сорбционного поглоще- 
ния, диффузии (с выравниванием концентрации среды в герметизаторе) и десорбции. Преобладание какой-либо 
стадии происходит при изменении энтропии и обусловлено температурой и давлением. 

Простейшим случаем фазового переноса в уплотнениях является течение вязкой несжимаемой жидко- 
сти в пористой среде. При этом определение величины утечки следует из закона Дарси. Однако в данном случае 
не учитывается, что уплотнительный стык, в отличие от объемного пористого тела, образован в результате кон- 
такта двух поверхностей. Поэтому с ростом нагрузки, а также с приработкой поверхностей отдельные пятна 
контакта могут образовывать изолированные объемы (эффективные каналы). Их число с ростом нагрузки 
уменьшается до момента образования гидродинамической пленки по периметру уплотнения [10-11]. 

Таким образом, принимая во внимание проводящую способность большинства герметизируемых сред, 
целесообразно исследовать обеспечение герметичности с позиций поляризации и внешнего индуцирования. 

Согласно [12], при прохождении среды через источник электрического тока или ионизатор молекулы 
получают заряд величины 4. Если заряженная частица, имеющая начальную скорость у, движется в однородном 
магнитном поле вдоль линий магнитной индукции В внешнего поля, то угол о между векторами у и В равен ну- 
лю. Магнитное поле на частицу не действует (сила Лоренца также равна нулю), и она движется равномерно и 
прямолинейно. 

Если же заряженная частица движется со скоростью у в магнитном поле перпендикулярно вектору В, то 


сила Лоренца Р =4[УВ] оказывается постоянной по модулю и перпендикулярной к траектории частицы, что 


исключает ее перемещение (рис. 1, а). 
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6) 


Рис. 1. Схемы герметизации соединений в управляемом магнитном поле: в однородном поле торообразного индуктора (а); 
в неоднородном поле соленоида (6); 1 — вал; 2 — втулка; 3 — магнит; 4 — заряженная частица среды; 5 — источник 
электрического тока (или ионизатор); 6 — блок автоматического управления соленоидом с генератором напряжения); 

у — орбита положительно заряженной частицы среды 


Предполагая направление движения среды, зная ее электрический потенциал, вязкость, плотность, про- 
никающую способность и изменяя напряженность поля и положение вектора магнитной индукции, возможно 
обеспечить силу Лоренца в величинах, достаточных для полного подавления движения. Тем самым реализуется 
эффект, реактивный гидростатическому давлению. 

Аналогично в неоднородном магнитном поле составляющая магнитной индукции В` создает силу ЕЁ", 
выталкивающую частицу в область слабого поля (см. рис. 1, 6). Изначально частица движется по радиусу по- 
стоянной величины со скоростью у. Однако собственное магнитное поле противодействует внешне наведенно- 
му, и этим объясняется крайняя неустойчивость его состояния. В результате частица выталкивается в область 
поля с пониженной напряженностью. 

Эффект герметизации основан на создании магнитным полем реактивного отталкивания предваритель- 
но заряженных молекул среды, проникающих на поверхность стыка. Причем степень герметизации зависит от 
магнитной индукции, величина которой устанавливается через силу тока с учетом относительной магнитной 
проницаемости и формфактора индуктора. 

Можно предположить, что условия вакуумизации стыка создаются при сближении поверхностей дета- 
лей сопряжения до проявления молекулярного (адгезионного) взаимодействия. При этом магнитная проницае- 
мость материалов деталей соединения представляется отношением магнитных потоков в материале Ф и в ваку- 
уме Фо (в = Ф/Ф)). При допущении Ф = Фу за единицу может быть принята магнитная проницаемость вакуума 
ф=1. 

В среде с неизменной магнитной проницаемостью индукция магнитного поля пропорциональна его 
напряженности. Следовательно, магнитный поток от внешнего поля, изначально определенный для вакуума, 
зависит от магнитной проницаемости среды. Ее присутствие на стыке соединения вследствие собственной маг- 
нитной проницаемости либо увеличивает магнитный поток (в парамагнитных средах), либо уменьшает его (в 
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диамагнитных). 

Таким образом, по изменению магнитного потока и, соответственно, магнитной индукции можно оце- 
нить и технологически скорректировать эффективность работы герметизатора в состоянии индуцирования. 

Результаты исследования. Описанные ниже кристаллографические направления могут быть приняты 
в качестве регламентированных при технологическом текстурировании деталей. 

1. Для разъемных соединений без полимерного герметизатора: [110] (для структур с объемноцентриро- 
ванной кубической (ОЦК) решеткой), [111] (с гранецентрированной кубической (ГЦК) решеткой решеткой), 
[1000] (с гексагональной плотноупакованной (ГПУ) решеткой). 

2. Для неразъемных соединений и соединений с герметизатором: [100]-[111] (с ОЦК-решеткой), [100]- 
[119] с ГЦК-решеткой), [0001]-[2110] (с ГПУ-решеткой). 

При установлении вектора магнитной индукции в указанных направлениях магнитострикция /. в кри- 
сталлах оказывается полеориентированной. Это связано с тем, что ее величина нелинейна с напряженностью 
поля и определяется преимущественно положением оси легчайшего намагничивания [13, 14]. 

При отсутствии предварительного текстурирования среднее изменение шероховатости для структур с 
решетками кубической сингонии представляется в виде: 


— 29 ‚ 
АЕ <[ Зена] , (3) 
5 5 
где Алоо, Лли — константы магнитострикции вдоль кристаллографических направлений [100] и [111]; “` — сред- 
ний размер дисперсных структурных составляющих. 


В случае текстурирования в направлениях [100] или [111] при ориентации поля по нормали к поверхно- 
сти (рис. 2) изменения шероховатости (К„) составят: 


0< ЛВ, < А/К, (4) 
ОА аль ИК, (5) 


где А — коэффициент потерь на размагничивание. 
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Рис. 2. Зависимость величины магнитострикции Лот напряженности Н магнитного поля 


Испытаны затворы с седлом из стали Х12М с текстурированным в направлении [111] покрытием фер- 
ритом кобальта СоЕе›О4 толщиной 20х10 м (К тах 2,4 мкм, "1 = 52510? и 72 = 50х10 3 м) и клапаном из латуни 
на величину потока утечки О, мкПа-м\/с гелия (Т = 77 К). Установлено, что при намагничивании соединения в 
поле 60 кА/м эта величина составляет 780-10240. Число циклов срабатывания затвора — 200-5000 соответ- 
ственно. При аналогичных условиях испытаний пары «седло — клапан» в отсутствие магнитного поля и пред- 
варительного текстурирования седел величина О составила 960—15340 мкПа-м\/с соответственно. 

При увеличении напряженности магнитного поля не установлены существенные изменения работоспо- 
собности соединения. По-видимому, это связано с тем, что эффективные магнитострикционные напряжения 
образуются в полях <40 кА/м (рис. 3). Данный вывод согласуется с результатами исследования магнитострик- 
ЦИИ. 
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Рис. 3. Зависимость величины магнитострикционных напряжений с от напряженности Н магнитного поля 


Изменение свойств герметизируемой среды, материалов герметизатора и температуры сказывается на 
проникающей способности и герметичности в целом. Проникающая способность сред уменьшается в ряду «газ 
— пар — жидкость». С учетом этого исследовалась работоспособность герметизаторов в легколетучих и газо- 
вых средах при их поляризации и намагничивании во внешнем поле. 

На торце неподвижного разъемного соединения вала с втулкой из меди крепился источник электриче- 
ского тока, один из полюсов которого фиксировался к торцам втулки. При этом второй полюс имел минималь- 
ный зазор с образующей вала. 

Соединение помещалось между полюсами тороидального электромагнита, подключенного к блоку ав- 
томатического управления с генератором напряжения, вольтметром и веберметром. Вектор магнитной индукции 
ориентировался перпендикулярно продольной оси соединения. 

Соединение устанавливалось в камере, заполненной диамагнитным бензолом. В ней создавалось давле- 
ние 30-35 кПа. Втулку жестко фиксировали, после чего включали источник электрического тока и индуктор и 
устанавливали напряженность магнитного поля 60 кА/м. Вал механически нагружали в радиальном и осевом 
направлении относительно втулки по симметричному знакопеременному циклу при величине нагрузки 1,5 Нс 
частотой 50 Гц. Испытания проводили до момента снижения магнитного потока, регистрируемого по показани- 
ям веберметра. Этот момент характеризовал нарушение герметичности, выражающееся в проникновении бен- 
зола на поверхность стыка. Блок автоматического управления индуктором корректировал напряжение в обмотке 
тороида, в связи с чем устанавливалась напряженность магнитного поля, равная изначально заданной. 

Установлено падение величины магнитного потока при работе соединения под нагрузкой в течение 
268 часов. Общая наработка соединения составила 1070 часов. Критерием завершения испытаний было сокра- 
щение интервала наработки до момента коррекции величины напряженности магнитного поля. Она должна 
превышать 60 кА/м. Это обусловлено необходимостью увеличения продолжительности наработки до значения, 
соответствующего предшествующим циклам. 

При ориентировании вектора магнитной индукции продольно оси вала при аналогичных режимах ис- 
пытания наработка до момента коррекции величины напряженности поля составила 87 часов. 

Имитирование условий эксплуатации соединения в среде газа проводили аналогичным образом. На вал 
устанавливали ионизатор, который крепили к торцам втулки. Соединение помещали в соленоидальный индук- 
тор, плотность витков которого увеличивалась от торца втулки. Благодаря этому создавалось неоднородное маг- 
нитное поле, напряженность которого увеличивалась от торца соленоида. 

Соединение помещали в камеру под давлением 30-35 кПа с парамагнитным азотом, после чего включа- 
ли ионизатор и индуктор. Вал нагружался на режимах, описанных выше, с периодическим фиксированием маг- 
нитного потока. Момент изменения его начального значения характеризовал нарушение герметичности и про- 
никновение азота на поверхность стыка. При этом блок автоматического управления индуктором обеспечивал 
снижение напряжения в обмотке соленоида до величины, соответствующей напряженности магнитного поля 
60 кА/м. 

Продолжительность работы соединения до момента корректировки величины Н составила 187 часов 
при общей наработке 935 часов. Большая продолжительность испытаний вызывала необходимость коррекции 
напряженности магнитного поля при интервале наработки менее 187 часов. 
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Обсуждение и заключения. Уменьшение утечек через герметизирующие устройства обеспечивается за 
счет управления плотностью стыка (удельным давлением р.) при снижении амплитуды и частоты относитель- 
ных смещений деталей соединения. Последнее достигается за счет увеличения площади фактического контакта 
и интенсивности взаимной адгезии контактирующих материалов, в частности благодаря магнитострикционным 
эффектам. При этом уровень магнитострикционных напряжений не должен приводить к пластическим дефор- 
мациям уплотнителя (т.е. р. = 110„) и может быть обеспечен в полях до 40—50 кА/м. При изменении вслед- 
ствие магнитострикции комплексного параметра шероховатости А отмечается увеличение площади опорной 
поверхности контакта, преимущественно за счет того материала, у которого модуль упругости оказывается 
меньше, а магнитострикция выше. 

Влияние увеличивающейся при магнитострикции площади опорной поверхности контакта на герме- 
тичность возрастает с уменьшением глубины сглаживания К, и отношения шагов поперечной и продольной 
шероховатости. В результате за счет компенсаторного эффекта со стороны магнитострикционных напряжений 
растет номинальное давление в соединениях. Магнитострикцию стимулирует кристаллографическая магнитная 
анизотропия, в частности при технологическом текстурировании. 

В ряду «газ — пар — жидкость» уменьшается проникающая способность герметизируемых сред. Она 
оказывается зависимой от температуры на контакте соединения. Разгерметизация диагностируется по измене- 
ниям магнитного потока под влиянием собственной магнитной проницаемости молекул герметизируемой сре- 
ды при их проникновении на поверхность стыка. Можно повысить герметичность, если подавить активность 
молекул. Для этого применяются ионизация и индуцирование в постоянном и переменном магнитном поле 
напряженностью <60 кА/м. 
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Введение. Современная наука признает актуальной про- 
блему адекватного моделирования физических процессов 
и объектов. Автоматизация работ в этой сфере повышает 
точность моделирования и экономит средства на натурном 
моделировании объектов изучения. Целью данного иссле- 
дования было создание компьютерной программы, кото- 
рая по исходному набору параметров автоматически фор- 
мирует безразмерные критерии подобия. 

Материалы и методы. В основе расчетного алгоритма 
лежит л-теорема теории подобия и основные положения 
теории размерностей. 

Результаты исследования. Создана компьютерная про- 
грамма, которая позволяет получить все теоретически 
возможные критерии подобия для моделирования изучае- 
мого физического процесса или объекта. 

Обсуждение и заключение. Полученные результаты могут 
быть использованы в промышленном и научном модели- 
ровании физических объектов исследования, расчете но- 
вых критериев подобия, решении задач описания сложных 
процессов и т. д. 


Ключевые слова: критерий подобия, теория подобия, 
теория размерностей, л-теорема, безразмерный комплекс, 
теорема Кирпичева — Гухмана. 
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Введение. Критерии подобия — это зависимости между безразмерными комбинациями параметров, 
описывающих объект или физический процесс. На основе таких критериев можно создавать математические 
модели изучаемых объектов [1—4]. Критерии подобия позволяют доказать соответствие модели и объекта изу- 


чения. Они незаменимы, если не сформировано математическое описание процессов по данным экспериментов 
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или наблюдений. Фактически критерии подобия являются безразмерными степенными комплексами наборов 
физических параметров, описывающих объекты и процессы. По форме они идентичны и численно равны для 
подобных объектов или процессов. 

Материалы и методы. Алгоритм формирования критериев подобия основан на второй теореме теории 
подобия, так называемой л-теореме. В [5] она сформулирована так: «Всякое полное уравнение физического 
процесса, записанное в определенной системе единиц, может быть представлено функциональной зависимо- 
стью между критериями подобия, полученными из участвующих в процессе параметров». 

Из этого следует, что, анализируя и комбинируя физические параметры, описывающие объект изуче- 
ния, можно составить критерии подобия данного объекта или физического процесса [5, 6]. 

При моделировании объектов исследования на основе теории подобия наиболее часто используется 
анализ размерностей физических параметров, описывающих объект в его наиболее важных точках [6, 7]. 

Исследователи выбирают определяющие физические параметры, описывающие объект изучения (в 
частности, давление, вязкость, плотность, скорость, температуру, сопротивление, силу тока) и, работая с их 
размерностями, получают теоретически возможное количество критериев подобия [8]. Эти критерии далее 
можно положить в основу математического или натурного моделирования 

Если у исходного объекта изучения и модели критерии подобия в критически важных зонах численно 
равны, то они физически подобны [6], согласно третьей теореме подобия (теорема Кирпичева — Гухмана) [9]. 
Следовательно, для создания, например, модели движения потока по контуру труб необходимо найти внешний 
вид и численное значение критериев подобия, описывающих реальный объект, а затем построить модель или 
лабораторную установку так, чтобы в ней значения критериев подобия были численно равны значениям этих 


же критериев на реальном объекте исследования. 


Результаты исследования. Программа М№опОпаСгиЕогтег 1.0 позволяет исследовать физические объ- 
екты и процессы. При этом можно выбрать существенные физические величины (рис. 1), которые определяют 
базовые параметры объекта изучения, и автоматически рассчитать все теоретически возможные критерии по- 
добия и безразмерные комплексы (рис. 2), сформированные на их основе. 
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Рис. 1. Выбор исходных физических параметров объекта исследования 
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Рис. 2. Рассчитанные критерии подобия 


Описание алгоритма программы. Представленный в программе алгоритм формирования критериев 
подобия базируется на размерностях выбранных физических параметров. Он реализуется в пять шагов. 

Шаг 1. Определение физических параметров, которые в значительной степени влияют на физические 
свойства объекта изучения. От правильности выбора параметров, отсечения несущественных зависит точность 
результатов. Большое количество выбранных параметров увеличивает число получаемых безразмерных ком- 
плексов, но не повышает точность сформированной на их основе модели. 

Шаг 2. Составление матрицы размерностей для выбранных физических величин. В ней по строкам 
расположены размерности отобранных параметров в соответствующих степенях. Размерность матрицы соот- 
ветствует количеству параметров и их единиц измерения, т. е. в матрице строк столько же, сколько и выбран- 
ных параметров, а столбцов столько же, сколько единиц измерения. Например, строка, соответствующая физи- 

мкг 





ческому параметру «сопротивление электрическое, №», который имеет размерность четыре единицы в 
СА 
матрице размерностей будет иметь вид |2 1-3 —2|. 

Шаг 3. Определение количества независимых параметров. Количество независимых параметров К со- 
ответствует рангу матрицы размерностей, т. е. наибольшему порядку неравного нулю определителя, который 
составлен на базе строк исследуемой матрицы. 

Шаг 4. Формирование групп независимых параметров. На этом шаге формируются сочетания по А 
штук из т параметров, которые составляют безразмерный комплекс. То есть единицы измерения, которые вхо- 
дят в размерность выбранных параметров, взаимно сокращаются, в результате данное сочетание выбранных 
физических параметров не имеет размерности. Говоря в общем, комплекс образуют такие параметры, частичная 
матрица размерностей которых имеет как минимум один определитель А-го порядка, не равный нулю. 

В соответствии с комбинаторной формулой размещений теоретически возможное количество таких 
комплексов [1]: 

с _ т! , 
К! (т-—К)! 
где т — общее число параметров, А — ранг полной матрицы размерностей. 
Далеко не все полученные комплексы выбранных параметров составляют комплекс с матрицей А-го 


порядка, поэтому весьма часто количество групп после проверки бывает меньше в Е 


Шаг 5. Определение выражений для критериев подобия на основе вычисленных групп независимых па- 
раметров. Выражения формируются по следующему алгоритму. 

1. Берем очередную группу параметров. Например, набор № 1 — Р\...Рь. 

2. Рассчитываем критерии подобия путем деления параметров очередной группы на параметры теку- 
щей. В частности, для набора № 1 могут быть составлены следующие критерии: 
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Простое деление параметра на набор не приводит к получению безразмерного выражения. Для этого 
нужно вычислить значения показателей степеней Х; 1... к для каждого критериального выражения. 


3. Размерности параметров в числителях и знаменателях критериальных комплексов должны быть рав- 
ны. Принимая во внимание этот факт, формируем и решаем систему линейных уравнений. Для создания такой 
системы приравниваем степени при одинаковых единицах измерения в числителях и знаменателях критериев, 
представляющих собой произведения сочетаний размерностей, выбранных в начале расчета физических пара- 
метров. В итоге получаем систему из 4 уравнений с ^ = 4 неизвестными. 

Получаем наборы по т-— критериев для каждой группы параметров. Максимальное количество безраз- 
мерных комплексов, которое можно рассчитать на основе т параметров, составит: 


п; = СА. (т-Ю). 
Эта величина может быть довольно большой. Например, для т =7 ид=3 п, =140. Однако рассчитан- 


ные критерии часто бывают одинаковыми и взаимообратными, поэтому общее количество отобранных крите- 
риев оказывается меньше, чем п;. Доля критериев, не прошедших проверку, зависит от исходного набора фи- 


зических параметров, отобранных для изучения. 

Интерфейс и функциональность программы. Программа Моп)ОиаСгиРогтег 1.0 написана на языке 
\У15иа1Ваз1с.Ме{. Она обладает интуитивно понятным \\!т4до\/5-интерфейсом и проста в эксплуатации. 

Пользователь выбирает физические параметры, существенные, на его взгляд, для моделируемого объ- 
екта или процесса (см. рис. 1). В результате выполнения команд «Сформировать все критерии» и «Убрать по- 
вторы критериев» на экран выводятся все критерии подобия (безразмерные комплексы), которые можно соста- 
вить на основе выбранных физических параметров (см. рис. 2). 

Команда «Найти стандартные критерии» позволяет установить наличие в списке известных (имеющих 
имена) критериев подобия. Кроме того, в МопОииСгиРогтег 1.0 предусмотрена возможность выводить на пе- 
чать весь список полученных критериев или его часть, сохранять список в виде отдельного файла, редактиро- 
вать базу стандартных критериев. 

Примеры использования программы 

Расчетный пример 1. В условиях постоянного напряжения найдем критерии подобия, описывающие 
прохождение тока в резисторе, дросселе и конденсаторе, подключенных последовательно [8]. 

Шаг 1. Выбираем физические параметры, существенные для моделируемого объекта. К ним можно от- 


в. 4. 42 2 
нести: силу тока, /[, А; индуктивность, Г, —_—>; емкость, С. >; сопротивление, В; — >; напряжение 
с^.А кг. м С.А 
кг. м2 
и, 5 и время (, с. Всего параметров т = 6, а базовых единиц измерения 4 = 4. 
С.А 
Шаг 2. Матрица размерностей указанных параметров имеет вид: 
м ею с А 
оо гоОо 
С |-2 4 2 
2 1 -2 —>2 
2 1-1 
1000 1 
|2 1-3 -—2 








Шаг 3. Ранг этой матрицы К = 3. Из этого следует, что число независимых параметров в группе будет 
равно трем, число критериев для каждой независимой группы т — А = 3. 
Шаг 4. Общее количество групп независимых параметров может составить Гы = г. =20. Комбинации 


параметров составляются автоматически описываемой программой. Алгоритм проверяет ранг частичной мат- 
рицы для каждой составленной комбинации. 
Например, для сочетания П.О частичная матрица будет выглядеть так: 
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000 1 
ПО =|2 1 -2 —2. 
2 1—3 - 


Ранг данной матрицы равен 3, следовательно, эта комбинация образует группу независимых парамет- 
ров. Не все из возможных 20 сочетаний являются безразмерными критериями подобия. Например, не является 
таким критерием сочетание СГ, потому что ранг матрицы не равен (. В контексте рассматриваемого исходного 
набора физических параметров только 15 можно считать критериями подобия. 

Шаг 5. Составление критериев подобия проходит через решение систем линейных уравнений, которые 
формируются на основе критериальных комплексов. Например, для сочетания 1/0 














ты 1 С Ра 
м= »„ 1о= › = 
[2 . [№3 „(а [22 . [№3 „(7724 [732 . [№3 „(4 
будут найдены значения степеней размерностей х и получены следующие критерии: 
1 С К 
= —, 2 = —, з=— з 
вы" РН" РЕ 


Для каждой из 15 возможных комбинаций исходных физических параметров будут получены по 3 кри- 
терия. Среди этих 45 критериев зачастую бывают повторяющиеся и взаимообратные. Для рассматриваемого 
примера МопОпаСгиЕогтег 1.0 выдает 11 уникальных критериев: 
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С05 . [05 С05 0 [05 
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О о 
7 , 3 , 9 = , 
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Расчетный пример 2. Рассмотрим процесс движения жидкости в горизонтальной трубе. Для этого 
сформируем критерии подобия, которые описывают падение давления в этом процессе [10]. 

Наиболее важными параметрами, которые описывают падение давления в трубе при прохождении по 
ней жидкости, можно назвать следующие: 


кг 
где Др —Щ давление, о ; 0 — скорость потока жидкости, 
м.с 


Др = Г (в, р. ь Г), 
м 
с 


кг 
; р — плотность жидкости, —; и — коэффи- 
м 


у кг 
циент динамическои вязкости, ——, Е— длина трубы, М. 
м:с 


Всего параметров т = 5, а базовых единиц измерения 4 = 3. Полная матрица размерностей выбранных 
параметров имеет вид: 


М кг с 
ит 1 
Др |- 1 -2 
Е, 09| 
р |-3 1 
о 0 -1 








Ранг приведенной матрицы К = 3. Число независимых параметров в группе — три и количество крите- 
риев для каждой независимой группы т — К = 2. Всего может быть "в = С =10 групп. 


Согласно расчетам программы, в данном примере из десяти возможных сочетаний независимыми 
группами являются только девять. Из этого следует, что всего можно получить 18 безразмерных комплексов. 
После удаления повторяющихся и взаимообратных программа выдает 6 уникальных сочетаний. 
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Важно отметить, что комплексы л) и п; являются выражениями критериев Эйлера и Рейнольдса соот- 
ветственно [3]. 

Обсуждение и заключение. Предложена компьютерная программа Моп)ииаСгиЕогтег 1.0. С ее помо- 
щью можно сформировать безразмерные комплексы — критерии подобия, описывающие объект изучения. Для 
этого предварительно выбираются физические параметры, которые определяют указанный объект в характер- 
ных пространственных и временных точках. Автоматизация научных изысканий с помощью предложенного 
решения повышает точность моделирования и экономит средства на натурном моделировании объектов изуче- 


НИЯ. 
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Введение. Статья посвящена исследованию криптографи- 
ческой системы МТКОЕпсгурь, расчету алгоритмической 
сложности разработки криптосистемы МТКОЕюсгур! и ее 
модификации. Цели исследования: разработка эффектив- 
ного постквантового криптографического алгоритма 
МТКОЕпсгурь обладающего высокой криптостойкостью к 
атакам с квантового компьютера, а также разработка мо- 
дификации предложенного алгоритма, анализ и экспери- 
ментальное доказательство его преимуществ. 

Материалы и методы. Предложено описание системы 
шифрования МТКОЕпсгуре. Изучена модификация рас- 
сматриваемого алгоритма, представлена блок-схема реа- 
лизации основанного на нем программного средства. При- 
веден пример работы программного средства и дана его 
характеристика. Достоверность результатов обоснована с 
помощью И-критерия Манна — Уитни. При проведении 
эксперимента использована сторонняя программная реали- 
зация криптографической системы КЗА. В исходный код 
всех трех программ МТКОЕпсгур, В$А, модификации 
МТВУЕпсгурЕ был внедрен элемент класса Зюру’аВ. Дан- 
ный класс предоставляет набор методов и свойств, которые 
можно использовать для точного измерения времени, затра- 
ченного на выполнение. Таким образом, появилась возмож- 
ность фиксировать результаты затраченного времени на всех 
трех основных этапах: создание ключей, шифрование и 
расшифрование сообщения. 

Результаты исследования. Доказаны преимущества раз- 
работанных криптосистем по характеристикам производи- 
тельности. Выполнено экспериментальное сравнение реа- 
лизованного алгоритма МТВОЕпсгур{ и его модификации. 
При этом обозначены все преимущества последней. 
Обсуждение и заключения. Экспериментально доказано 
преимущество использования модификации алгоритма 
МТКОЕпсгуре. Новое приложение на 25 % быстрее выпол- 
няет общую работу по генерации ключей, шифрованию и 
расшифрованию. Помимо этого оптимизируется исполь- 
зование внутренней памяти за счет уменьшения веса ис- 


* : 
Работа выполнена в рамках инициативной НИР. 
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ходного файла программы и размера секретного ключа. Кеу. \!Веп абетрИп? №ю сгасК а сарБенехь сгурюзтарЫс 
При попытке взлома шифротекста проявляется крипто- гобизпез$ ап4 сотр ехИу оЁ изше адпапиит аеогИ из аге 
графическая стойкость и сложность использования кван- зВо\уп. 


товых алгоритмов. 
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Введение. Толчком для разработки новой криптографической системы послужила статья [1]. В ней по- 
казано, что квантовые компьютеры потенциально угрожают взломом всем широко используемым криптогра- 
фическим алгоритмам. 

Представленная в данной статье программная разработка обладает криптостойкостью по отношению к 
возможным квантовым атакам и превосходит аналоги (например, криптосистему КЗА) по характеристикам 
скорости работы алгоритма и по количеству затрачиваемых ресурсов [2]. Этим обусловлена актуальность рабо- 
Ты. 

Объектом исследования является криптосистема МТВОЕпсгури. 

Предмет исследования — алгоритмическая сложность разработки криптосистемы МТКОЕпсгур и ее 
модификации. 

Цели исследования: разработка эффективного постквантового криптографического алгоритма 
МТКОЕпсгурф обладающего высокой криптоскойстью по отношению к атакам с квантового компьютера, а так- 
же создание модификации предложенного алгоритма, анализ и экспериментальное доказательство его преиму- 
шеств. 

В соответствии с поставленной целью были определены следующие задачи. 

1. Исследовать алгоритм работы криптосистемы МТКОЕпсгури. 

2. Разработать алгоритм модификации МТВОЕпсгуре. 

3. Реализовать программные средства криптосистемы МТВОЕпсгур! и ее модификации. 

4. Проанализировать и сравнить две программы друг с другом и с их аналогом — криптосистемой ВЗА. 

Материалы и методы. Рассмотрим описание системы шифрования МТКОЕпсгур". Криптографическая 
система с открытым ключом МТВОЕпсгур! использует операции над кольцом Й [Л] / (Х\ - 1) многочленов сте- 
пени, не превосходящей М- 1 [3]: 

а =а +а, *Х! + а *Х? + -- + аул * ХМ 1, 
где ао, а1, а2 ... ам_1 — целые числа. 
Операции сложения и умножения производятся как обычно, за исключением того, что ХМ заменяется 


№+2 
,Х 


на 1, Х\+1 заменяется на Х1 заменяется на Х? ит. д. 


Криптосистема определяется рядом параметров, основные из которых: №, ри д. Для сохранения стой- 
кости алгоритма необходимо, чтобы параметры р и 4 были взаимно простыми. 
Чтобы обеспечить высокую стойкость алгоритма к различным атакам, рекомендуется использовать 
следующие параметры (рис. 1): 
Обозначение М 9 9: | 99 аг Гарантированная стойкость 
МТКУ167:3 167 128 
МТКУ251:3 251 | 128 


р 
3|61 |20 181 Умеренный уровень стойкости 
3 
МТВУ503:3 503 | 256 |3 216 | 72 |55 Высочайший уровень стойкости 
2 
2 
2 


50 |24 | 16 | Стандартный уровень стойкости 


МТКУ167:2 167 | 127 
МТКУ251:2 251 | 127 
МТВУ503:2 503 | 253 


45 35 | 18 | Умеренный уровень стойкости 


35 |35 |22 | Стандартный уровень стойкости 
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155 100 | 65 | Высочайший уровень стойкости 


Рис. 1. Рекомендованные параметры 
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Результаты исследования 
Генерация ключей. Боб хочет передать сообщение Алисе. Для этого ему необходимы открытый и за- 


крытый ключи. Поэтому он случайным образом выбирает два малых полинома / и © из кольца усеченных мно- 
гочленов А. Малость полиномов означает, что относительно произвольного полинома по модулю 4, в котором 
коэффициенты равномерно распределены, у малого полинома они будут много меньше 4 [4]. Для определения 
малости полиномов используются числа Хи 4®, которые Боб выбирает самостоятельно. 

Полином / будет иметь коэффициентов, равных единице, (4/- 1) коэффициентов, равных минус еди- 
нице, и остальные, равные нулю. 

Полином © будет иметь 4® коэффициентов, равных единице, столько же коэффициентов, равных минус 
единице, и остальные, равные нулю. 

Боб должен хранить выбранные полиномы в тайне, так как любой, кому они станут известны, сможет 
расшифровать сообщение. 

Следующим шагом Боб вычисляет обратные полиномы /, и {, по модулю р и 4 соответственно, такие, 
что: 

Дх }» = 1(тоа р) и} х {. = 1(тоа 4). 

Если случайно эти обратные полиномы не существуют, то Боб возвращается назад и заново выбирает 
полином { 

Секретный ключ — это пара (}, /»), а открытый ключ й вычисляется по формуле: 

В =рхДхд (тоа 4). 

Шифрование. Алиса хочет отправить сообщение Бобу с помощью открытого ключа й. Для этого Алисе 
нужно представить свое сообщение в виде полинома т с коэффициентами по модулю р, выбранными из диапа- 
зона (-р/2, р/2]. Затем Алисе необходимо выбрать другой малый полином 7, который называется «ослепляю- 
щий», и вычислить шифротекст: 

е = (гхА+т)(тоа 4). 

Расшифрование. Боб получает от Алисы зашифрованное сообщение е и хочет его расшифровать. Пер- 
вым делом, используя свой секретный ключ, Боб вычисляет: 

а = {хе (тоа 4). 

Так как Боб вычисляет значение а по модулю числа 4, он должен выбрать его коэффициенты из диапа- 
зона (-4/2, 4/2] и затем вычислить: 

р =а (тоар). 
Наконец, Боб, используя вторую часть секретного ключа, получает исходное сообщение от Алисы: 
с = хЬ (тоар). 

Модификация алгоритма МТК ОЕпсгург. Как видно из описания алгоритма, полином / должен отвечать 
следующим требованиям: 

— полином Гявляется обратимым по модулю р, 
— полином Гявляется обратимым по модулю 4, 
— полином Гявляется малым полиномом. 

В самом алгоритме обратимость по модулю р и 4 гарантировалась следующим образом. Если вдруг ге- 
нерировался полином /, не обратимый по одному из модулей, то он отбрасывался и генерировался следующий 
— итак до тех пор, пока необходимый полином не будет найден. 

Предлагаемая модификация заключается в том, чтобы заменить полином Гполиномом вида: 

=1+рР, (1) 
где Е’ — малый полином. 

Данный подход имеет следующие преимущества. 

1. Из самого выражения (1) видно, что полином / всегда обратим по модулю р. Этот факт ускоряет ге- 


нерацию ключей, так как не требуется дополнительно ВЫЧИСЛЯТЬ },. 
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2. Поскольку / ' = 1 тоа р, то при расшифровании не требуется дополнительное умножение на /` ', что 
ускоряет сам процесс расшифровки. Закрытый ключ в этом случае будет представлять собой не пару ({, /›), а 
(Г). 

Генерация ключей. Как и в оригинальном алгоритме, вначале Боб выбирает параметры шифрования М, р, 4 
и числа 4% 4. Затем он случайным образом выбирает два малых полинома Ри © из кольца усеченных многочленов 
К. 

Вычисляет модифицированный полином /по формуле (1). 

Следующим шагом Боб вычисляет обратный полином /, по модулю 4: 
хр = 1(тоа 4). 
Если случайно обратный полином не находится, Боб возвращается назад и заново выбирает полином /. 
Секретный ключ — это полином Х а открытый ключ й вычисляется следующим образом: 
В =рхДхд (тоа 4). 

Шифрование. Шифрование остается без изменений, все в точности как в оригинальном алгоритме 
МТКОЕпсгуре. 

Алиса хочет отправить сообщение Бобу с помощью открытого ключа й. Для этого Алисе нужно пред- 
ставить свое сообщение в виде полинома т с коэффициентами по модулю р, выбранными из диапазона (-р/2, 
р/2]. Затем Алисе необходимо выбрать другой малый полином г, который называется «ослепляющий», и вы- 
числить шифротекст: 

е = (гхА-+т)(тоа 4). 

Расшифрование. Боб получает от Алисы зашифрованное сообщение е и хочет его расшифровать. Пер- 
вым делом, используя свой секретный ключ, Боб вычисляет: 

а = {хе (тоа 4). 

Так как Боб вычисляет значение а по модулю числа 4, он должен выбрать его коэффициенты из диапа- 
зона (-4/2, 4/2] и затем вычислить: 

Ь =а (тоар). 

Всб, на этом вычисления закончены, мы получили исходное сообщение от Алисы: В =т [5]. 

Доказательство работы модифицированного алгоритма. Для доказательства работы алгоритма рас- 
смотрим сам процесс расшифрования. 

Зашифрованное сообщение от Алисы имеет вид: 

е = (гхйА+т)(тоа 4). 
Боб использует свой закрытый ключ — полином 
а = {хе (тоа 4) = (фх (гхй + т))(тоаа) = (фх (гхрдх9+ +т))(тоа 4). 
В результате: 
а = (фт хад+тх ХУ (тоа а). 


Следующим шагом Боб получает полином Ь путем уменьшения коэффициентов полинома а по моду- 


Ь = а(тоар) = тх {(тоар) = (т + т хрх Ю(тоар) = т (тоар). 

Таким образом, мы проверили и доказали, что полином Ь действительно является исходным сообщени- 
ем т. 

Реализация алгоритма. В качестве языка программирования использован объектно-ориентированный 
язык программирования С#, относящийся к семье языков с С-подобным синтаксисом. Среда разработки — 
Мисгозой У15па1 Зшато 2015 Ещегризе. Главное преимущество данного программного продукта — приложение 
с графическим интерфейсом, который позволяет пользователю быстро разобраться в устройстве и схеме работы 
данного программного продукта. 


На рис. 2 приведена обобщенная блок-схема работы программного средства. 
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файла данных 


Запись результата в файл 













Расшифрование |-го блока 
файла шифротекста 


Запись результата в файл 


Рис. 2. Блок-схема алгоритма реализации программного средства 


Пример работы программного средства. МТВОЕпсгур! использует три постоянных параметра: №, р, 4. 
Они вводятся пользователем на панели «Параметры системы». 

Далее располагается панель «Генерация ключей». На вход подаются параметры, введенные пользова- 
телем на предыдущем этапе. Затем случайным образом формируются два полинома Ги = — такие, чтобы: 
— количество коэффициентов, равных единице и нулю, было равно числу, которое заранее определено в про- 
грамме; 
— степень полиномов соответствовала введенному параметру №. 

Затем используется теорема Евклида о нахождении наибольшего общего делителя для многочленов и 
ее обратный ход для вычисления двух обратных к / полиномов по модулю р и 4 соответственно. Вычисляется 


открытый ключ Й. 
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Пользователь вводит двоичный вид сообщения т на панели «Шифрование» в окне «Сообщение». Сто- 
ит отметить, что исходное сообщение необходимо разбить на блоки по Мбит, каждый из которых будет обраба- 
тываться отдельно и преобразовывать каждый блок в полином с коэффициентами {-1, 0, 1} (в данной програм- 
ме вместо значения «—1» используется «2». 

Для операции шифрования программе необходимо выполнить предварительную подготовку — сгене- 
рировать один ослепляющий полином. Он формируется по такому же принципу, что и полиномы /, ©. 

Полученный открытый ключ и полином 7 дают возможность зашифровать сообщение т, используя со- 
ответствующую формулу. Затем результат проверяется на принадлежность кольцу усеченных многочленов сте- 
пени, не превосходящей №- 1, и выводится на экран. 

Следующий блок реализует механизм дешифрации. Для этого последовательно выполняются действия 
с учетом кольца ДХИ(ХМ — 1). 

Результат проведенных операций отображен на рис. 3. 


В*] МТВУЕпсгур! = 
































Параметры Генерация ключей Шифрование 
криптосистемы 
ЬЙ Демонстрация одного раунда работы 
№М= |11 4-4. 9-3 Битовое сообщение: |10011101001 
и - Сгенерировать Г= (х^0) +к^ 1х” 2х” 3+7 6) +(х^10) 
а= 7 Зашифровать 











+ = 4х^0 + (-3х^1) + (3х^2) + (-3х^5) +Зх^7 +3х^11 
9=х^0+(х^1) +(х^3) +к^4+(х^9) +к^10 


Дешифрация 
Е = 333-3-3230233 Шифротекст: 310-1-31-211-31 
Выход В = 1102-3312-202 
Расшифровать 


Открытый ключ: (В ); 


Закрытый ключ: ($). 10011101001 


Рис. 3. Внешний вид программы МТВОЕпсгур! 


На рис. 4 отображены результаты работы программного средства в виде исходного текста, зашифро- 
ванного и расшифрованного. 


Я т ех+54 — Блокно] ] зн тех — Бло Мани = | х 
Файл Правка Формат Файл Правка Форм Файл Правка Формат Вид Справка 


Текст это: письмен! АБ р Ри ЕАГ,СТекст это: письменная или печатная фиксация ^ 
высказывания или <! °АБ&р МиВ ЁАБ р речевого высказывания или сообщения в 

реализации; Выражен| СА... + ШАВП А..,а П] противоположность устной реализации; Выраженная и 
знаков (независимо! [и@ АБ ^ иВ ЕАГ закрепленная посредством языковых знаков 
реализации) чувств! ТАГ"`У Ми ЕР.. | (независимо от письменной или устной формы их 
т.ч. литературного [°АБ`У Еи@: „(4 реализации) чувственно воспринимаемая сторона 
коммуникации, обла, ЕАГ а УаВ ЕР....р речевого, в т.ч. литературного, произведения; 
(целостностью) и о! иВИЕРЬ"АУМи@ АБ Минимальная единица речевой коммуникации, 

Текст во втором зн! Ои@ А... `УШЕиВу Е обладающая относительным единством 

в третьем значении @и@ Уд. (+У (целостностью) и относительной автономией. 
отправной точкой в! ВПЕР..$р АВ м Текст во втором значении, являющийся одним из 
установление на ос! АЁ"^ | иВПЕР..,`У аспектов текста в третьем значении (а именно, 
отношений между ег! $`У ГиВШЁА ОУ ЕаЁ его «планом выражения»), служит отправной точкой 
Нередко под тексто! ЕРЬ$а Ъ ПИ "а всех филологических процедур. Его установление 
от отдельной морфа Р....р ПАи@ЁА `УШМна основе критического изучения истории текста и 
употребление размы! и@ЕА,,.` Еи@ЕА.. | отношений между его вариантами составляет задачу 
исследовательских | (БЕ. аУ Еи@Е *% текстологии. Нередко под текстом понимается 
Текст в третьем зн, Р..$4У и@ АГ. |также любой речевой отрезок — от отдельной 
«речевое произведе|! и@ПУАБ р Н ВЕРЬ {морфемы до литературы в целом, однако такое 
«теории целого тек Фа ТАБ а употребление размывает специфику понятия и ведет 
своих частей, или { \иВШР,, ру к нечеткости исследовательских результатов. 
обладает известной А.. `УЦЕаВЦУА.."рУ |Текст в третьем значении («сообщение», «целый 
специфической орга|й АГ"^ Е ВИТРЬЗа} текст» › «речевое произведение») — предмет 
организацией этих | Еи@: Г"ЕУ и А, особого раздела филологии, «теории целого 
внутренней структу|Р Ра@ А,,(4 ^ Втекста», или лингвистики текста. В отличие от и 








Рис. 4. Результат работы программного средства 


Эксперименты. Цель экспериментальное исследования — показать и доказать преимущества по харак- 
теристикам производительности разработанных криптосистем в сравнении с аналогом. (В рамках данной рабо- 
ты в качестве альтернативной была выбрана криптосистема с открытым ключом В $А [6].) Кроме того, необхо- 
димо было сравнить реализованный алгоритм МТКОЕпсгур! и его модификацию, чтобы отметить преимущества 
последней и привести экспериментальные доказательства. 

Для обоснования достоверности результатов использовался статистический И-критерий Манна — Уит- 
ни, который используется для сравнения двух независимых выборок по уровню какого-либо признака, изме- 
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ренного количественно. Критерий позволяет определить степень различия между выборками и является более 
мощным по сравнению с критерием Розенбаума. 

Для проведения эксперимента была использована сторонняя программная реализация криптографиче- 
ской системы В$А. В исходный код всех трех программ МТКОЕпсгурь ВА, модификации МТКОЕпсгур! был 
внедрен элемент класса Зюру’асв. Данный класс предоставляет набор методов и свойств, которые можно ис- 
пользовать для точного измерения времени, затраченного на выполнение. Таким образом, появилась возмож- 
ность фиксирования результатов затраченного времени на всех трех основных этапах: создание ключей, шиф- 
рование и расшифрование сообщения. Результаты эксперимента представлены в табл. 1. 








































































































Таблица 1 
Результаты эксперимента 
№ Время работы Время работы модификации Время робу В ЗА, вех 
МТКОЕйпсгурь, сек МТКОЕпсгурь сек 
Генерация ключей 
1 0,0239676 0,0190014 0,0964144 
2 0,0149743 0,0109994 0,1023557 
3 0,0129915 0,0099936 0,1372658 
4 0,0170503 0,0099772 0,0491137 
5 0,0139885 0,0099773 0,0869555 
6 0,0139880 0,0099782 0,0587503 
7 0,0129761 0,0099936 0,0986608 
8 0,0139931 0,0109775 0,0707417 
9 0,0139918 0,0099773 0,0595015 
10 0,0169748 0,0099777 0,0517874 
Шифрование 
1 0,0069961 0,0069770 0,0453772 
2 0,0049961 0,0059796 0,0598658 
3 0,0059954 0,0049965 0,0265461 
4 0,0049966 0,0069966 0,0402407 
5 0,0069765 0,0060854 0,0783415 
6 0,0049961 0,0069989 0,0097463 
7 0,0059795 0,0059954 0,0247876 
8 0,0069961 0,0059954 0,0272320 
9 0,0089956 0,0049798 0,0494907 
10 0,0059954 0,0049744 0,0100916 
Расшифрование 

1 0,0009843 0,0009843 1,0624476 
2 0,0010212 0,0010026 0,3692022 
3 0,0010193 0,0009974 0,7651423 
4 0,0009989 0,0009975 0,9922103 
5 0,0030012 0,0009988 0,3056423 
6 0,0010017 0,0009975 0,3757183 
7 0,0009989 0,0009988 1,2844272 
8 0,0010021 0,0009989 0,5902735 
9 0,0010026 0,0009984 0,9932694 
10 0,0009844 0,0009970 0,8740392 
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К результатам, приведенным в табл. 1, применим О-критерий Манна — Уитни. Сравним скорость гене- 
рации ключей алгоритма МТКОЕпсгур! и его модификации. 

Следуя алгоритму И-критерия Манна — Уитни, распишем вычисление критерия по шагам. 

1. Составить общий ранжированный список из обеих выборок, присвоив меньшему рангу меньшее зна- 
чение. 

2. Разделить общий ранжированный список на два, состоящих из элементов первой и второй выборок. 

3. Подсчитать сумму рангов, приходящихся на первую и вторую выборки по отдельности, как показано 
на рис. 5. 


ЕЕ И ЕЕ ЕЕ 
ПЕ ПО ЕЕ ОИ ПОЗ ООО ООС 
ПЕ ПОС ЗО ПОЗОИ ПОООСЕООО ПООЕОИ 
ОСИ ПОТЕ ОО ПО ЗИ ООО СООО ПООЕИ 
ПЕ ПОТЕ ООО ПОЗ ООО СООО ООС 
р и | о | 
О ОС ООО ПОС ПОХ ПОС 
О О О О ОС 


Рис. 5. Второй шаг вычисления И-критерия Манна — Уитни 





4. Вычислить значение О-критерия Манна — Уитни: И = 9. По таблице для уровня статистической зна- 
чимости следует определить критическое значение критерия для данных и] ип? (рис. 6). 





Рис. 6. Четвертый шаг вычисления И-критерия Манна — Уитни 


Так как значения п и и2 для всех экспериментов одинаковы, то эта таблица будет использоваться в 
каждом вычислении. 

Отсюда можно сделать вывод, что полученное эмпирическое значение И = 9 находится в зоне значимо- 
сти (рис. 7). Следовательно, наблюдается значительное различие между скоростью работы программы 
МТКОЕпсгур! и ее модификацией. Если этот факт выразить в процентном отношении, получится, что модифи- 
кация МТВОЕпсгур{ быстрее самой программы на 28 %. 


Ось значимости: 





Рис. 7. Ось значимости 


Аналогично, опуская подробности вычислений, приведем результаты И-критерия Манна — Уитни и 
процентное превосходство для остальных случаев. 

Генерация ключей 

1. Модификация МТКОЕпсгуре быстрее программы МТВОЕпсгурй на 28 %. Полученное эмпирическое 
значение И = 9 находится в зоне значимости. 
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2. Программа МТВОЕпсгур быстрее ВЗА на 80 %. Полученное эмпирическое значение И = 0 находится 
в зоне значимости. 

3. Модификация МТВОЕпсгур( быстрее ВЗА на 86 %. Полученное эмпирическое значение И = 0 нахо- 
дится в зоне значимости. 

Шифрование сообщения 

1. Полученное эмпирическое значение И = 47 находится в зоне незначимости. 

2. Программа МТВОЕпсгур быстрее ВЗА на 80 %. Полученное эмпирическое значение И = 0 находится 
в зоне значимости. 

3. Модификация МТКОЕпсгур! быстрее ВЗА на 84 %. Полученное эмпирическое значение 0 = 0 нахо- 
дится в зоне значимости. 

Расшифрование 

1. Модификация МТКОЕпсгур! быстрее программы МТКОЕпсгуре на 17 %. Полученное эмпирическое 
значение И = 24 находится в зоне неопределенности. 

2. Программа МТВОЕпсгур быстрее ВЗА на 99 %. Полученное эмпирическое значение И = 0 находится 
в зоне значимости. 

3. Модификация МТКОЕпсгур! быстрее ВЗА на 99 %. Полученное эмпирическое значение И = 0 нахо- 
дится в зоне значимости. 

Обсуждение и заключения. В рамках данной работы созданы: 

— программное обеспечение, реализующее работу криптосистемы МТВОЕпсгур; 
— программное средство, реализующее модификацию этой криптосистемы. 

Из преимуществ МТКОЕпсгур! перед аналогом — криптосистемой ВЗА можно указать более высокую 
скорость работы. Выполнение операций шифрования и расшифрования требует О(и?) операций, в отличие от 
О(?) у того же В$А. По экспериментальным данным программа МТВОЕпсгур! значительно выигрывает по ско- 
рости работы алгоритма в сравнении с ВЗА. Кроме того, отмечается небольшое увеличение стойкости при фак- 
тически такой же длине ключа. Недостаток системы — необходимость использования рекомендованных пара- 
метров. 

Что касается стойкости МТКОЕпсгурь, то после создания квантовых компьютеров будут решены задачи 
быстрой факторизации и дискретного логарифмирования [7]. В этом случае ВЗА, ОЗА и подобные им алгорит- 
мы станут бесполезными. Актуальность МТКОЕпсгур{ сохранится: он будет вполне применим и в «посткванто- 
вую» эпоху, т. к. не существует алгоритма, решающего задачу кратчайшего вектора решетки. 

Экспериментально доказано преимущество использования модификации алгоритма МТКОЕгпсгуре. Раз- 
работанное приложение на 25 % быстрее выполняет общую работу по генерации ключей, шифрованию и рас- 
шифрованию. Кроме того, оптимизируется использование внутренней памяти за счет уменьшения веса исход- 
ного файла программы и размера секретного ключа. При попытке взлома шифротекста проявляется криптогра- 
фическая стойкость и сложность использования квантовых алгоритмов. 
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Оптимизация двухстадийного режима метаногенеза на основе принципа максимума 


Понтрягина" 


С. А. Королев', Д. В. Майков?” 


' Ижевский государственный технический университет им. М. Т. Калашникова, г. Ижевск, Российская Федерация 
? Ижевский торгово-экономический техникум, г. Ижевск, Российская Федерация 
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Введение. Статья посвящена решению задачи оптимального 
управления процессом получения биогаза при непрерывном 
режиме его переработки в двух метантенках. Цели работы: 
представить математическую модель данного процесса и 
разработать алгоритм выбора оптимального управления. 
Материалы и методы. Созданная математическая модель 
описывает получение биогаза из отходов животноводства 
при последовательной переработке субстрата в двух метан- 
тенках. Рассматриваются случаи одинаковых и различных 
температурных сред (мезофильной и термофильной). Для 
данной модели сформулирована задача оптимального 
управления в виде задачи Лагранжа. Ее управляющими 
параметрами являются скорости поступления субстрата в 
метантенки. Предложен алгоритм решения данной задачи, 
основанный на чис 

ленной реализации принципа максимума Понтрягина. При 
оптимизации применялся гибридный генетический алго- 
ритм с дополнительным поиском в окрестности лучшего 
решения методом сопряженных градиентов. 

Результаты исследования. Разработана новая математиче- 
ская модель, описывающая процесс получения биогаза при 
последовательной переработке субстрата в двух метантен- 
ках. Предложен и программно реализован численный алго- 
ритм решения задачи оптимального управления.Численные 
исследования показали, что для термофильной среды ско- 
рость образования биогаза практически вдвое выше, чем 
для мезофильной. Установлено, что последовательная пе- 
реработка субстрата в двух метантенках с одинаковыми 
температурными средами позволяет вдвое увеличить ско- 
рость образования биогаза. Если температурные среды в 
метантенках различны, то в первом из них следует исполь- 
зовать мезофильную среду, а во втором — термофильную. 
При этом скорость образования биогаза несколько ниже по 
сравнению со случаем, когда в каждом из метантенков ме- 
зофильная среда, однако степень переработки субстрата 
выше на 10—15 %. 


* Работа выполнена в рамках инициативной НИР. 
“* Е-пай: з@шай.га, Маукоур@уапдех.га 
“* ТВе гезеагсВ 15 Чопе ут Ше йаше оЁе шерепдене В&). 
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Обсуждение и заключения. Полученные результаты могут 
быть использованы для расчета и конструирования биога- 
зовых установок, а также при разработке соответствующего 
программного обеспечения. 


Ключевые слова: метаногенез, биогаз, метантенк, перера- 
ботка отходов животноводства, математическая модель, 
система дифференциальных уравнений, численное реше- 
ние, методы оптимизации, оптимальное управление, прин- 
цип максимума Понтрягина. 
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Введение. При работе животноводческих предприятий, в частности, птицефабрик, свинокомплексов, 
ферм крупного рогатого скота (КРС), образуется большое количество отходов. В результате их сбраживания 
(ферментации) в специальных резервуарах (метантенках) можно получить топливный газ (биогаз) и ценные 
органические удобрения [1]. Этот процесс называется метаногенезом. Он может протекать при периодическом 
и непрерывном поступлении субстрата. В первом случае метантенк наполняется однократно и полностью опо- 
рожняется по завершении ферментации. Во втором случае одновременно и непрерывно происходят два процес- 
са: поступление субстрата и удаление его переработанной порции. 

Обычно для метаногенеза используется мезофильная среда (при температуре 25-38 °С) или термо- 
фильная (45—60 °С). Оптимальная температура для мезофильной среды — 37 °С, для термофильной — 56 °С. 
Продолжительность ферментации для этих сред равна 25 и 12 суткам соответственно. 

Экономическая эффективность производства биогаза зависит от различных факторов: вид и количество 
сырья, климатические условия [2] и т. д. Кроме того, на выработку биогаза существенно влияет скорость по- 
ступления субстрата в метантенк. Значение данного параметра зависит от объема метантенка и вида сырья. 
Чтобы найти оптимальное значение указанной величины, требуется решить задачу оптимального управления. 

Следует отметить, что при непрерывном режиме ферментации субстрат не успевает пройти полную пе- 
реработку. Для повышения выработки биогаза необходимо использовать два метантенка, чтобы субстрат по- 
следовательно перерабатывался в каждом из них. Различные аспекты такого процесса изучены в работах [3, 4], 
а его технические реализации отражены в патентах [5-7]. 

Во работах, посвященных математическому моделированию метаногенеза, например [8, 9], не пред- 
ставлен поиск оптимального управления либо находится его асимптотическое значение [10-12]. Аналитическое 
решение задачи оптимального управления получено в [13]. Однако описанная там математическая модель су- 
щественно отличается от представленной в данной статье. Предлагаемый здесь численный метод решения за- 
дачи оптимального управления применим к широкому классу моделей. 

Материалы и методы. Схема последовательной переработки субстрата представлена на рис. 1. 


Метантенк 1 Метантенк 2 


То, Чо Ч 4> 
Хь Е Х,, Г, 


Рис. 1. Схема двухстадийного режима метаногенеза 





Пусть [0 — концентрация питательных веществ в субстрате, поступающем в первый метантенк 
(кг/м°); [ и Г› — концентрация питательных веществ в субстрате для первого и второго метантенков соответ- 
ственно (кг/м’); Хр и Х› — концентрация метанообразующих бактерий в первом и втором метантен- 
ках (кг/м°); О; — объем субстрата в первом метантенке (м), О› — во втором; 40 — скорость поступления 
субстрата в первый метантенк (м”/сут), 41 = "2 и 42 = скорости выделения субстрата из первого и 


второго метантенков соответственно. 
Потерь субстрата не происходит, поэтому 
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40 = 41 = 92 =Ч4. 


Величина Чо определяется по объему отходов, количеству и производительности биогазовых устано- 








ВОК. 
Относительные скорости поступления субстрата для метантенков составляют: 
та 1 а 
те и ‚902 _ 92 
Оо @ О Ш 0 
Объемы субстратов в метантенках равны 
Ч 4 _ щ О, 
О =— и 9 =— = | 
и и> и> 


Для описания метаногенеза при последовательной переработке субстрата используется математическая 
модель, основанная на математической модели популяционной динамики метанообразующих бактерий для од- 
ностадийного режима метаногенеза [11, 14]: 


41 _| ВтаА Ни 














.Х,, 
и Ч] + [ Ы + [ 
Ч т Г Ри 
с =щ: (1-2) 
АХ. Нта2 72 Ь 
2 а 7 Ита202 ца 55 
ф а2 + [> Ь> яв [2 
4 у В>Нте272 > 
= =и2 ‘(М -12) | 
ИИ а2 + [2 
Здесь нижние индексы переменных ( Х,, [,) и параметров (Нис г, Нид г, а, В, Вг, и; ) соответствуют 
номеру метантенка Ге {1, 2}. Параметрами модели служат Ии.; И Иша; — максимально возможные относи- 
тельные скорости соответственно прироста и отмирания бактерий (сут '); В; — безразмерный коэффициент 
усвоения субстрата; а; и Ь; — эмпирические коэффициенты (м\/кг); ^ — параметр, равный нулю, если темпе- 


ратурные среды в метантенках различны, и единице, если эти среды совпадают. 
Значения параметров модели определяются в соответствии с выбранным температурным режимом ме- 
таногенеза (мезофильным ИЛИ термофильным). Если температурные среды в метантенках совпадают, то 


Итет = Ите2 = Ито › Ита1 = Ита2 = Ита ИТ. Д. 
Начальные условия: 
Х! (0) = Х, (0) = Хо, 11 (0) = 15 (0) = 1%, 
где Хо — естественная концентрация метанообразующих бактерий в исходном сырье; Гу равно концентрации 


питательных веществ в непереработанном субстрате. 

Модель (1) построена в предположении, что в метантенках поддерживается оптимальная и неизменная 
температура процесса. Если температурные среды различны, то при поступлении субстрата из одного метан- 
тенка в другой должен происходить нагрев или охлаждение до температуры во втором метантенке. 

Скорость образования биогаза (м”/сут) в -м метантенке равна 


УИ те ИХ 
и; =—— , 
а; +2; 
где у; == коэффициент, характеризующий скорость преобразования питательных веществ субстрата В биогаз 


(м. м/кг). 
Для получения задачи оптимального управления необходимо дополнить систему уравнений (1) крите- 
риальным функционалом 





Т 
&Х [5 Х. 
= | УИ 1 72Нте? 2—2 @-> тах, (2) 
0 


“1 +Ё 92 +2 
определяющим общий выход биогаза с 1 м* субстрата в первом и втором метантенках за время Г. Оптимизи- 
руемыми параметрами задачи служат относительные скорости поступления субстрата в метантенки и и и В 


системе уравнений (1). 
Системы (1)-(2) являются задачами Лагранжа. В общем случае задачи оптимального управления име- 
ется система дифференциальных уравнений вида 


(хе), (0), 0), го, т жек вая”, (3) 
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ее начальные условия: 


х (0) =хХо. 
Требуется найти оптимальное управление и (г), доставляющее максимум критериальному функциона- 
лу 
Т 
7=[ю(<(9, ч(1), 6) тах. (4) 
0 
Для численного решения задачи на промежутке [0,7] вводится разностная сетка с узлами 1 =0, 4, 
2, В, И, -. и =Т с постоянным шагом й=1;.1 —& [15]. 


Численное решение системы дифференциальных уравнений (3) осуществляется методом Рунге — Кутты 
четвертого порядка: 





Хх =Х; + +2К> +2Кз +К4), 1=0, 4—1 › 
К = м, в), 


й й 
к [ме мьи+ (5) 


й й 
К: =Ё х,+—.К, м, Е +- |, 
3 [ 2 2 "] 


Для решения задачи удобно ввести 4^-мерный полный вектор управлений Ч = (и;), 1=а4. При этом 


разностной аппроксимацией критериального функционала (4) является выражение: 
И 
1=/(0)=У 8 (хр, ц,,)-й > шах. (6) 
1=0 


Задача решается в приведенной ниже последовательности. 

1) Задается полный вектор управлений Ц. 

2) Численно решается исходная система дифференциальных уравнений (3) с помощью соотношений (5) 
и вычисляется значение критериального функционала (4) с помощью разностной аппроксимации (6). 

3) Численно интегрируется система сопряженных уравнений (в направлении «справа налево») согласно 
соотношениям: 





ОР : 
рр =рни+й- +0; рн" 8=Ъа-1, (7) 


1 


ра =0. 





В ОЕ (х,, и, й,) 
Здесь р = р(#) — двойственные переменные принципа максимума Понтрягина; С; = Эх матри- 
1 
ца Якоби, составленная для системы (3). 

4) Выполняется процесс оптимизации по вектору 1. В данной работе в качестве метода оптимизации 
использовался генетический алгоритм с вещественным кодированием и дополнительным поиском в окрестно- 
сти лучшего решения с помощью метода сопряженных градиентов. 

Для рассматриваемой задачи оптимизации метаногенеза вектор фазовых переменных равен 


х = соот (ХТ, [, Х>, 12) ‚ вектор управления и = союп (м, и> ) ‚ а матрица Якоби имеет вид: 




















Ни Ники я ата ; Ант 0 0 
а+и В+ (+) (+8) 
У Ва тет 1 
= Е. 0 0 
& И (а+д) 
Хи 0 Ново Миша Ч2Мтаа Хо, БортизЖ 
ы 2 
92 + [> [и + [> (а> + [> } (> + > } 
ь у В>Нте2.72 , Ва Ите2 Хо 
2 А 2 
а> +15 (а+ы) 
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Результаты исследования. Фиксируются теоретические значения переменных, найденных в результа- 
те численного решения системы дифференциальных уравнений математической модели метаногенеза. Учиты- 
ваются отклонения указанных значений от экспериментальных значений относительно вектора параметров [10, 
16] (табл. 1). Значения параметров модели метаногенеза оцениваются путем минимизации суммы квадратов 
этих отклонений. 

Таблица 1 


Оценки значений параметров модели метаногенеза для мезофильной / термофильной сред 






































Источник сырья 
Коэффициент 

Птицефабрики Свинокомплексы Фермы КРС 

Ито 0,482 / 0,821 0,346 / 0,783 0,297 / 0,563 

Ита 0,353 / 0,528 0,291 / 0,423 0,254 / 0,351 

а 34,781 / 43,875 7,242 / 21,653 5,013 / 8,733 

Ь 116,457 / 14,674 37,347 / 9,278 18,722 / 5,455 

В 2,344 / 3,189 1,495 / 2,084 1,413 / 1,983 

у 1,463 / 1,963 1,373 /1,907 1,299 / 1,813 








Рассматриваются следующие способы переработки субстрата. 
1. Используется единственный метантенк. 


П. Используются два метантенка с последовательной переработкой субстрата, температурные среды в 


которых совпадают. 


При этом для каждого метантенка возможны два варианта среды: мезофильная ИЛИ термофильная 



































(табл. 2). 
Оптимальные параметры метаногенеза 
ое ристиаа Источник сырья 
Птицефабрики Свинокомплексы Фермы КРС 
Т. Один метантенк с мезофильной / термофильной средой 
и”, сут" 0,149 / 0,402 0,128 / 0,341 0,112 /0,268 
и’, м’/сут 7,25 / 24,51 3,58 / 10,76 1,81 /5,42 
В: 56 / 43 51/38 52/41 
П. Два метантенка с мезофильной / термофильной средой в каждом 

и, сут ' 0,168 / 0,432 0,144 / 0,362 0,117/0,287 
и> ‚сут 0,491 / 2,448 0,412 / 1,858 0,367 / 1,424 
и, м’/сут 15,19 / 50,21 7,79 /22,19 3,85 / 10,92 
12/1, % 45 /35 38 / 33 38 / 33 














Ш. Два метантенка, в первом — мезофильная, во втором термофильная среда / 


в первом — термофильная, во втором — мезофильная среда 




















и] ‚сут ' 0,149 / 0,402 0,128 / 0,341 0,112 / 0,268 
из, сут" 0,312 / 0,051 0,222 / 0,059 0,183 / 0,051 
и, мУсут 15.35 /25,12 6,61 / 11,42 3,26 / 5,42 
2 / Ц, % 30/34 29/35 26/32 
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* * 
Здесь и —_ оптимальные значения относительной скорости поступления субстрата; иж — соответ- 


* 


ствующий ИМ среднесуточный выход биогаза; —_ —_ полнота переработки питательных веществ субстрата. 
0 


Графики скорости образования биогаза показаны на рис. 2. 


Скорость образования биогаза: 
— — в схеме с одним метантенком 
в схеме с двумя метантенками 











Скорость образования биогаза: 
— —_ в схеме с одним метантенком 
в схеме с двумя метантенками 








0 20 40 60 80 100 
+ сут 
6) 
Рис. 2. Скорость образования биогаза при использовании отходов свинокомплексов для двух метантенков 
с мезофильной (а) и термофильной (6) средой в каждом 


Ш. Используются два метантенка с последовательной переработкой субстрата, температурные среды в 
которых отличаются. 

При этом рассматриваются два варианта. 

1) В первом метантенке мезофильная среда, а во втором — термофильная. 

2) В первом метантенке термофильная среда, а во втором — мезофильная. 

Оптимальные значения параметров см. в табл. 2. 

На рис. 3 приведен график скорости образования биогаза. 
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у, Скорость образования биогаза: 
— — в схеме с одним метантенком 
м` (сут в схеме с двумя метантенками 

















0 20 40 60 80 100 
Ь сут 
а) 
и’ Скорость образования биогаза: 
— — в схеме с одним метантенком 
м` /сут в схеме с двумя метантенками 


Винни 

ИГ. 
ДТ 
Е 


0 20 40 60 80 100 
1, сут 











6) 
Рис. 3. Скорость образования биогаза при использовании отходов свинокомплексов и двух метантенков: в первом 
метантенке мезофильная, во втором — термофильная среда (а); в первом метантенке термофильная, 
во втором — мезофильная среда (6) 


Использование одного метантенка при непрерывном режиме ферментации позволяет переработать пи- 
тательные вещества субстрата на 40-45 % для мезофильной среды и на 50—55 % — для термофильной. 

При использовании двух метантенков с одинаковыми температурными средами доля переработки суб- 
страта увеличивается на 10 % для мезофильной и на 5 % для термофильной сред. При этом среднесуточная вы- 
работка биогаза возрастает практически вдвое. Скорость образования биогаза для термофильной среды практи- 
чески вдвое выше, чем для мезофильной. 

Два метантенка, эксплуатируемые в различных температурных средах, демонстрируют различные ре- 
зультаты. Так, при использовании первого варианта доля переработки субстрата примерно на 25 % выше, чем 
при использовании одного метантенка с мезофильной средой. При втором варианте доля переработки субстрата 
на 5—10 % выше, чем при использовании одного метантенка с термофильной средой. Кроме того, организация 
процесса по первому варианту обеспечивает практически двукратное увеличение скорости образования биогаза 
(если сравнивать с одним метантенком с мезофильной средой). Во втором варианте значимый прирост скорости 
не отмечается (если сравнивать с термофильной средой). Из вышесказанного следует, что использование второ- 
го варианта нецелесообразно. 

Первый вариант несколько проигрывает в скорости образования биогаза случаю использования двух 
метантенков с мезофильной средой в каждом, хотя обеспечивает более высокую (на 10-15 %) степень перера- 
ботки субстрата. Это связано с тем, что при использовании двух метантенков с мезофильной средой оптималь- 
ные значения относительных (а следовательно, и абсолютных) скоростей поступления субстрата выше, т. е. в 
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единицу времени перерабатывается больше субстрата. Скорость образования биогаза для второго варианта 
вдвое ниже, чем для случая использования двух метантенков с термофильной средой. 

Обсуждение и заключения. Разработаны математические модели, описывающие метаногенез при по- 
следовательной переработке субстрата в двух метантенках, когда температурные среды в них совпадают и от- 
личаются. Этим моделям соответствует задача оптимального управления процессом метаногенеза на основе 
принципа максимума Понтрягина. В работе представлен алгоритм ее решения. Управляющими параметрами 
служат относительные скорости поступления субстрата в метантенки. 

В результате численных исследований обнаружено, что последовательное использование двух метан- 
тенков с одинаковыми температурными средами позволяет увеличить степень переработки субстрата на 5-— 
10 %. При этом скорость образования биогаза возрастает вдвое. Степень переработки субстрата для мезофиль- 
ного режима также на 5—10 % выше, чем для термофильного. При этом скорость образования биогаза в термо- 
фильной среде практически вдвое выше, чем в мезофильной. Это связано с более высокой интенсивностью 
процесса (выше оптимальное значение относительной скорости поступления субстрата). 

Также показаны результаты эксплуатации метантенков с различными температурными средами. В пер- 
вом варианте субстрат сначала перерабатывается в мезофильной среде, а затем поступает в метантенк с термо- 
фильной средой. Во втором варианте, наоборот, субстрат из метантенка с термофильной средой поступает в 
метантенк с мезофильной средой. При первом варианте перерабатывается больше субстрата. Если сравнивать с 
применением одного метантенка с мезофильной средой, преимущество составляет около 25 %. Если сравнивать 
с использованием одного метантенка с термофильной средой — 5-10 %. Для первого варианта скорость обра- 
зования биогаза практически вдвое выше по сравнению со случаем одного метантенка с мезофильной средой. 
Во втором варианте скорость образования биогаза почти такая же, как при использовании одного метантенка с 
термофильной средой. Таким образом, применение второго варианта нецелесообразно. 

Для первого варианта скорость образования биогаза несколько ниже, чем в случае применения двух 
метантенков с мезофильной средой в каждом, но в то же время степень переработки субстрата выше на 10- 
15 %. Скорость образования биогаза во втором варианте вдвое ниже, чем при использовании двух метантенков 
с термофильной средой в каждом, что еще раз показывает неэффективность второго варианта. 
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